Papers for Friday, Jan 23 2026

We present OmniSpectra, the first native-resolution foundation model for astronomy spectra. Unlike traditional models, which are limited to fixed-length input sizes or configurations, OmniSpectra handles spectra of any length at their original size, without resampling or interpolation. Despite the large-scale spectroscopic data from diverse surveys fueling the rapid growth of astronomy, existing foundation models are limited to a fixed wavelength range and specific instruments. OmniSpectra is the first foundation model to learn simultaneously from multiple real-world spectra surveys with different configurations at a large scale. We achieve this by designing a novel architecture with adaptive patching across variable lengths, sinusoidal global wavelength encoding, local positional embeddings through depthwise convolution, and validity-aware self-attention masks. Allowing us to learn multi-scale spatial patterns while skipping attention for invalid patches. Even with a limited training example, OmniSpectra demonstrates excellent zero-shot generalization compared to methods tailored for specific tasks. This transfer learning capability makes this model the state-of-the-art across various astronomy tasks, including source classification, redshift estimation, and properties prediction for stars and galaxies. OmniSpectra reduces the need for training individual models for different tasks from scratch, establishing itself as the next-generation astronomy foundation model.

LISA, the Laser Interferometer Space Antenna, due to launch mid-2035, is a large class space mission by the European Space Agency (ESA). In partnership with NASA and ESA-member states, ESA is on track to launch what is expected to be the first space-based gravitational wave detector. By hosting detectors in space, one gains access to a lower frequency band of gravitational wave sources and with them, a plethora of new science. To maximise this scientific gain, ESA and NASA selected 20 scientists for the LISA Science Team, to carry out and/or lead necessary actions on the run up to LISA launch. We give a short overview and update of the LISA mission, some of its science objectives and related waveforms, as well as the work of the LISA Science Team as of December 2025.

A new, iterative algorithm is presented to calculate the Embedded Element Pattern (EEP) tranformation from a set of patterns computed for a uniform antenna port loading (scaled identinty matrix) to a set of those computed for a non-uniform one (arbitrary diagonal matrix). This method proves particularly useful when inverting the computations to derive the non-uniform entries of the arbitrary load, given the minimum number of EEPs necessary, which disposes of the redundancy of other matrix-based computations and leads to numerically stable impedance fault calculation. As the EEPs are envisioned to be obtained primarily through measurement, our method is also tested with the inclusion of various noise components and its convergence is evaluated, suggesting the minimum SNR and fading level of the measurement apparatus, as well as the optimal choice of reference antenna to minimise the estimation error.

We study the kinematics of condensations in 30 fields forming high-mass stars with ALMA at a high-resolution of ~0.08'' on average (~230 au). The presence of disks is important for feeding high-mass stars without feedback halting growth as their masses increase. In the search for velocity gradients resembling rotation that can reveal the presence of disks, we analyze the emission of gas tracers in 49 objects using CH$_3$OH, CH$_3$CN, and tentative detections of HNCO and cis-HCOOH. Most of the velocity distributions show velocity gradients indicative of rotation. We reveal a total of 32 disk candidates, the largest sample to date that has been uniformly analyzed at a few hundred au scales in the high-mass regime. Their position-velocity maps are generally asymmetric with one side brighter than the opposite. We successfully fit a power law to the position-velocity maps of the disk candidates and find indices between -0.5 (Keplerian rotation) and -1 (rotation under specific angular momentum conservation) with a median of -0.7. Under Keplerian rotation assumption, we estimate central masses, uncorrected for inclination, ranging between 7 to 45 M$_\odot$. Excluding outliers, the disk candidates are relatively more compact (<200 au) and less massive (<5 M$_\odot$) than previous results at coarser angular resolution. We calculate an average Toomre-$Q$ parameter and find that most are gravitationally unstable (median of 0.5). We conclude that these observations offer the first opportunity to separate the disk and envelope components of hot cores on a statistically significant sample, and confirm that anisotropic collapse plays an role in feeding high-mass (proto)stars.

Stellar streams, the debris of tidally disrupted satellites, trace their host's gravitational potential and thus probe dark matter halo structure. While six-dimensional phase-space data of Galactic streams enable precise dark matter halo modelling in the Milky Way, streams around external galaxies are typically available only as low surface brightness features without kinematics (i.e. two-dimensional photometric data), providing only weak constraints when considered individually. We present a hierarchical Bayesian framework that infers the population distribution of halo flattening using only projected stream tracks. Streams are forward-modelled in StreaMAX, a new JAX-accelerated particle-spray package that achieves orders of magnitude faster stream generation when compared to traditional methods. For each stream we fit an axisymmetric dark matter halo model and obtain a posterior on the flattening. These posteriors are then combined through hierarchical reweighting to constrain the population distribution. Using mock data, we show that individual fits recover the correct flattening with modest precision and exhibit projection-induced multi-modalities. Nevertheless, aggregating these fits yields accurate and confident constraints on the underlying population distribution of dark matter halo morphologies, clearly distinguishing between oblate, spherical, and prolate populations. The total computational cost scales linearly with sample size. Our results demonstrate that ensembles of purely photometric streams carry sufficient information to constrain dark matter halo shapes in external galaxies at the population level. With the forthcoming samples from Euclid and Rubin/LSST, this approach offers a practical path to population-level inferences of halo morphology without any kinematic measurements.

Radio observations offer a dust-unobscured view of galaxy star formation via the radio continuum-star formation rate (RC--SFR) relation. Emerging evidence of a stellar mass dependence in the RC--SFR relation raises the broader question of how other galaxy properties may influence this relation. In this work, we study the dependence of the global RC--SFR relation on galaxy properties in local ($z\,\leq$\,0.3) star-forming galaxies (SFGs) using the second data release of the LOFAR Two-Metre Sky Survey (LoTSS-DR2). Employing a non-parametric decision-tree regression algorithm, we identify the most important galaxy properties for estimating the radio luminosity using a sample of 18,828 emission-line-classified SFGs based on spectroscopic data from the SDSS-DR8. Along with the spectroscopically obtained SFRs and stellar mass values, we also use SFRs and stellar masses derived using photometric SED-fitting from the \textit{GALEX}--SDSS--\textit{WISE} Legacy Catalogue (GSWLC) for the same sample. We find that a galaxy's SFR is most important for predicting the radio luminosity, followed by the stellar mass, at $>5\sigma$ significance. Complementing the LoTSS catalogue 150\,MHz flux densities with aperture photometry for the rest of the emission-line classified sample (35,099 galaxies in total), we obtain a new calibration of the RC--SFR relation, which does not change significantly whether we use spectroscopic or photometrically derived SFRs and stellar masses, despite the fact that the methods probe star formation on different characteristic timescales. Our results highlight the utility of decision-tree algorithms for handling censored radio-selected galaxy samples, which will be useful for future spectroscopic surveys of radio sources.

Cosmic voids have increasingly emerged as a powerful cosmological probe. However, their large spatial extent and intrinsically underdense environments make their identification highly sensitive to shot noise, redshift-space distortions (RSD), and observational systematics, particularly for topological and density-based void definitions. We introduce the Back-In-Time Void Finder (BitVF), a novel dynamical and physically motivated algorithm that identifies cosmic voids as regions of negative divergence of the Lagrangian displacement field reconstructed from the present-day tracer distribution. The reconstruction relies on an optimized discrete optimal transport algorithm that recovers the backward-in-time dynamics of tracers, naturally accounting for tracer bias without relying on cosmological assumptions. We validate BitVF against the widely used topological void finder REVOLVER using high-resolution N-body simulations, showing that it produces void catalogs with smoother and more physically motivated density profiles, as well as abundances that are more stable under tracer subsampling and shot noise. We further apply it to realistic DESI-like mock light-cone galaxy catalogs, demonstrating that it intrinsically mitigates redshift-space systematic effects, preserving real-space void size functions more faithfully than topological methods. Modeling RSD, the reconstruction can be combined with a fiducial cosmology and an assumed tracer bias within a bias-corrected Kaiser framework, yielding reconstructed-space void catalogs consistent with real-space statistics across redshift. Its performance is characterized as a function of the main internal parameters, showing an optimal balance between accuracy, computational efficiency, and applicability to stage IV galaxy surveys. BitVF will be publicly released within the CosmoBolognaLib.

Simulations predict that supermassive black hole binaries (SMBHBs) will exhibit periodic brightness variations that may exceed the stochastic variability intrinsic to active galactic nuclei (AGN). In this paper, we simulate SMBHBs with damped random walk (DRW) AGN variability and an added sinusoidal signal from the orbital motion, and test three methods -- the Generalized Lomb Scargle Periodogram (GLSP), the nested Bayesian sampler (NBS), and the Weighted Wavelet Z-Transform (WWZ) -- to determine which is best at recovering the periodicity. Our simulated light curves follow the properties of the Catalina Real-Time Transient Survey (CRTS), Legacy Survey of Space and Time (LSST), and Zwicky Transient Facility (ZTF) to best inform current and future SMBHB searches. We map a broad range of parameter space and identify which DRW-only light curves best mimic periodicity and pass each method's model selection. The NBS performs best at detecting periodicity and filtering out DRW-only light curves. Combined candidate selection with both the NBS and GLSP significantly reduces false positive rates with marginal impact to true positive rates. With this joint model selection pipeline, we find the lowest false positive rates in ZTF-like simulations and the highest detection rates in LSST-like simulations. Using a modified computation of the False Alarm Probability (FAP) with GLSP, we efficiently triage LSST AGN light curves (~10^7 light curves in ~10-30 hours) and achieve true- and false- positive rates of ~40% and ~0.5%, respectively.

Recent work has shown that optical-radio position offsets and radio position variability are inversely correlated with the photometric variability of active galactic nuclei (AGN). A key prediction of these findings is that a reference frame constructed using highly photometrically variable AGN should be more stable than a frame that does not account for variability and that variability can be used to optimally weight all sources in order to maximize frame stability. Using ICRF3 matched to Gaia DR3, we employed a bootstrap method to estimate the multi-epoch stability of frames constructed using AGN selected at varying levels of photometric variability. We fit vector spherical harmonics to the coordinate differences between the three ICRF3 frames (S/X, K, and X/Ka) and Gaia and quantified the statistical dispersion as a function of blazar-like (high variability), quasar-like (low variability), and intermediate-variability class. An S/X reference frame constructed using blazars exceeds the stability of a frame constructed with quasars by a factor of 6 and is twice as stable as the ICRF3 defining sources. At K and X/Ka, a blazar-based frame matches or exceeds the stability of the defining sources by a factor of 1.4 in the case of X/Ka and exceeds the stability of a frame based on quasars by over a factor of 2 in both cases. The smaller improvement at K and X/Ka is likely because sources selected at higher frequency are more likely to be blazars. We derived a variability-based astrometric covariance scaling method that results in factor of 2 reduction in frame distortions and instabilities between S/X and Gaia, with a mild improvement for K but no difference for X/Ka, which is dominated by known distortions. Our results confirm the prediction that an optimal weighting of the link between the optical and radio celestial reference frames is enabled by accounting for photometric variability.

The nature of radio outflows in radio-quiet AGN remains poorly understood. In this study, we present kpc-scale polarization observations of the Seyfert galaxy NGC\,4151 using the Karl G. Jansky Very Large Array (VLA) in B-array at 3 and 10 GHz. We find that the inferred magnetic (B-) field structures show a stratification: a `spine-sheath'-like structure, with fields perpendicular to the jet direction in the `spine' and parallel in the `sheath', is observed in the higher resolution (0.5 arcsec) image at 10 GHz. In addition, a `wind'-like component with B-fields perpendicular to the radio outflow is observed in the 3 GHz image (resolution 2 arcsec); this feature is prominent along the `receding' (eastern) jet direction. Rotation measure (RM) ranges from $-230$ to 250 rad m$^{-2}$ over the polarized regions, indicating a low-electron-density ($10^{-2}-10^{-3}$ cm$^{-3}$) tenuous medium surrounding the source causing Faraday rotation. A {tentative} RM gradient of $+75$ to $-25$ rad m$^{-2}$ is observed transverse to the northern `wind' component, while a similar gradient with opposite sign is seen across the southern `wind' component, suggestive of a helical magnetic field threading the outflow. Based on an analysis of the available radio and X-ray data, we conclude that the stratified radio outflow in NGC 4151 is magnetically-driven. The bi-conical radio `wind' is found to be massive ($1050-3200 M_\odot$) with a high mass outflow rate ($0.01-0.03$ M$_\odot$ yr$^{-1}$) but low in kinetic power ($<0.01$% of L$_{\rm{bol}}$), making it less impactful for galactic-scale feedback. Our study suggests that radio-quiet AGN may also host magnetically dominant jets and winds, even while their jets are smaller and weaker compared to radio-loud AGN.

We present serendipitous simultaneous radio and X-ray light curves of the dwarf Seyfert galaxy UGC 6728 spanning 5 years. The X-ray light curve exhibits a flaring period, followed by a gradual rise and decline. Throughout these events, the X-ray hardness ratio and spectrum do not change significantly. The radio flux is constant, as far as can be determined from its sparse sampling, until the end of the X-ray flare, then decreases by a factor of two by the midpoint of the gradual X-ray rise before returning to baseline at the end of the X-ray decline. We interpret this behavior in light of a similar event recently reported in NGC 2992, in which there is a temporary obscuration of the radio source by a blob of plasma ejected by a magnetic reconnection in the accretion disk. The energetics of the X-ray flare are consistent with those expected from magnetic disk activity. As in NGC 2992, the X-ray spectrum does not evolve during the obscuration event. We also discuss the possibility that the observed phenomena are due to normal AGN coronal flaring and variability, which is plausible but unlikely given the lack of spectral variation.

Compact galaxy groups are ideal laboratories for studying the effects of interactions between AGN and multiple nearby galaxies. Recent JWST observations of the nearby compact group Stephan's Quintet highlight tidal flows between the interacting galaxies as well as outflows from the active galaxy NGC 7319. To study the kinematics on a large scale throughout the group, we obtained spatially-resolved long-slit spectra of Stephan's Quintet at multiple slit positions with Apache Point Observatory's Kitt Peak Ohio State Multi-Object Spectrograph. We fit multiple Gaussians to the H$\alpha$ $\lambda$6563 Å and [N II] $\lambda\lambda$6548, 6583 Å emission lines to isolate the different kinematic components. We used the kinematics to develop the first biconical outflow model of the narrow-line region of NGC 7319. Using a combination of galactic rotation models, biconical outflow models, and kinematic maps of the ionized gas, we disentangled the outflows, rotation, and tidal flows in the group. We found outflow radial velocities up to 550 km s$^{-1}$ peaking at 2.6 kpc from the central supermassive black hole, and a transition from AGN-powered outflows to gravitationally-powered tidal flows at a projected distance between 2.4 -- 6.3 kpc. We performed a line ratio analysis and determined the gas shows Seyfert-like ionization out to 6.3 kpc (projected), which supports our finding that gas outside this radius is predominantly powered by tidal flows. Our separation of kinematic components in Stephan's Quintet will enable future studies of the physical conditions and dynamical forces in the ionized gas to better quantify the feeding and feedback processes of AGN in compact groups.

[Abbreviated] Context. Solar twins, stars whose stellar parameters (Teff, log g, and [M/H]) are very close to the Solar ones, offer a unique opportunity to investigate Galactic archaeology with very high accuracy and precision. However, most previous catalogs of Solar twins contain only a small number of objects (typically a few tens), and their selection functions are poorly characterized. Aims. We aim at building a large catalog of Solar twins from Gaia DR3 GSP-Spec, providing model-driven, rather than data-driven, stellar parameters including ages, together with a well-characterized selection function. Methods. Using stellar parameters from the Gaia DR3 GSP-Spec catalog, we selected Solar-twin candidates whose parameters lie within +- 200 K in Teff, +- 0.2 in log g, and +- 0.1 dex in [M/H] of the Solar values. Candidates unlikely to be genuine Solar twins were removed using Gaia flags and photometric constraints. We determined accurate ages for individual twins with a Bayesian isochrone-projection method, considering three combinations of parameters: Teff, [M/H], and either log g, M_G, or M_Ks. We also constructed a mock catalog to characterize the selection function. Results. Our final GSP-Spec Solar-twin catalog contains 6,594 stars. The mock catalog consisting of 75,588 artificial twins well reproduces the main characteristics of the observed catalog, especially for ages determined with M_G or M_Ks. To demonstrate the usefulness of our catalog, we compared chemical abundances [X/Fe] with age. We statistically confirmed the age--[X/Fe] relations for several species (e.g., Al, Si, Ca, and Y), demonstrating that trends previously identified in small but very high-precision samples persist in a much larger, independent sample. Conclusions. Our study bridges small high-precision Solar-twin samples and large data-driven ones, enabling demographic studies of Solar twins.

Accreting systems can launch powerful outflows which interact with the surrounding medium. We combine new radio observations of the accreting neutron star X-ray binary (XRB) Circinus X-1 (Cir X-1) with archival radio observations going back 24 years. The $\sim3$ pc scale wide-angle radio and X-ray emitting caps found around Cir X-1 are identified as synchrotron emitting shocks with significant proper motion and morphological evolution on decade timescales. Proper motion measurements of the shocks reveal they are mildly relativistic and decelerating, with apparent velocity of $0.14c\pm0.03c$ at a propagation distance of 2 pc. We demonstrate that these shocks are likely powered by a hidden relativistic ($\gtrsim0.3c$) wide-angle conical outflow launched in $1972\pm3$, in stark contrast to known structures around other XRBs formed by collimated jets over 1000s of years. The minimum time-averaged power of the outflow required to produce the observed synchrotron emission is $\sim0.1L_\text{Edd}$, while the time-averaged power required for the kinetic energy of the shocks is $\sim40 \left(\frac{n}{10^{-2} \text{cm}^{-3}}\right)L_\text{Edd}$, where $n$ is the average ambient medium number density. This reveals the outflow powering the shocks is likely significantly super-Eddington. We measure significant linear polarisation up to $52\pm6\%$ in the shocks demonstrating the presence of an ordered magnetic field of strength $\sim200~\mu\text{G}$. We show that the shocks are potential PeVatrons, capable of accelerating electrons to $\sim0.7~\text{PeV}$ and protons to $\sim20~\text{PeV}$, and we estimate the injection and energetic efficiencies of electron acceleration in the shocks. Finally, we predict that next generation gamma-ray facilities may be able to detect hadronic signatures from the shocks.

Establishing ages for young clusters is key for properly tracking the star formation history of a region. In this paper we investigate a new approach to estimating ages for young populations, based on the well-founded assumption that the initial mass function is the same throughout a star forming cloud. We trial this method for six young clusters in the Perseus star forming region. For all six clusters, we construct new member samples in a homogeneous way using Gaia DR3. We estimate masses by comparing 2MASS photometry to theoretical isochrones, including Monte Carlo simulations to propagate the errors. We compare the mass distributions of the clusters for a range of plausible ages, looking for a combination of ages that results in indistinguishable mass distributions across the region. We find the best fit for ages of 1 Myr for NGC1333+Autochthe, 2 Myr for IC348, 2-3 Myr for Heleus, 3-4 Myr for Mestor, 4-5 Myr for Electryon+Cynurus, and 5-8 Myr for Alcaeus. All other combinations of ages are ruled out by this criterion. The established age sequence is consistent with the relative ages inferred from disc fractions, and broadly aligns with the age sequence determined in previous studies using isochrone fitting. We suggest that this approach can be a useful complement and cross-check to established methods to estimate ages in young populations.

The dynamics of the solid crust + magnetic field lines of pulsars is a much debated issue, and remains unsettled after 50 years. Some pieces of evidence have emerged to complete and confirm theoretical calculations and expectations. We discuss in the present work an interpretation of the behavior of the ''Crab Twin'' pulsar PSR B0540-69 in terms of the evolution of the magnetic field/quakes, connecting the behavior of the braking index with the underlying platelet drift and sudden discontinuous rearrangement (fast-slip) and long-term ones (slow-slip events), suggested by analogy with existing theoretical picture observed in the Earth's crust. The relationship of this scenario with permanent torque-changing glitches seen in the Crab and other young pulsars, and a set of similar events in the same object and others is addressed. We conclude that this physical approach is in principle consistent with all these sudden events, and point out future work to clarify the whole picture.

We present an upgraded detailed numerical calculations of supernova (SN) interactions with significantly aspherical circumstellar matter (CSM), primarily formed as a disc or bipolar lobes. The circumstellar disc can arise as a result of, for example, mass transfer in a binary, while bipolar lobes can be the result of a violent pre-explosive ejection of matter, similar to the iconic cases of luminous blue variable stars (LBVs). We numerically simulate the radiation-hydrodynamic (RHD) behaviour of interaction processes using a 2D cylindrical version of the RHD code CASTRO. We then calculate light curves, spectral patterns, and polarisation profiles, all up to a relatively long time of two years after an SN shock breakout and from different directions, using the multidimensional Monte Carlo radiation transfer (MC-RT) codes SEDONA and SIROCCO. We calculated a total of five models for the two aforementioned configurations of the surrounding CSM, for stratified density levels, comparing the simulated hydrodynamic behaviour and differences in their observable properties. RHD models exhibit similar behaviour to previous adiabatic models, but with a significantly slower expansion velocity. The calculated light curves show a relatively smooth evolution in SN-disc interaction, and declines and brightening in SN-lobes interaction. Comparing models with real events with a presumed similar physical process provides guidance for selecting a more accurate CSM configuration when simulating real situations.

Observational studies have made substantial progress in characterizing quenching as a function of stellar mass and environment, but they are often limited in their ability to constrain quenching timescales and to determine the dominant environmental process responsible for the shutting down of star formation. To address this, we combine recent Sloan Digital Sky Survey (SDSS) observations with the SHARK v2.0 semi-analytic model to study the quenching of satellite galaxies in groups and clusters. We generate mock SDSS-like observations to calibrate the hot halo and cold interstellar medium (ISM) gas stripping prescriptions against observed satellite quenched fractions, finding that the previously adopted stripping prescriptions in SHARK v2.0 are too aggressive and overestimate the quenched fraction of satellite galaxies. Reducing the efficiency of both hot and cold gas stripping yields excellent agreement with observations for low- and intermediate-mass satellite galaxies. We use the calibrated model to investigate quenching timescales and find that satellites quench more quickly in clusters compared to groups, with timescales that generally decrease with increasing stellar mass. The long (>2 Gyr) timescales we measure favour hot halo gas removal as the dominant driver of satellite quenching.

Magnetars are now known to be among the most strongly polarized celestial sources in X-rays. Here we report on the $500\,\mathrm{ks}$ observation of the magnetar 1E 1547.0-5408 performed by the Imaging X-ray Polarimetry Explorer (IXPE) in March 2025. The IXPE spectrum is well reproduced by a single thermal component with blackbody temperature $kT_\mathrm{BB}\sim 0.67\,\mathrm{keV}$ and emission radius $R_\mathrm{BB}\sim 1.2\,\mathrm{km}$. The source exhibits a high linear polarization degree in the $2$--$6\,\mathrm{keV}$ band ($\mathrm{PD}=47.7\pm2.9\%$) with polarization angle $\mathrm{PA}=75^\circ.8 \pm 1.^\circ8$, measured West of celestial North. While $\mathrm{PA}$ does not appear to vary with energy, there is some evidence (at the $1\sigma$ confidence level) of a minimum in $\mathrm{PD}$ between $3$ and $4\,\mathrm{keV}$, compatible with what is expected by partial mode conversion at the vacuum resonance in a magnetized atmosphere. Phase-resolved spectral and polarimetric analyses reveal that X-ray thermal radiation likely originates from a single, fairly small hot spot with a non-uniform temperature distribution. Fitting the phase-dependent $\mathrm{PA}$ measured by IXPE with a rotating vector model (RVM) constrains the source geometry and indicates that both the dipole axis and line-of-sight are misaligned with respect to the spin axis. Under these conditions, the high polarization of the source cannot be regarded as compelling evidence for the presence of vacuum birefringence in the star magnetosphere. Nevertheless, the fact that the RVM successfully reproduces the modulation of the X-ray polarization angle and the behavior of $\mathrm{PD}$ with the energy hint once more to the presence of QED effects in magnetars.

We present a multi-epoch test of f(T) gravity with nonminimal torsion-matter coupling, combining early- and late-Universe observations. At the MeV scale, Big-Bang Nucleosynthesis constrains the fractional variation of the weak freeze-out temperature, |{\delta}{\tau}_f/{\tau}_f|, thereby mapping light-element abundances into limits on deviations from the standard expansion history. At low redshift, we confront the model with type Ia supernovae, baryon acoustic oscillations, and cosmic-chronometer data, which respectively probe distances, the late-time standard ruler, and the Hubble rate. Independent analyses highlight the complementary roles of each dataset, while a joint SNe Ia + BAO + CC fit breaks degeneracies and yields the tightest combined bounds. As an illustration, we examine two representative torsion-modified gravity scenarios: BBN strongly limits large departures from standard cosmology, whereas late-time probes remain compatible with a near-{\Lambda}CDM background. This unified approach demonstrates the power of linking early-Universe nuclear physics with precision cosmological observables in assessing torsional extensions of gravity.

Bouncing cosmologies offer an alternative to inflation by resolving the initial singularity through a contracting phase followed by a bounce into expansion. In many such models, the contracting phase is dominated by a single matter component, typically pressureless dust, which leads to an almost scale-invariant spectrum of scalar cosmological perturbations with a slight blue tilt, so that generating the observed red-tilted spectrum within this framework was challenging. In this work, we consider a more realistic scenario in which the contracting phase includes both matter and radiation, as required on physical grounds. We show that the presence of radiation can naturally induce a red tilt in the spectrum of curvature perturbations seeded by quantum vacuum fluctuations in the remote past of the contraction. Since the perturbations of the two fluids are coupled via gravity, vacuum initial conditions must be carefully defined. We demonstrate that, without fine-tuning, the resulting entropy perturbations are subdominant with respect to curvature perturbations. This suggests that a minimal two-component bounce model, involving only ordinary matter and radiation, can connect to the standard expanding cosmology with observationally viable initial conditions.

We present a comprehensive catalog of 251 potential counterparts for 250 unassociated hard X-ray sources detected in the Swift Burst Alert Telescope (BAT) 150-month hard X-ray survey. Over 150 months of observation, BAT has detected 2339 sources in the 15-150 keV energy range. Among these, 344 do not have a previously identified low-energy counterpart. Our study focuses on the analysis of soft X-ray observations at energies below 10 keV, spatially overlapping with these new Swift-BAT hard X-ray sources. Such observations were taken with Chandra, Swift-XRT, eROSITA, and XMM-Newton. Within the sample of 251 potential counterparts, 94 (37 percent) are identified as active galactic nuclei and 58 (23 percent) as galaxies. The remaining 99 sources (40 percent) include pulsars, cataclysmic variables, and unclassified soft X-ray counterparts in the 0.5-10 keV band. Redshift information is available for 139 out of the 251 sources, and its distribution is in close agreement with the redshift distribution of previous BAT catalogs. We also present the results of a small optical spectroscopy campaign of 9 out of 58 galaxies. The majority of these are classified as Seyfert 2 galaxies at redshifts slightly larger than the median of the BAT AGN sample.

We discuss a potential survey to map dark matter and the emergence of large-scale structure to redshift z ~ 1.5 (baseline) or z~3.5 (with near-IR extension) using a massively multiplexed spectrograph on a 10m-class telescope, such as the proposed Widefield Spectroscopic Telescope.

Next-generation surveys such as the Euclid survey, the Legacy Survey of Space and Time (LSST), and the China Space Station Telescope (CSST) survey are expected to discover ~10^5 galaxy-galaxy scale strong gravitational lenses. This motivates the development of scalable and robust lens modeling approaches that can efficiently and reliably learn from wide-field survey datasets before high-resolution follow-up. We design a scalable, Bayesian, Lenstronomy-based pipeline and apply it to a sample of sixteen lens candidates observed with the Parallel Imager for Southern Cosmology Observations (PISCO) on the Magellan telescope. PISCO provides four-band imaging (z, i, r, g) with colours, depth and seeing conditions comparable to LSST. To fully exploit the constraining power of this dataset, our pipeline performs simultaneous multi-band modeling, using a common mass profile across all four bands while allowing independent light profiles in each. This approach leverages color information to provide joint constraints on the lens mass and yields reduced uncertainties compared to single-band analyses. Fifteen out of sixteen PISCO lens candidates are successfully recovered with interpretable lensing configurations, including DESJ0533-2536, the first reported hyperbolic-umbilic galaxy-galaxy scale strong lensing candidate. We further assess how much model complexity can be reliably constrained given the resolution and seeing of PISCO-like data. Overall, our results demonstrate that scalable, multi-band lens modeling of ground-based data can extract meaningful constraints on mass and source morphology, providing a practical pathway to maximize the scientific return from large samples in upcoming surveys.

Post-common envelope binaries (PCEBs) consisting of a white dwarf (WD) plus a main-sequence (MS) star can constrain current prescriptions of common envelope evolution (CEE) and calibrate theoretical models of binary formation and evolution. Most PCEBs studied to date have typical orbital periods of hours to a few days and can be well explained by assuming inefficient CEE to expel the envelope. However, there are currently several systems with relatively wide orbital periods ($>$18 days). To explain these wide PCEBs, additional sources of energy have been suggested to be taken into account. Here, we present the discovery and observational characterization of a compact object ($M\,\geq\,0.58\,\rm M_{\odot}$) with a K-type star companion in the solar neighborhood ($d\sim 112$ pc) and an orbital period of $P_{\rm orb}\sim 14$ days. The compact object binary is likely to be a system consisting of a WD and a barium dwarf. Such a system with an orbital period within the gap between tight and wide binaries provides a test of whether additional energy sources are required to explain its formation. Using binary evolution models, we investigate the evolutionary history of this wide PCEB system and find that the observed properties of this source can be explained without invoking any extra energy source.

A couple of milliseconds after the merger of a binary system of neutron stars can play a fundamental role in amplifying the comparatively low initial magnetic fields into magnetar strengths. The basic mechanism responsible for this amplification is the Kelvin-Helmholtz instability (KHI) and we here report the first systematic study of the impact of rapid rotation on the KHI-amplification process exploiting general-relativistic magnetohydrodynamic simulations at very high-resolutions of $35\,{\rm m}$. Concentrating on four different spinning configurations, we find that aligned, anti-aligned, and mixed (aligned/anti-aligned) spin configurations lead to markedly different growth rates of the electromagnetic (EM) energy, field topologies, and vortex properties when compared to the irrotational case. These differences arise from intrinsic variations in the system dynamics, such as tidal deformation, collision strength, and contact surface area, with the anti-aligned configuration producing the largest vorticity and growth in EM energy. Importantly, while different spin configurations lead to significantly different initial growth rates of the poloidal/toroidal components, all systems converge to a specific topological partition. Our simulations are confined to a short window in time, but the different EM energies produced as a result of spin will imprint the EM emission at merger and provide information on the spinning state at merger.

In this paper, we review the extended halo material and the circumgalactic medium (CGM), including both dust and gas, and discuss promising science cases that could be realized using the KASI Deep Rolling Imaging Fast Telescope (K-DRIFT). Scattered starlight from cirrus clouds in our Galaxy poses one of the major challenges to studying the low surface brightness features of extragalactic sources. Therefore, it is essential to investigate how to discriminate extragalactic sources from the cirrus cloud features. At the same time, interstellar dust clouds themselves are fundamental to understanding dust properties and the interstellar radiation field, both of which are essential for studies of chemical evolution and star formation in our Galaxy. Measuring the reddening of background sources, such as quasars, with K-DRIFT, which benefits from its broad field of view and accurate background subtraction, allows for effective detection of extended dust in galactic halos, the CGM, and intracluster space. Observations of the H-alpha emission lines can be used to identify signatures of star formation activity within galaxies, as well as the environmental effects acting on them. Galactic winds driven by active galactic nuclei and starbursts can be traced through H-alpha emission. Strong ram pressure stripping effectively removes the interstellar medium (ISM) from galaxies. The stripped ISM becomes ionized or dissociated through mixing with the hot intracluster medium (ICM), forming H-alpha tails. The surface brightness of these H-alpha tails correlates not only with the presence of star formation in the tails but also the mixing stage of the stripped ISM and ICM. The H-alpha survey with K-DRIFT will enable the investigation of the evolutionary stages of ram pressure stripped galaxies in cluster environments, as well as the multiphase gas reservoir around galaxies and in the CGM.

For more than half a century, the puzzling W-type phenomenon in contact binaries has challenged astrophysicists. In these systems, the less massive component exhibits a higher surface temperature than its more massive companion, which is a reversal of the typical A-type configuration, where the more massive star is hotter. This counterintuitive temperature inversion defies the basic stellar physics and still lacks a widely accepted explanation. In this study, we assembled a sample of over 3,000 extensively observed contact binaries and derived their complete set of physical parameters. Our statistical analysis revealed a strong positive correlation between the occurrence of W-type contact binaries and the intensity and frequency of magnetic activities. This result strongly supports the hypothesis that magnetic activities are the primary driver of the W-type phenomenon and offers a compelling explanation for the observed transitions between the W-type and A-type.

Massive stars ($M>8\mathrm{M_\odot}$) emit neutrinos known as pre-supernova (pre-SN) neutrinos through thermal and nuclear interactions for cooling the stellar core during the final stage of stellar evolution. Real-time monitoring of their pre-SN neutrino interaction rate offers a crucial opportunity to issue an early warning to a core-collapse supernova. Some neutrino detectors, including KamLAND and Super-Kamiokande already operate pre-SN alarm systems based on a statistically significant excess of the observed event rate over the expected background. To improve alarm sensitivity, we propose an alarm method which incorporates the time evolution of the observed pre-SN neutrino event rate. The method uses a log likelihood ratio test that references multiple theoretical stellar-evolution models and treats the core collapse time as a nuisance parameter to be profiled over. The performance of the proposed method was evaluated using simulated data for the KamLAND, Super-Kamiokande with dissolved Gadolinium (SK-Gd) and their combined analysis. The results demonstrate a significant improvement in the warning time compared to the conventional rate-only method, while maintaining the same false alarm rate.

Transmission spectroscopy serves as a valuable tool for probing atmospheric absorption features in the terminator regions of exoplanets. Stellar surface heterogeneity can introduce wavelength-dependent contamination that complicates the interpretation of planetary spectra. We aim to investigate the atmosphere of the warm sub-Saturn NGTS-5b through optical transmission spectroscopy. Two transits were observed with the low-resolution Optical System for Imaging and low-Intermediate-Resolution Integrated Spectroscopy (OSIRIS) on the 10.4 m Gran Telescopio Canarias (GTC). Chromatic transit light curves were modeled to derive optical transmission spectra and multiple Bayesian spectral retrievals were performed to characterize the atmospheric properties. Model comparisons provide strong evidence for contamination from unocculted stellar spots. A joint retrieval of the transmission spectra, assuming equilibrium chemistry, indicates a relatively clear atmosphere with a sub-solar C/O ratio of $<$0.22 (90% upper limit) and a low metallicity of $0.10^{+0.34}_{-0.05} \times$ solar. Retrievals assuming free chemistry yield strong evidence for the presence of $\rm H_2O$, with its abundance constrained to $\log X_{\mathrm{H_2O}} = -0.79^{+0.14}_{-0.17}$. However, the abundances of other species remain unconstrained due to the limited optical wavelength coverage. The discrepancies between the two NGTS-5b transit spectra can be attributed to varying levels of stellar contamination. NGTS-5b thus appears to host a relatively clear, water-rich atmosphere, pending confirmation from additional observations of molecular bands in the infrared.

T Coronae Borealis (T CrB) is one of the eleven known recurrent novae in our Galaxy. It was observed in outburst in 1866 and 1946, with additional likely eruptions recorded in 1217 and 1787. Given its predicted recurrence period of approximately 80 yr, the next outburst is anticipated to occur imminently, thus motivating a thorough examination of the main characteristics of this system. We present new hydrodynamic models of the explosion of T CrB for different combinations of parameters (i.e., the mass, composition, and initial luminosity of the white dwarf, the metallicity of the accreted matter, and the mass-transfer rate). We show that mass-accretion rates between 10-8 - 10-7 Msun yr-1 are required to trigger an outburst after 80 yr of accretion of solar-composition material onto white dwarfs with masses about 1.30 - 1.38 Msun. For lower white dwarf luminosities, less massive white dwarfs, or reduced metallicity in the accreted material, higher mass-accretion rates are required to drive an explosion within this timescale. A decrease in metallicity or initial white dwarf luminosity leads to higher accumulated masses and ignition pressures, resulting in more violent outbursts. These outbursts exhibit higher peak temperatures, higher ejected masses, and greater kinetic energies. Models computed for different white dwarf masses but identical initial luminosities reveal significant differences in the elemental abundances of a wide range of species, including Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca, and Sc. These compositional differences offer a potential diagnostic tool for constraining the parameter space and discriminating between the various T CrB models reported in this study.

Recent high-resolution imaging observations of strong lens systems reveal dense low-mass perturbers. We propose a soliton core, whose central density is boosted by a supermassive black hole (SMBH), in the fuzzy dark matter (FDM) model as an efficient perturber in strong gravitational lensing. The higher central density makes it less efficient in the tidal mass loss, and leads to the higher impact in gravitational lensing. We show that the mass profile of a $\sim 10^6M_\odot$ perturber in JVAS B1938+666, which does not resemble any known astronomical object, can be wel explained by a soliton core in the FDM model with the mass of $4\times 10^{-21}$eV hosting an SMBH with the mass of $4\times 10^5M_\odot$. The high mass of the SMBH may be explained by several scenarios that predcit heavy SMBH seeds such as the direct collapse black hole formation and primordial black holes.

Most efforts in the modeling of recurrent novae have centered on the initial phases of the explosion and ejection, overlooking the subsequent interaction of the ejecta, first with the accretion disk orbiting the white dwarf and ultimately with the secondary star. To address this gap, a series of 3D smoothed-particle hydrodynamics simulations was conducted. These simulations explored the dynamic interactions between the nova ejecta, accretion disk, and stellar companion within the framework of the recurrent nova system U Sco. Notably, the simulations incorporate rotation around the system's center of mass. The primary goal of these simulations was to qualitatively examine the impact of various model parameters, including ejecta mass, velocity, and density, as well as the mass and geometry of the accretion disk. Simulations reveal complete disruption and sweeping of the accretion disk orbiting the white dwarf star for models with flared disks and Mejecta/Mdisk larger than 1. In contrast, V-shaped disks with a (constant) high initial density and Mejecta/Mdisk < 1 partially survive the impact with the nova ejecta. A very minor chemical contamination of the secondary star is anticipated in the U Sco case based on the limited impact of nova ejecta particles on the subgiant in all simulations. Minor mass ejection from the subgiant's outer layers is observed during the late-stage collision with ejecta and disk material, with some particles ejected from the binary system and some accreted by the white dwarf.

GRB~250704B/EP250704a, identified as a short gamma-ray burst (sGRB), exhibited prolonged X-ray emission following the prompt phase and, in optical and infrared (IR) bands, an unusual one-day plateau succeeded by a rapid decline. This sGRB was observed by multiple satellites and ground-based observatories across the electromagnetic spectrum. This study presents temporal and spectral analyses from radio to gamma-ray frequencies, spanning several observation periods beginning after the trigger and continuing for nearly 2 days. The results of the temporal and spectral analyses of the prompt episode, the extended X-ray component, and the afterglow phase are consistent with a millisecond magnetar undergoing accretion. The long-lasting X-ray emission is attributed to the internal energy dissipation of the magnetar spin-down power, governed by the magnetization parameter; the extended optical/IR plateau to synchrotron afterglow emission with energy injection; and the steep decay to changes in microphysical parameters during the post-jet break phase. The X-ray observations are consistent with the superposition of spin-down luminosity and synchrotron afterglow scenario. These findings suggest that the compact-object remnant is most likely a long-lived magnetar.

Nuclear astrophysics aims at unraveling the cosmic origins of chemical elements and the physical processes powering stars. It constitutes a truly multidisciplinary field, that integrates tools, advancements, and accomplishments from theoretical astrophysics, observational astronomy, cosmochemistry, and theoretical and experimental atomic and nuclear physics. For instance, the advent of high-energy astrophysics, facilitated by space-borne observatories, has ushered in a new era, offering a unique, panchromatic view of the universe (i.e., allowing multifrequency observations of stellar events); supercomputers are also playing a pivotal role, furnishing astrophysicists with computational capabilities essential for studying the intricate evolution of stars within a multidimensional framework; cosmochemists, through examination of primitive meteorites, are uncovering tiny fragments of stardust, shedding light on the physical processes operating in stars and on the mechanisms that govern condensation of stellar ejecta into solids; simultaneously, nuclear physicists managed to measure nuclear reactions at (or close to) stellar energies, using both stable and radioactive ion beam facilities. This paper provides a multidisciplinary view on nucleosynthesis accompanying stellar explosions, with a specific focus on thermonuclear supernovae, classical novae, and type I X-ray bursts.

Motivated by recent claims of a possible redshift dependence in late-Universe determinations of the Hubble constant ($H_0$), we test the robustness of this behaviour using multiple cosmological probes. We perform a joint redshift-binned analysis of $H_0$ across eight bins using late-Universe probes -- Pantheon+ SNe~Ia, DESI BAO, cosmic chronometers, and water megamasers -- under three cosmological frameworks: flat $\Lambda$CDM, CPL, and Padé cosmography. Under a common baseline scheme, all three models show a qualitatively similar, low-amplitude variation in the per-bin $H_0$ estimates. A simple Fourier-like parametrization captures this behaviour, but the amplitude differs from zero only at a marginal significance of about $1.71$--$1.94\,\sigma$, with similar behaviour observed across all three cosmological frameworks. We then investigate the robustness and possible origin of this feature. Alternative binning schemes preserve its qualitative form, whereas single-probe per-bin fits (SNe-only, CC-only, BAO-only) yield ratios $H_{0,i}/H_{0,\mathrm{global}}$ mostly consistent with unity and do not reproduce the pronounced drift seen in the joint baseline constraints. Finally, by comparing different global versus piecewise-constant configurations for $\{H_0,\Omega_m,M,r_d\}$, we find that a baseline-like oscillatory pattern re-emerges only when multiple degenerate parameter combinations are allowed to vary across bins, while it is strongly suppressed when only $H_0$ is bin-dependent. Taken together, these results indicate that the apparent oscillatory behaviour of $H_0(z)$ in late-time arises from known parameter degeneracies and does not constitute robust evidence for a genuine redshift evolution.

Deuterated molecules are a useful diagnostic tool to probe the evolution and the kinematics in the earliest stages of star formation. Due to the low temperatures and high densities in the centre of pre-stellar cores, the deuterium fraction is enhanced by several orders of magnitude. We study the distribution of the emission and the deuteration of the two carbon chains HC3N and CH3CCH throughout the pre-stellar core L1544. We analyse emission maps of CH3CCH, CH2DCCH, CH3CCD, HC3N, HCC13CN, and DC3N, observed with the IRAM 30m single-dish radio telescope. We use non-LTE radiative transfer calculations, combined with chemical modelling of the molecular abundances, to constrain physical parameters of the observed species. Following this, we derive the column density and deuteration maps. We find D-fractions of N(DC3N)/N(HC3N)=0.04-0.07, N(CH2DCCH)/N(CH3CCH)=0.09-0.15, and N(CH3CCD)/N(CH3CCH)=0.07-0.09. The deuteration of HC3N appears homogeneous across the core, with widespread D-fraction values above 0.06, tracing intermediate-density gas in the outer layers of the core. CH3CCD is most efficiently formed in the higher-density regions towards the core centre, while the D-fraction of CH2DCCH traces a local density enhancement in the north-east of the core, coinciding with the CH3OH emission peak. The results suggest that gas-phase reactions dominate the formation and deuteration of both HC3N and CH3CCH in L1544, with spatial variations driven by physical structure, density and external radiation. The significantly higher D-fraction of CH2DCCH compared to CH3CCD and a tentative gradient with higher values in the north suggest different deuteration mechanisms for the two functional groups. Similarities between the CH2DCCH emission and CH2DOH might indicate an additional deuteration pathway of CH3CCH on the surfaces of dust grains, as observed for H2CO.

Accurate determination of fireball direction is essential for retrieving trajectories and velocities. Errors in these measurements have significant implications, affecting the calculated pre-impact orbit, influencing mass estimates, and impacting the accuracy of dark flight simulations, where applicable. Here we implement a new atmospheric refraction correction technique that addresses a significant aspect previously overlooked in the field of meteor science. Traditional refraction correction techniques, originally designed for objects positioned at infinite distances, tend to overcompensate when applied to objects within the Earth's atmosphere. To rectify this issue, our study introduces the concept of the atmospheric refraction delta z correction technique, involving the artificial elevation of the observer site height above sea level. We utilize analytically derived formulas for the delta z correction in conjunction with commonly used refraction models, validating these results against a numerical solution that traces light rays through the atmosphere. This ray-tracing model is applied to finely meshed atmospheric layers, yielding precise correction values. We evaluate multiple sources of error in order to quantify the achievable accuracy of the proposed method. Our approach (1) enables the determination of fireball positions with improved astrometric accuracy, (2) removes the explicit dependence on the fireball distance from the observer or its height above Earth's surface within the limits imposed by realistic atmospheric variability, and (3) simplifies meteor data processing by providing a robust framework for analyzing low-elevation fireball observations, for which atmospheric refraction is significant and is automatically corrected by the method. As a result of this work, we provide open, publicly accessible software for calculating the delta z correction.

In this study, we investigate the post-destabilization evolution of a filament in a gravity-balanced model. We adopt the filament model proposed by Solov'ev (2010), in which a dense filament is supported against gravity by the repulsive force between the filament current and its sub-photospheric image. We first performed an analytical investigation of this model. For the numerical study, we use a two-dimensional magnetohydrodynamic (MHD) model that solves the MHD equations with the Lare2d numerical code. Results: In this filament model, analytical expressions are derived for the electric current density, plasma density, and their spatial distributions as functions of the model parameters. The total electric current and the filament weight are also calculated. For the numerical simulations, we constructed an equilibrium filament characterized by a magnetic field of $B_0$ = $10^{-3}$ T, mass density $\rho_0$ ~ 1.3 x $10^{-9}$ kg m$^{-3}$, and temperature T ~ 13000 K. The system was destabilized either by increasing the currents or by reducing the filament density, and its evolution was computed. In both destabilization regimes, the filament was ejected, then halted at a certain altitude, and subsequently fell back, repeating this cycle with a period of about 600 s. The maximum filament ejection velocity was approximately 80 and 40 km $s^{-1}$, respectively. Beneath the ejected filament a current sheet forms, where magnetic reconnection occurs. The maximum ejection altitudes were determined as functions of both the destabilizing currents and the degree of filament plasma dilution. Finally, we compared results of this MHD model with those of an ideal vacuum model and discussed all results.

Reconstructing images of the radio sky from incomplete Fourier information is a key challenge in radio astronomy. In this work, we present a method for radio interferometeic image reconstruction using a data-driven prior for the radio sky based on denoising diffusion probabilistic models (DDPMs). We first train a DDPM on radio galaxy observations from the VLA FIRST survey. We create simulated VLA, EHT, and ALMA observations of radio galaxies, then use an unsupervised posterior sampling method called Denoising Diffusion Restoration Models (DDRM) to reconstruct the corresponding images, using our DDPM as a prior. Our approach is agnostic to the measured radio interferometric data and naturally incorporates the physics of the measurement process. We are able to reconstruct images with very high fidelity PSNR>60, a marked improvement over CLEAN and similar image reconstruction methods using conditional DDPMs

Intrinsic alignments are a major astrophysical contaminant for next generation large-sky surveys like Euclid and LSST. Large hydrodynamic simulations are crucial for informing the alignment modelling for these surveys. We measure position-position and position-shape correlations of a Luminous Red Galaxy sample from the FLAMINGO suite of hydrodynamical simulations, measuring the alignment signal for more than 4.9 million galaxies at redshift 0. We jointly model the clustering and alignment correlations to provide the tightest constraints on the alignment amplitude to date from a hydrodynamic simulation. We find that both the Non-Linear Alignment (NLA) and the more complex Tidal Alignment Tidal Torquing (TATT) models provide good fits to the data. We compare the measured $A_1$ amplitude to observational data and find good agreement. We measure the dependence of the NLA and TATT free parameters on halo mass. We also introduce a mass-dependent TATT model, TATT-M, by finding empirical relations between the halo mass and the TATT parameters. This allows us to fit TATT with only one parameter, $A_1$, with $A_2/A_1$ being a constant and $A_{1\delta}/A_1$ being a function of halo mass. Using a Bayesian approach, we find that TATT-M is very strongly preferred by the data over NLA. Using the baryonic feedback variations of the FLAMINGO simulation suite, we test whether the TATT parameters are sensitive to feedback. Variations in AGN and supernova feedback do not significantly change the alignment amplitude beyond the change associated with the dependence of galaxy stellar mass on the strength of feedback. Our results inform the IA modelling for upcoming surveys by providing guidance on model choices, priors and sensitivities to feedback.

The determination of the solar oxygen abundance remains a central problem in astrophysics, as its accuracy is limited not only by models but also by systematics. While many of these factors have been thoroughly characterized, the effect of the solar activity cycle has so far remained unexplored. Due to its relative strength and accessibility, the O I infrared triplet is typically the primary choice for abundance studies. However, previous investigations have shown that abundances inferred from this triplet tend to be higher than expected on active stars, whereas such an overabundance effect is not observed for the much weaker forbidden O I 6300 Å line. This raises the question of whether a similar trend can be found for the Sun. To address this question, we analyze two decades' worth of synoptic disk-integrated Sun-as-a-star datasets from the FEROS, HARPS-N, PEPSI, and NEID spectrographs, focusing on the infrared triplet (7772, 7774, 7775 Å) and the forbidden O I 6300 Å line. The excellent signal-to-noise ratio of the PEPSI observations allows us to detect a weak but significant variation in the equivalent widths of the infrared triplet, corresponding to about 0.01 dex difference in abundance between activity minimum and maximum. This value is significantly smaller than the typical uncertainties on the solar oxygen abundance. Due to higher scatter, no comparable trend is found in the other data sets. Based on these results, we conclude that within the typical uncertainties presented in other works, we can assume the inferred solar oxygen abundance to be stable across the solar cycle, but that this effect may be significant for other, more active stars.

Machine learning is starting to be used in almost every industry and academic research, and solar physics is no exception. A newly developed machine learning model named MagNet helps us to tackle some of the most serious challenges in data mining by generating transverse fields of solar active regions. Being trained on line-of-sight magnetograms from Michelson Doppler Imager at Solar and Heliospheric Observatory (SOHO/MDI), H{\alpha} maps from Big Bear Solar Observatory and Kanzelhohe Solar Observatory and vector magnetograms from Helioseismic and Magnetic Imager at Solar Dynamic Observatory (SDO/HMI), this model provides vector magnetograms in active regions for SOHO/MDI data covering the strong solar cycle 23. In this study, we performed out-of-sample validation of the MagNet model with data from Imaging Vector Magnetograph (IVM) at Mees Solar Observatory, which was not included in the training process. Our results show good correlation between the AI generated data and the observed vector magnetograms and therefore strengthen the confidence of implementing MagNet to the entire SOHO/MDI archive and future scientific analysis of the AI generated data.

JWST observations provide for the first time evidence for an atmosphere on a rocky exoplanet - 55 Cnc e. The atmosphere of 55 Cnc e is hot with $\text{T}_{\text{eq}}>2000$K and shows strong variability, for which cloud formation above a molten crust could be one possible explanation. The composition of the atmosphere of 55 Cnc e is still unknown but suggests the presence of volatiles. We have run cloud formation models on a grid of N-dominated, O-dominated, C-dominated and H-dominated atmospheres to investigate which type of cloud we could expect on hot super-Earths and hot sub-Neptunes ($1000$K $<$ T $<$ $3000$K). Our models combine radiative transfer with equilibrium chemistry of the gaseous and condensed phases, vertical mixing of condensable species, sedimentation, nucleation and coagulation. We find that the condensability of species is highly dependent on the oxygen abundance of an atmosphere. Oxygen poor atmospheres can be heated by UV and optical absorbers PS, TiO and CN which create temperature inversions. These inhibit condensation. Oxygen rich atmospheres are colder without temperature inversions, and are therefore more favourable environments for cloud formation. The major expected cloud component in O-dominated atmospheres with solar refractory abundance is TiO$_2$(s). Spectral features of clouds in these worlds are stronger in transmission than in emission, in particular at short wavelengths. We find a lack of optical data of solid species in comparison to the variety of stable cloud components which can form on hot, rocky planets.

The thermal Sunyaev-Zel'dovich (tSZ) effect provides a powerful probe of the thermal pressure of ionised gas in galaxy clusters and the cosmic web; constrained simulations reconstruct the mass and velocity fields of the local Universe. We explore how these two may be mutually informative: the tSZ signal provides a benchmark for assessing the fidelity of constrained simulations, and constrained simulations contribute information on the positions, total masses and density profiles of cosmic web structures for use in tSZ studies. We focus on cluster predictions in the Bayesian Origin Reconstruction from Galaxies (BORG) paradigm, introducing CSiBORG-Manticore, a new state-of-the-art suite of digital twins -- data-constrained posterior simulations whose initial conditions are inferred via Bayesian forward modelling. We develop a framework for scoring constrained simulations on their ability to match measured Compton-$y$ maps from Planck for cluster cutouts, and use it to demonstrate improvement from previous BORG reconstructions. We further validate halo masses against weak-lensing-calibrated X-ray masses from eROSITA. We also show how high-fidelity digital twins offer a practical route to extracting additional information from tSZ data through a novel calibration of the mass-observable relation, and provide a complementary framework to purely statistical analyses of Compton-$y$ maps. This paves the way for integrating the large-scale structure information inherent in constrained simulations into the study of CMB secondary anisotropies.

We present deep, high spatial resolution interferometric observations of 0.88 mm continuum emission from the TYC 5709-354-1 system, hereafter WISPIT 2, obtained with the goal of detecting circumplanetary emission in the vicinity of the newly discovered WISPIT 2b planet. Observations with the most extended baseline configuration offered by ALMA, achieving an angular resolution of $25 \times 17$ mas ($3.3 \times 2.2$ au), revealed a single, narrow ring with a deprojected radius of 144.4 au and width of 7.2 au, and no evidence of circumplanetary emission within the cavity. Injection and recovery tests demonstrate that these observations can rule out point-like emission at the location of WISPIT 2b brighter than $\approx 45$ $\mu$Jy at the $3\sigma$ level. While these data can rule out PDS 70c like circumplanetary emission, the upper limit is consistent with empirical mass-flux relationships extrapolated from the stellar regime. Visibility modeling of the continuum ring confirms that WISPIT 2b lies significantly interior to the mm dust ring, raising doubts about the ability of WISPIT 2b to be the only driver of the dust structure. Possible solutions include either another lower mass companion, residing between WISPIT 2b and the cavity edge, likely in the gap seen by SPHERE at $\sim130$ au, or that WISPIT 2b is either substantially more massive than IR-photometry based estimates ($\sim15$ $M_{\rm Jup}$) or on a moderately eccentric orbit. The combination of observations sensitive to the gas and dust distributions on larger spatial scales and dedicated hydrodynamical modeling will help differentiate between scenarios.

We present the Near Infrared Camera (NIRCam) imaging products of the fifth data release (DR5) of the James Webb Space Telescope (JWST) Advanced Deep Extragalactic Survey (JADES). The JADES survey is one of the most ambitious programs yet conducted on JWST, producing deep infrared imaging and multiobject spectroscopy on the GOODS-S and GOODS-N extragalactic deep fields in order to explore galaxies to the earliest epoch. Here we describe the NIRCam data reduction procedures that result in deep and well-characterized mosaics in up to 18 filters covering 469 arcmin$^2$, with 250 arcmin$^2$ having at least 8 filters of coverage. This release contains the full NIRCam imaging of JADES, over 800 JWST mission hours, as well as co-reductions of 19 other programs in these two premier deep fields. We perform detailed tests on the final data products, thereby characterizing the photometric properties, point-spread function, and astrometric alignment. We release mosaics for individual programs (or epochs, depending on scheduling) and the mosaics combining data from all programs in order to facilitate photometric variability studies and the deepest possible photometry.

Medium to ultra-deep mid-infrared imaging surveys with the James Webb Space Telescope (JWST)'s Mid-Infrared Instrument (MIRI) are reframing our view of the early Universe, from the emergence of ultra-red dusty and quiescent galaxies to the epoch of reionization to the first galaxies. Here we present the MIRI coordinated parallels component of the JADES program, which obtained ultra-deep (155 ks) imaging at $7.7 \mu$m over $\sim10$ arcmin$^2$ as well as medium depth ($\sim5-15$ ks) imaging at $7.7, 12.8$, and $15 \mu$m over $\sim36$, 25, and 22 arcmin$^2$, respectively, in the GOODS-S and GOODS-N fields. This paper describes the data reduction, which combines the official JWST Calibration Pipeline with custom steps to optimize flagging of warm/hot pixels and optimize background subtraction. We further introduce a new step to address artifacts caused by persistence from saturating sources. The final, fully reduced JADES/MIRI mosaics are being released as part of JADES Data Release 5, along with prior-based forced photometry using NIRCam detection images, providing critical rest-frame near-infrared and optical constraints on early galaxy populations.

JADES Data Release 5 (DR5) photometric catalogs and describes the methodologies used for source detection, deblending, photometry, uncertainty estimation, and catalog curation. The catalogs are constructed from 35 space-based imaging mosaics obtained with JWST/NIRCam, JWST/MIRI, HST/ACS, and HST/WFC3, combining approximately 1250 hours of JADES imaging with extensive additional public JWST and HST observations in the GOODS fields. Sources are identified using custom signal-to-noise-based detection and deblending algorithms optimized for the depth, resolution, and complex point-spread-function structure of JWST imaging. Source centroids, shapes, and photometric apertures are determined using a new fast two-dimensional Gaussian regression method applied to detection-image profiles. We provide forced circular-aperture photometry, ellipsoidal Kron photometry, and curve-of-growth measurements for every source in every band. We introduce a new pixel-level regression framework to model photometric uncertainties as a function of aperture size and local mosaic properties, accounting for correlated noise in heterogeneous JWST mosaics. Photometric redshifts are computed using template-based fitting applied to both small-aperture photometry on unconvolved images and Kron photometry on common-PSF mosaics. The JADES DR5 catalogs supersede previous JADES photometric releases, and are publicly released through the Mikulski Archive for Space Telescopes and an interactive web interface.

We present morphological parameters and their uncertainties for all sources detected in JWST/NIRCam imaging in GOODS-N and GOODS-S from the JWST Advanced Deep Extragalactic Survey (JADES) catalogs. We model the surface brightness profiles of these sources with single-component Sérsic profiles, performing Bayesian inference of galaxy structural parameters. We fit each of the $>10^5$ sources with every available JWST/NIRCam wide-band filter individually, amounting to over 3 million Sérsic profiles computed. We provide catalogs of this morphological information, building one of the largest extragalactic morphological datasets to date, which we share alongside imaging and photometry from the JADES Data Release 5. With this information, we analyze the rest-frame optical redshift evolution of the effective radius and the surface luminosity density within a radius of 1 kiloparsec, $\Sigma_{\text{1 kpc}}$, for 24,692 galaxies at $z>1$. We find $r_{\text{eff}} \propto (1+z)^{-0.635 \pm 0.013}$ kpc, while $\Sigma_{\text{1 kpc}}$ is relatively constant across time. Additionally, we explore bulge-disk decomposition on a subset of 8,390 galaxies in the JADES deep imaging covering the Hubble Ultra Deep Field, finding the effective radius of the bulge-components to increase marginally with time, whereas the disk-component sizes evolve as $r_{\text{eff,disk}} \propto (1+z)^{-1.091 \pm 0.043}$. Future work modeling multi-component surface brightness profiles will enable further analysis of the morphological evolution of galaxies across cosmic time.

Wisps are among the most prominent scattered light artifacts in JWST/NIRCam imaging. They often appear in certain regions of the detectors and contaminate observations at surface-brightness levels relevant for faint-source photometry. We introduce a new subtraction method that uses the non-negative matrix factorization (NMF) algorithm to model and remove wisps. Using deep NIRCam observations from the JWST Advanced Deep Extragalactic Survey (JADES) and other programs, we construct multi-component, filter- and detector-specific wisp templates that capture the wisp structures and their exposure-to-exposure morphological variations. Wisps in individual exposures are represented as non-negative linear combinations of these templates, consistent with their additive nature and reducing degeneracies relative to single-template scaling. Compared to existing approaches, our method delivers lower residual root mean square in wisp-affected regions and reduces photometric bias and scatter to levels consistent with clean detector areas. The NMF wisp templates are readily applicable to other datasets and are publicly released to support future NIRCam extragalactic surveys.

We present a sample of 2081 sources selected at photometric redshift $z_{\mathrm{phot}} > 8$ across the JADES DR5 data release in GOODS-S and GOODS-N over a total area of 469 square arcmin. These sources range from $M_{\mathrm{UV}} = -22$ to $M_{\mathrm{UV}} = -16$, with 19 objects at $z_{\mathrm{phot}} > 14$. We estimate the UV slopes for the full sample from fits to the photometry and find evidence for a steepening of the relationship between the UV continuum slope and $M_{\mathrm{UV}}$ to higher redshifts, a result that differs from prior analyses of brighter samples in the literature. We provide evidence that over one quarter of our sources have evidence for being morphologically extended, with many galaxies showing multiple bright knots or clumps even out to $z \sim 13 - 14$, an indication of how galaxies at Cosmic Dawn are growing and evolving. We discuss JADES-GN+189.15982+62.28899, a GOODS-N F200W dropout galaxy at $z_{\mathrm{phot}} \sim 15 - 18$ which has been observed spectroscopically with JWST/NIRSpec in prism mode, resulting in a very low signal-to-noise spectrum that is consistent with the photometry and rules out a number of low-redshift solutions for the source. Finally, we use a subsample of 123 objects in our sample with spectroscopic redshifts to explore the usage of alternate fitting templates and a prescription for Ly-$\alpha$ damping wing absorption, finding that both produce significant improvements to the estimated photometric redshifts.

We report a galaxy overdensity candidate at $z\approx 10.5$ in the JWST Advanced Deep Extragalactic Survey (JADES). This overdensity contains 18 galaxies with consistent photometric redshifts and robust F115W dropouts within 8 comoving Mpc in projection. The galaxy number density is four times higher than the field expectation, accounting for one-third of comparably bright galaxies and nearly 50% of the total star formation rate at $10<z_\mathrm{phot}<12$ in the GOODS-S field. Two compact members of the overdensity show potential Balmer breaks suggestive of evolved stellar populations or little red dots (LRDs). One-third of galaxies have close companions or substructures within 1 kpc at consistent photometric redshifts, implying more frequent interactions in an overdense environment. Most galaxies have stellar masses of 0.6-3$\times10^8$ $M_\odot$, half-light radii of $\sim$200 pc, and star formation rates of $\sim$5 $M_\odot \mathrm{yr^{-1}}$, with no significant deviation from typical high-redshift scaling relations. We find tentative evidence for a spatially varying Ly$\alpha$ transmission inferred photometrically, consistent with an emerging ionized bubble. This overdensity provides a rare opportunity for probing the environmental impact on galaxy evolution and the onset of cosmic reionization within the first 500 Myr.

Using JWST NIRCam imaging and grism spectroscopy from the JWST Advanced Deep Extragalactic Survey (JADES) Origins Fields, we report spectroscopic redshift measurements of 1,445 emission-line galaxies at $z=0-9$. Within this sample, we identify two prominent galaxy protoclusters at $z = 3.47$ and 3.69, each anchored by massive dusty star-forming galaxies (DSFGs). In the vicinity of these systems, we discover seven background galaxies at $z=3.6 - 6$ that simultaneously exhibit strong rest-frame optical emission lines (e.g., [O III] and H$\alpha$) and unusually reddened UV-to-optical continua. We attribute this reddening to dust extinction arising from the circumgalactic medium (CGM) of the foreground DSFGs at projected separations of 7-30 kpc. We infer a high dust column density ($\gtrsim 10^{-1}$ Msun/kpc^2), substantially exceeding those measured in low-redshift halos and those predicted by hydrodynamical simulations like IllustrisTNG and FIRE-2. The steep extinction curves, comparable to or steeper than that of the SMC, indicate a dominant population of small dust grains in the high-redshift CGM. We conclude that DSFGs at this epoch host large reservoirs of dusty CGM enriched to solar metallicity. These extended dust components are largely invisible to (sub-)millimeter interferometers such as ALMA because of their low surface brightness. We discuss the physical processes in dust transport that might be key to reproducing our observations, including galaxy mergers, cool-phase gas outflows, dust shattering, sputtering and radiation pressure. Finally, we caution that foreground CGM dust extinction may redden background galaxies at intermediate redshifts to mimic Lyman-break galaxies at $z\gtrsim10$.

In the local Universe, supermassive black hole (SMBH) masses strongly correlate with their host-galaxies' stellar masses ($M_{*}$), but galaxies hosting faint AGN recently found by JWST may deviate from this relation. To constrain the M$_{\text{BH}}$-M$_{*}$ relation at high redshift, we performed AGN-host image decomposition for 17 low-luminosity AGN galaxies at $z$\,$\sim$\,4--6 using NIRCam images in the JADES GOODS-N field. These sources are identified as AGNs from broad H$\alpha$ emission lines detected by the CONGRESS and FRESCO surveys. We used \textsc{galfit+MCMC} to fit spatial profiles in 7 wide-band images and detected extended emission in 9 sources out of 17. The close spatial alignment between the extended components and the AGN centers indicates that this emission likely originates from the host galaxies. These sources are extended at 0.9--2.0~$\mu$m, suggesting significant host-galaxy light in the rest-frame UV. For the sources with the host detection, the stellar mass inferred based on image decomposition result can be 1-2 dex lower than the results without image decomposition. The BH-to-stellar mass ratio spans $M_{\text{BH}}/M_\ast$\,$\sim$\,0.01--1.48, placing them well above the local $M_{\text{BH}}$--$M_\ast$ relation. In contrast, the host-galaxy size--mass relation broadly agrees with previous measurements. Our results suggest that the host galaxies of these faint AGN are either genuinely under-massive compared to their black hole masses, or too compact to be spatially resolved.

Dust attenuation modifies the observed spectral energy distribution (SED), leading to biases in the properties inferred from integrated SED fitting. As spatially resolved SED modeling becomes feasible for large high-redshift samples, it is increasingly important to assess how dust attenuation affects resolved mass estimates. We evaluate the impact of dust attenuation on stellar mass estimates derived from integrating spatially resolved SED fitting results. We perform spatially resolved and integrated SED fitting on a sample of 3408 galaxies at $3 \leq z < 9$ from the GOODS South field, combining deep NIRCam from the JWST Advanced Deep Extragalactic Survey (JADES) and HST/ACS imaging from GOODS and CANDELS. We compare galaxy-integrated properties derived from fitting the summed SED with those obtained from spatially resolved SED modeling. Using a two-component dust attenuation model with a variable slope, we investigate how the dust attenuation slope, A(V), and stellar population properties contribute to discrepancies in the resulting stellar mass estimates. Resolved stellar masses are systematically higher than integrated estimates, with a median offset of +0.24 dex. Resolved analyses recover higher dust attenuations ($\Delta A(V)\approx +0.08$ mag), lower birth cloud fractions ($\Delta\mu \approx -0.28$), and grayer attenuation curves ($\Delta\delta_{\mathrm{ISM}} = +0.08$), arising from preferential sampling of compact star-forming regions. Integrated fits underestimate stellar ages by $\sim23\%$ at $z < 5$ and 31$\%$ at $z \gtrsim 5$. The stellar mass offset correlates strongly with the age difference and the attenuation slope difference, indicating that age-dependent outshining and spatially varying dust geometry are primary drivers of the discrepancy between resolved and integrated stellar masses.

We present nitrogen abundance measurements based on the low-ionisation [NII]6583 emission line for 588 galaxies between 1.5<z<7.0 from the JWST Advanced Deep Extragalactic Survey (JADES). We detect the temperature-sensitive [OIII]4363 auroral line in 40 galaxies in our sample, affording $T_e$-based abundances for this subset. We find that the average N/O abundance ratio in our low-metallicity sample is at least 0.1 dex higher than z ~ 0 samples. In particular, we find significant scatter toward high N/O, with five galaxies being identified with enhanced nitrogen abundances (log(N/O)>-1.1) at low-metallicity (12+log(O/H)<8.0) from $T_e$-based measurements. Meanwhile, applying strong-line abundance measurements to the remainder of our sample reveals a further 14 candidate galaxies passing these abundance cuts, implying that around 13 % of 12+log(O/H)<8.0 galaxies at these redshifts are nitrogen-enhanced at this level. We find that N/O abundance in low-metallicity systems correlates with SFR, surface density of SFR, and surface density of stellar mass at high redshift, while only in high-metallicity systems does a correlation with stellar mass emerge. Despite healthy representation of these `moderately nitrogen-enhanced' galaxies (-1.1<log(N/O)<-0.6), no galaxies in our low-metallicity sample are identified as having log(N/O)>-0.6, abundances that are typical of high-redshift NIII]- and NIV]-emitters. This demonstrates that the extreme nitrogen enhancements seen in some NIII]- and NIV]-emitters are only attained during the most extreme starbursts. This suggests that these elevated abundances are caused by enrichment from young massive stars in extreme environments and that the impact of this enrichment pathway is milder, though still important, for high-redshift systems on the star-forming main sequence.

Massive star-forming clumps are a prominent feature of high-redshift galaxies and are thought to trace gravitational fragmentation, feedback, and bulge growth in gas-rich disks. We present a statistical analysis of clump-like structures in $\sim$3600 galaxies spanning $2 \lesssim z \lesssim 8$ from deep JWST/NIRCam imaging in the JADES GOODS--South field. Clumps are identified as residual features after subtracting smooth Sérsic profiles, enabling a uniform, rest-frame optical census of sub-galactic structure. We characterize their physical properties, size--mass relations, and spatial distributions to constrain models of sub-galactic structure formation and evolution. We find that clumps in our sample are typically low-mass ($10^{\sim7-8}M_\odot$), actively star-forming, and show diverse gas-phase metallicity, dust attenuation, and stellar population properties. Their sizes and average pairwise separations increase with cosmic time (toward lower redshift), consistent with inside-out disk growth. The clump mass function follows a power law with slope $\alpha = -1.50_{-0.17}^{+0.19}$, consistent with fragmentation in turbulent disks. We find a deficit of relatively young clumps near galaxy centers and a radial transition in the size--mass relation: outer clumps exhibit steeper, near-virial slopes ($R_{\rm e}\propto M_*^{\sim 0.3}$), while inner clumps follow flatter trends ($R_{\rm e}\propto M_*^{\sim 0.2}$), consistent with structural evolution via migration or disruption. These results provide new constraints on the formation, survival, and dynamical evolution of clumps, highlighting their role in shaping galaxy morphology during the peak of cosmic star formation.

The microlensing event Gaia20fnr is a long-duration, non-caustic-crossing binary-lens event at high Galactic latitude. Triggered by a photometric rise detected by the Gaia space mission, the event was followed up with observations from multiple ground-based facilities and four space telescopes: Gaia, NEOWISE, Swift, and TESS. We characterize the Gaia20fnr microlensing system by determining the physical and orbital properties of the binary lens, the nature of the luminous source, and the kinematics of both the source and the lens. We employed a binary-lens microlensing model including full Keplerian orbital motion and annual microlens parallax to fit the photometric data. The event is best explained by a K2 giant source at $D_{\rm S} = 3.10 \pm 0.10\,\mathrm{kpc}$ lensed by a stellar binary composed of $M_{\rm L,1} = 0.46 \pm 0.06\,M_\odot$ and $M_{\rm L,2} = 0.52 \pm 0.06\,M_\odot$ at a distance of $D_{\rm L} = 0.54 \pm 0.05\,\mathrm{kpc}$. The light curve exhibits strong signatures of orbital motion and requires a full Keplerian model with a period of $P = 0.67 \pm 0.04\,\mathrm{yr}$ and a radial-velocity semi-amplitude of $K_1 = 16.9 \pm 0.9\,\mathrm{km\,s^{-1}}$. Gaia20fnr is one of the few microlensing events for which a complete Keplerian binary-lens solution has been derived. The model can be tested with follow-up radial-velocity and high-resolution imaging observations as well as forthcoming Gaia DR4 and DR5 astrometric time-series data. Its long duration, multi-peak structure, and extensive coverage make it a benchmark for studying faint nearby low-mass binaries through microlensing.

We present the discovery of HE0144-4657, an ultra metal-poor, CNO-enhanced star dynamically associated with the Helmi Stream disrupted dwarf-galaxy remnant. This star was first identified as a carbon-enhanced, metal-poor star candidate from the Hamburg/ESO objective-prism survey, then followed up with medium- and high-resolution spectroscopy. At [Fe/H]=-4.11, HE0144-4657 is the lowest metallicity star found in a stellar stream to date. Its chemistry is consistent with field halo stars in the same metallicity regime, and the light-element (atomic number Z<=30) chemical abundance pattern suggests that HE0144-4657 is a bona-fide second-generation star with a possible Population III progenitor in the 50Msun mass range with low explosion energy. One possible scenario for the origin of HE0144-4657 is that it was formed in an ultra-faint dwarf galaxy accreted by the Helmi Stream progenitor system before merging with the Milky Way. This discovery provides further evidence for the extragalactic origin of carbon-enhanced ultra metal-poor stars in the Milky Way and for the specific environments conducive to their formation.

IRAS23226-3843 has been identified as a highly variable Seyfert galaxy and even as a changing-look active galactic nucleus based on optical spectra. Here we present follow-up observations - taken over the past five years - for examining the ongoing photometric and spectral variations in this remarkable galaxy. We carried out SWIFT observations of IRAS23226-3843 together with new optical spectra taken in 2023 and 2024. In parallel we investigate ASAS-SN photometric data from 2014 till 2025. IRAS23226-3843 stayed on a high continuum flux level in the X-ray as well as in the optical since a historic outburst in 2019. However, it shows strong short-term variations on timescales of a few months. Densely sampled ASAS-SN V-band continuum data from 2014 till 2025 confirm that behavior. IRAS23226-3843 switched from a clear Seyfert 1 type in December 2019 to a Seyfert 1.9/2 type in July 2020 based on its optical spectra. Afterward, it again became a Seyfert 1 type with symmetric broad single-peaked Balmer line profiles in January 2023. These spectra prove the repeating changing-look character of the galaxy.IRAS23226-3843 exhibits extreme high Balmer decrements Ha\Hb based on their broad line components. The Balmer decrement values are on the order of 10. IRAS23226-3843 successively showed all types of broad line Balmer profiles during the past 25 years over periods of many years: asymmetric single-peaked, double-peaked, as well as single-peaked and symmetric profiles in addition to its Seyfert 1.9/2 transition. These variations are not clearly correlated with continuum and line intensity variations.

The presence of multi-phase outflows in low ionisation nuclear emission-line regions (LINERs) has been confirmed to be frequent, but the mechanisms that launch them are still under study. We aim to explore the connections between the ionised gas outflow, radio continuum structures and X-ray emission detected in the LINER NGC4438. We analyse L, C and X-band images (from 1.4 to 12 GHz) of the LINER NGC4438, combining high-resolution data from enhanced Multi Element Radio Linked Interferometer Network (e-MERLIN) and Karl G Jansky Very Large Array (VLA). We produce radio flux, spectral index maps, and an energetic model that allows us to characterise the source. We incorporate optical integral field spectroscopy (IFS) data (GTC/MEGARA) and Chandra X-ray data, with comparable resolution, to better trace the outflow, the AGN and their potential connection. We present new L, C, and X-band high-resolution, high-sensitivity radio images and spectral-index maps that probe $\sim$ 25 pc scales in NGC 4438. These data reveal a close morphological correspondence between the radio structures and the ionised gas bubble. Using a spatially resolved energetic model based on radio flux and spectral index, we disentangle the compact AGN emission from the extended bubble for the first time, establishing their distinct physical origins. We measure a kinetic power of $\sim 5\times 10^{44}$ erg s$^{-1}$ for the radio bubble, exceeding the power of the ionised outflow by more than three orders of magnitude. Our multi-wavelength analysis indicates that NGC 4438 is undergoing jet-mode feedback, where a low-luminosity, weakly collimated jet impacts the dense northern interstellar medium. This interaction drives shock-ionised gas, produces a moderate velocity outflow that removes material from the region, and generates thermal X-ray emission coincident with the radio and H$\alpha$ cavity.

On December 6, 2024, 1ES 1959+650, one of the X-ray brightest blazars known, underwent a high-amplitude X-ray outburst detected by SVOM, the first such discovery with this mission. The source was subsequently monitored with SVOM and Swift from December 2024 to March 2025. We report the detection and multi-wavelength follow-up of this event, and describe the temporal and spectral evolution observed during the campaign. Data from SVOM/MXT, SVOM/ECLAIRs, and Swift/XRT were analyzed with log-parabola models to track flux and spectral variability. The source was detected in a bright state over the 0.3-50 keV range. During the three months of monitoring, the X-ray flux varied significantly, showing episodes of spectral hardening at high flux levels. The spectral curvature evolved more irregularly and did not show a clear trend with flux. A shift of the Spectral Energy Distribution (SED) synchrotron peak to higher energies is seen when the flux increases. This constitutes the first blazar outburst discovered in X-rays by SVOM. The coordinated follow-up with Swift provided continuous coverage of the flare and highlights the strong complementarity of the two missions for time-domain studies of blazars. The flare shows no clear signatures of either Fermi I or Fermi II acceleration, suggesting a mixed Fermi I/II scenario.

Although galaxy evolution is governed by the interplay between baryonic physics and dark matter halo assembly, how halo properties shape observed galaxies remains unclear. With current challenges in measuring halo properties, the stellar-to-total dynamical mass relation is introduced as an alternative metric sensitive to the dark matter content within galaxies. We explore how spatially resolved stellar population properties vary across this relation using optical IFS data and photometry from 265 CALIFA galaxies. Spatially resolved ages and metallicities, [M/H], are derived using a Bayesian framework fed with a library of model spectra based on stochastic star formation and metallicity histories and dust attenuation. We study these properties in terms of both stellar and total dynamical mass, with the latter being enclosed mass within three effective radii from Jeans dynamical modeling. We find that ages and [M/H] measured at different annuli depend on both stellar and total mass, yet showing distinct radial trends. While the dependence of age on total mass is more prominent in the outskirts, that of [M/H] is significant in the inner parts. This behavior is reflected in the stellar population profiles and gradients, more strongly for age and connected to morphology. Intermediate-mass early-types have higher stellar-to-total mass ratios and flatter age profiles with older ages, and steep negative [M/H] profiles, whereas later-types have lower stellar-to-total mass ratios, negative age profiles with younger ages and shallower negative [M/H] profiles. Moreover, at fixed stellar mass galaxies have more negative age gradients and shallower [M/H] ones as total mass increases. Our results show that total dynamical mass is linked to systematic variations in stellar populations and radial gradients at fixed stellar mass, suggesting a relevant role of dark matter halos in shaping galaxy properties

Unresolved stellar companions can cause both under-estimations in the radii of transiting planets and over-estimations of their detectability, affecting our ability to reliably measure planet occurrence rates. To quantify the latter, we identified a control sample of 198 Kepler stars with sensitivity to Earth-like planets if they were single stars, and imaged them with adaptive optics. In 20% of systems, we detected stellar companions that were close enough to go unresolved in Kepler observations. We calculated the distribution of planet radius correction factors needed to adjust for these observed companions, along with simulations of undetected companions to which our observations were not sensitive. We then used these correction factors to optimize an occurrence rate model for small close-in planets while correcting Kepler's detection efficiency for the presence of unresolved companions, and quantified how this correction affects occurrence estimates. Median occurrence rates for small planets between $2-100$ days increased by an average factor of $1.08-1.19$ (depending on statistical treatments), with the largest differences found for smaller planets at larger orbital periods. We found that the frequency of Earth-sized planets in the habitable zone ($\eta_\oplus$) increased by a factor of ${1.18}_{-0.66}^{+0.43}-{1.46}_{-0.83}^{+0.53}$ when accounting for the effect of unresolved companions on Kepler's detection sensitivity.

We present measurements of the cool molecular gas mass around the nuclei of two gas-rich mergers, III Zw 035 and IRAS F01364-1042, whose enclosed masses (M$_\mathrm{enc}$) within the central 40-80 pc would be overmassive if attributed entirely to the supermassive black hole mass (SMBH) and compared to SMBH-galaxy scaling relations. Our gas mass measurements are derived from Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 long-baseline observations of CO(J=2-1) and 230 GHz continuum emission at 14-20 pc resolution, which probes below the resolving limit of the previous black hole mass measurements. Subtracting molecular gas mass from these enclosed masses is not enough to reconcile with BH-galaxy relationships, but independently measuring M$_\mathrm{enc}$ using the cold CO(2-1) gas does shift the black holes down to their expected values. Still, these ALMA data reveal respective molecular gas masses of $\sim$3$\times$10$^7$ to $\sim$6$\times$10$^8$ M$_\odot$ within 70 pc of these black holes, which could challenge some black hole accretion models that assume nuclear gas like this has no angular momentum.

We use radiation-hydrodynamical simulations to investigate the formation and synthetic X-ray emission of hot bubbles within planetary nebulae (PNe) driven by the powerful winds of H-deficient, [Wolf-Rayet]([WR])-type stars. Our models, based on {\sc mesa} stellar evolution tracks for 1--3 M$_{\odot}$ progenitors, adopt a recent mass-loss rate prescription for [WR] stars and incorporate the enhanced radiative cooling of their C-rich material, comparing the results against standard H-rich PN models. The enhanced mass-loss in the [WR] models leads to an accelerated post-AGB evolution and a subsequent delay in hot bubble formation compared to their H-rich counterparts, as suggested by a previous work. By computing synthetic X-ray spectra that account for the mixed H-rich and H-deficient gas phases, we find that models incorporating [WR] winds exhibit significantly higher X-ray luminosities ($L_\mathrm{X}$) than their H-rich counterparts, but the emissivity-weighted plasma temperature of the X-ray-emitting gas converge to values of $T_\mathrm{X} = [1-3] \times 10^{6}$~K, regardless of whether the system follows a [WR]-type or an H-rich post-AGB evolutionary path. Our results reinforce previous suggestions that mixing is a key mechanism in generating the observed soft X-ray emission even for PN hosting [WR] central stars.

Tunnelling reactions of molecules embedded on cryogenic noble-gas matrices are being used in fundamental studies of how reactivity varies with the nature of the supposedly inert matrix as well as pointers to the chemistry occurring in the interstellar medium on ice-grains. To these ends we present chemical kinetic rate constants for the \textit{cis} to \textit{trans} isomerisation of seleno-, thio- and monomeric formic acids and that of their three dimeric species, based on multidimensional calculations in the gas-phase, from 10~K to 300~K as a guide to the matrix reactions.

Magnetars are extreme neutron stars powered by ultra-strong magnetic fields ($\sim10^{14}$ Gauss) and are compelling engines for some of the most powerful extragalactic transients such as Super Luminous Supernovae, Gamma-Ray Bursts, and Fast Radio Bursts. Yet their formation rate relative to ordinary neutron stars remains uncertain, often precluding direct comparisons with the rates of these extragalactic transients. Furthermore, magnetars have been recently shown to be evolutionarily related to other neutron star classes, complicating the estimate of the exact magnetar fraction within the neutron star population. We study the magnetar birth fraction in core-collapse supernovae using pulsar population synthesis of all isolated neutron star classes in our Galaxy, incorporating self-consistently the Galactic dynamical evolution, spin-down and magneto-thermal evolution. This approach allows us to derive strong constraints from small close-to-complete observational samples. In particular, looking at the age-limited young ($<$2 kyr) neutron star population in the Milky Way we find 24 detected young neutron stars, with only 10 of them (41%) being classical rotational powered pulsars, while the others (59%) are either magnetars or central compact objects, the latter believed to be equally magnetically powered. We further compare the results with the nearby volume-limited class ($<$500 pc) of X-ray Dim Isolated Neutron stars, old nearby magnetars. We conclude that the observed population of isolated neutron stars in the Galaxy can be reproduced only by assuming a core-collapse supernova rate larger than two, and a larger magnetar fraction than previously inferred. By assuming a bimodal initial magnetic field ($B_0$) distribution at birth, we find that the magnetar class peaks between $B_0\sim 1-2.5\times10^{14}$ Gauss and represents on average $\sim50$% of the entire neutron star population.

We present cosmological parameters measurements from the full combination of DESI DR1 galaxy clustering data described with large-scale structure effective field theory. By incorporating additional datasets (photometric galaxies and CMB lensing cross-correlations) and extending the bispectrum likelihood to smaller scales using a consistent one-loop theory computation, we achieve substantial gains in constraining power relative to previous analyses. Combining with the latest DESI baryon acoustic oscillation data and using cosmic microwave background (CMB) priors on the power spectrum tilt and baryon density, we obtain tight constraints on the $\Lambda$CDM model, finding the Hubble constant $H_0=69.08\pm 0.37~\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, the matter density fraction $\Omega_m=0.2973\pm 0.0050$, and the mass fluctuation amplitude $\sigma_8 = 0.815\pm 0.016$ (or the lensing parameter $S_8\equiv\sigma_8\sqrt{\Omega_m/0.3}=0.811\pm 0.016$), corresponding to $0.6\%$, $1.7\%$, and $2\%$ precision respectively. Adding the Pantheon+ supernova sample (SNe), we find a preference of $2.6\sigma$ for the $w_0w_a$ dynamical dark energy model from low-redshift data alone, which increases to $2.8\sigma$ when exchanging the SNe with Planck CMB data. Combining full-shape data with BAO, CMB, and SNe likelihoods, we improve the dark energy figure-of-merit by $18\%$ and bound the sum of the neutrino masses to $M_\nu<0.057$ eV in $\Lambda$CDM and $M_\nu<0.095$ eV in the $w_0w_a$ dynamical dark energy model (both at 95\% CL). This represents an improvement of $25\%$ over the background expansion constraints and the strongest bound on neutrino masses in $w_0w_a$CDM to date. Our results suggest that the preference for the normal ordering of neutrino mass states holds regardless of the cosmological background model, and is robust in light of tensions between cosmological datasets.

FarView is a proposed low frequency radio interferometer for deployment on the lunar far side, enabled by the Moon's radio quiet environment. Operating over 1-50 MHz inaccessible from Earth, FarView will open a new observational window and promote discovery class science in cosmology, heliophysics, Galactic and exoplanet astrophysics. The primary science is measurement of the redshifted 21 cm signal from the Cosmic Dark Ages (z=30-100), identified by the Astro2020 Decadal Survey as a priority cosmology discovery area. FarView will deliver 3D tomographic measurements and precision power spectra of neutral hydrogen in a largely linear regime, enabling tests of inflationary initial conditions, primordial non Gaussianity, dark matter properties, neutrino masses, and early dark energy. The reference design consists of 100000 crossed dipole antennas in a dense core-halo configuration spanning 200 sq km. A compact 4 km core with 83000 dipoles maximizes sensitivity to large scale cosmological modes, while 20000 halo elements extending to 14 km provide angular resolution and calibration for foreground characterization. Sensitivity forecasts indicate a 10-sigma detection of the Dark Ages 21 cm power spectrum at z=30 over five years of half duty cycle lunar night observations. An FFT-based EPIC beamformer is identified as an efficient signal processing architecture. Beyond cosmology, FarView will enable interferometric imaging of low frequency solar radio bursts, advancing space weather studies. Additional capabilities include stellar space weather observations, Galactic cosmic ray tomography via free-free absorption, and searches for auroral radio emission from exoplanet magnetospheres, a probe of exoplanet habitability. FarView represents a flagship class opportunity to establish the Moon as a platform for foundational astrophysics while delivering unique observational capabilities.

Dynamical modelling is a fundamental tool for measuring galaxy masses and density profiles in the era of large integral-field spectroscopic surveys and Bayesian inference. Solutions based on the Jeans equations are popular due to their robustness and computational efficiency. However, traditional semi-analytic Jeans solvers often require restrictive assumptions about the velocity anisotropy to remain computationally tractable. This paper presents a new spectral solver for the axisymmetric Jeans equations designed to overcome these limitations. I first illustrate, using orbit integrations in realistic potentials, that spherical alignment of the velocity ellipsoid is a physically well-motivated approximation for galaxy modelling. The new method employs a spectral technique to solve the Jeans partial differential equations directly. Two design choices are critical for accuracy and speed: (i) solving for the slowly-varying velocity dispersion rather than the rapidly varying pressure, and (ii) imposing a Robin boundary condition to enforce the asymptotic decay on a finite domain. This formulation supports arbitrary anisotropy distributions beta(r, theta) while simultaneously increasing computational speed by orders of magnitude compared to standard high-accuracy quadratures. Validated against exact analytic benchmarks, the solver recovers intrinsic moments with sub-percent accuracy. The implementation will be included in the public JamPy package and is structured to be optimally suited for massive parallelization on specialized hardware such as GPUs, enabling the rigorous exploration of complex parameter spaces.

As the recurrent nova T Coronae Borealis (T CrB) approaches its next predicted thermonuclear eruption, it is currently exhibiting a "super-active state" (SAS) characterized by enhanced multiwavelength emission similar to the behavior recorded prior to the 1946 outburst. We present a multiwavelength analysis of the SAS and the subsequent "faint state" using observations from HST, Swift, NuSTAR, and XMM-Newton. Our results indicate that the SAS was driven by an increase in the mass accretion rate, which caused the accretion disk's boundary layer to become optically thick. A weighted least squares regression analysis quantifies the evolution of the accretion components, displaying a highly significant (4.5$\sigma$) increase in the luminosity of the optically thin cooling flow (L$_{cf}$) and a marginal (2.58$\sigma$) decrease in the optically thick boundary layer luminosity (L$_{bb}$) as the system transitioned into the faint state. We find that this dimming is consistent with an intrinsic change in the accretion flow rather than dust obscuration, supported by the lack of infrared excess and the stability of the 2175 Å feature. Additionally, a time-series analysis using autoregressive modeling to account for correlated red noise revealed no significant periodicities, thereby disputing the previously reported $\sim$6000 s signal. These findings suggest that the pre-outburst evolution of T CrB is characterized by significant changes in the accretion disk structure and boundary layer, providing a self-consistent physical framework for the system's behavior as it approaches eruption.

Analyses of type Ia supernovae have helped us shed light on the existence and nature of dark energy. Most of these analyses have relied on Bayesian techniques. In this work, we rely on resampling techniques to analyse supernova data. In particular, we use the generalised least squares method together with Jackknife and Bootstrap techniques to estimate parameters of $\Lambda$CDM, flat $\Lambda$CDM, $w$CDM, flat $w$CDM, and flat $w_0\,w_a$CDM models from the recent PantheonPlus and SH0ES data. For completeness, we also perform Bayesian analysis using Markov chain Monte Carlo (MCMC) and nested sampling algorithms, and compare the results. We note that resampling techniques can help highlight the limitations of the data. For instance, we see that the Jackknife method estimates a strong positive correlation between $h$ and $M$ and higher standard deviations for both. This may have significant implications for the Hubble tension. We conclude with a discussion of our results.

We present Atacama Large Millimeter/Submillimeter Array observations of multi-wavelength dust emissions at 3.1\,mm and 1.3\,mm; along with molecular line emissions of CO(2--1), CO(3--2), \mbox{$^{13}$CO(3--2)}, and C$^{18}$O(3--2) at spatial resolutions of 7--45 AU towards the protoplanetary system AS 205. The dust emissions exhibit two distinct components of AS 205 N and AS 205 S, separated by 1.3 arcsec. While gas kinematics within the dust disk regions are dominated by Keplerian rotation, the more extended gas emission displays complex morphology and kinematics strongly affected by the binary gravitational interaction in the outer regions. The stellar masses of AS 205 N and AS 205 S are estimated at $0.78\pm0.19$\,M$_\odot$ and $1.93\pm0.86$\,M$_\odot$, respectively. Azimuthal variation is observed in the spectral index distribution of both disks. In AS 205 N, the spectral index minimum in the southwest is coincident with the peaks of CO($2-1$), CO($3-2$), and $^{13}$CO($3-2$) integrated intensity and the relative position of its southern counterpart. On the other hand, the spectral index distribution in \ass~exhibits two prominent maxima, with the one in the northeast aligning with the peak of $^{13}$CO($3-2$), and the peak in the south coinciding with local maxima in CO($2-1$) and CO($3-2$) azimuthal profiles. These results suggest a correlation between dust grain size and/or optical depth with the gas distributions. Dust-trapping along the spiral arms possibly contributes to the spectral index minima in AS 205 N; however, the observed asymmetry across both disks suggests the involvement of additional mechanisms.

Context. Solar activity, dominated by the 11-year cyclic evolution, has been observed directly since 1610. Before that, indirect cosmogenic proxy data are used to reconstruct it over millennia. Recently, the precision of radiocarbon measurements has improved sufficiently to allow reconstructing solar activity over millennia. Aims. The first detailed reconstruction of solar activity, represented by annual sunspot numbers, is presented for 1-969 CE. Methods. The reconstruction of sunspot numbers from D14C was performed using a physics-based method involving several steps: using a carbon-cycle box model, the 14C production rate, corrected for the geomagnetic shielding, was computed from the measured data; The open solar magnetic flux was computed using a model of the heliospheric cosmic-ray modulation; Sunspot numbers were calculated using a model of the evolution of the Sun's magnetic field. The Markov Chain Monte Carlo approach was used to account for different sources of uncertainty. Results. Annual sunspot numbers were reconstructed for the first millennium CE. This period includes one extreme solar event of 774 CE and one Grand solar minimum of 650-730 CE. We could identify 91 solar cycles, of which 26 were well-defined, while 24 and 41 were reasonably and poorly defined, respectively. The mean cycle length was 10.6 years, but the lengths of individual cycles vary between 8 and 15 years. The existence of empirical Waldmeier's relations remains inconclusive. No significant periodicities were found beyond the 11-year cycle. Conclusions. This work fills the gap in the solar cycle statistics between the previously reconstructed first millennium BCE and the second millennium CE, providing vital constraints for the solar dynamo and irradiance models. A consistent 3-millennium-long reconstruction of sunspot numbers, based on a composite multi-proxy cosmogenic record, is pending.

Inspired by string theory and cosmological constant problem, it is plausible that the Universe's vacuum structure is characterized by a landscape of metastable vacua. The existence of dark matter and dark energy further suggests that the dark sector may inhabit its own "dark landscape". If the dark vacuum is metastable, bubbles of lower-energy phases can nucleate at an approximately constant rate. Because the Hubble expansion rate is monotonically non-increasing with cosmic time, such nucleation can eventually lead to percolation and completion of a dark-sector phase transition. In this work, we investigate the phenomenological consequences of this transition, focusing on the resulting stochastic gravitational-wave background and the potential formation of primordial black holes. We find that the gravitational wave spectrum peaks at $k_{\mathrm{peak}}=3.1 H_{\mathrm{PT}}$, with an amplitude $\Omega_{\mathrm{GW}}^{\mathrm{peak}}\simeq1.5 \Omega_\gamma(\Delta\rho/\rho_{\mathrm{tot}})^2$. Furthermore, the formation of primordial black holes is suppressed due to $\Delta N_{\mathrm{eff}}$ constraint.

A dense neutrino plasma can exhibit collective flavor evolution caused by neutrino--neutrino refraction. Recently, a new class of exact nonlinear inhomogeneous solutions was discovered: single-wave (SW) solutions of the fast flavor system. The key property is that the flavor occupation numbers remain homogeneous, whereas the field of flavor coherence varies spatially with a single wave vector. The equations of motion for this structure resemble those of a collection of classical spins, in analogy with the homogeneous slow and fast flavor cases. In contrast, the SW system is not integrable (it does not possess Gaudin invariants) so that, while two-beam pendulum solutions are inevitable, they do not extend to a multi-angle system. We develop a taxonomy of all known nonlinear collective flavor solutions, explaining the overlap between categories and their differences.

Inspired by regular black holes (RBHs) sourced by dark matter halos, we generalize the anisotropic energy-momentum tensor by relaxing the $P_r = -\rho$ condition between radial pressure and density. We demonstrate that while RBHs are a unique special case, a broader class of relations yields horizonless compact stars. Under specific parameter limits, these objects approach a ``frozen state," mimicking black hole features without an event horizon. These compact star solutions could satisfy weak energy conditions and provide a robust mechanism for dark matter-sourced black hole mimickers.

This paper examines how, in the tenth century, medieval Iranian scholars Abu Rayhan al-Biruni and Abu al-Wafa al-Buzjani determined the difference in geographical longitude between the cities of Baghdad and Khwarezm through simultaneous observation of a lunar eclipse. Brief academic biographies of these scholars are presented, with emphasis on their contributions to mathematics and astronomy. The study discusses the importance of determining geographical coordinates - especially longitude - in the science of the 10th-11th centuries, provides an overview of the methods of coordinate determination available at the time, and highlights the problem of synchronizing remote observations prior to the advent of electronic communication. Particular attention is devoted to a detailed analysis of the method based on observing a lunar eclipse to simultaneously measure longitude differences: the necessary conditions and organization of the experiment, the instruments employed, the mathematical calculations, and error estimates are described. The longitude difference obtained by al-Biruni and al-Buzjani is compared with modern values. The conclusion discusses the scientific significance of this method for the history of science and astronomy.

We study parity-violating birefringence induced by an asymmetric neutrino background in plasmas associated with active galactic nuclei (AGN). We derive a directionality factor arising from the relative bulk motion between the neutrino medium and plasma, and show that it can produce an anomalous frequency dependence of the polarization-rotation angle, distinct from the $\omega^{-2}$ scaling of Faraday rotation. This anomalous scaling can occur either at the resonance plasma frequency condition $\omega \simeq \omega_p$, or when $E_\nu^{0}\simeq m_\nu \omega/\omega_p$ lies within the range of the neutrino energy spectrum. We estimate the effect for three scenarios: jets propagating through the cosmic neutrino background (C$\nu$B), jets with an internal flux of high-energy neutrinos, and accretion-disk plasma permeated by the C$\nu$B. Of the three scenarios, the latter gives the largest rotation angle $\phi_{\rm d} \sim 10^{-35}\,\mathrm{rad}$, at X-ray frequencies. Although the predicted rotation angles are below current polarimetric sensitivity, the identified spectral signatures provide a theoretical framework for probing neutrino asymmetries and AGN plasma properties independent of magnetic field models.

Planetary surfaces are typically analyzed using high-level semantic concepts in natural language, yet vast orbital image archives remain organized at the pixel level. This mismatch limits scalable, open-ended exploration of planetary surfaces. Here we present MarScope, a planetary-scale vision-language framework enabling natural language-driven, label-free mapping of Martian landforms. MarScope aligns planetary images and text in a shared semantic space, trained on over 200,000 curated image-text pairs. This framework transforms global geomorphic mapping on Mars by replacing pre-defined classifications with flexible semantic retrieval, enabling arbitrary user queries across the entire planet in 5 seconds with F1 scores up to 0.978. Applications further show that it extends beyond morphological classification to facilitate process-oriented analysis and similarity-based geomorphological mapping at a planetary scale. MarScope establishes a new paradigm where natural language serves as a direct interface for scientific discovery over massive geospatial datasets.

Compact astrophysical systems such as neutron stars and black holes provide powerful laboratories for testing feebly coupled dark photons (DPs). We investigate light DPs kinetically mixed with the visible photon that need not be the dark matter, focusing on resonant photon-DP oscillations in magnetized, modeled plasma environments. We show that realistic non-monotonic plasma density profiles generically enhance resonant conversion relative to monotonic models, leading to substantially stronger constraints on the photon-DP kinetic mixing parameter ($\epsilon$). Using spectral data from the supermassive black hole (SMBH) M87*, extending to the LOFAR band, we derive a bound $\epsilon \simeq 7\times10^{-6}$ at the DP mass $m_{A'} \simeq 5\times10^{-7}\,\mathrm{eV}$ for oscillation distance $3r_{\rm ph}$, where $r_{\rm ph}$ denotes the photon sphere radius. From the Crab pulsar-wind Nebula, we obtain an even stronger constraint, $\epsilon \simeq 8\times10^{-7}$ at $m_{A'} \simeq 4\times10^{-9}\,\mathrm{eV}$ for oscillation baselines of order $10^{3}\,\mathrm{km}$, surpassing existing astrophysical limits in realistic plasma backgrounds. While laboratory and cosmological bounds remain slightly stronger at comparable masses, observation of compact objects with larger surface magnetic fields and measurements of photon spectra at lower frequencies would enhance the limits on the photon-DP coupling by orders of magnitude.