Stellar masses are a fundamental property to understand models of pre-main sequence evolution, but their values derived from Hertzsprung-Russell (HR) diagrams are strongly model dependent. We benchmark pre-main sequence stellar evolutionary tracks using stellar masses dynamically estimated by fitting a parametric model to ALMA observations of the $^{12}$CO $(J=3-2)$ line transition emitted by the disks orbiting 20 sources in the old ($4-14$ Myr) Upper Scorpius star forming region. We derive stellar masses from HR diagram fitting for ten different stellar evolutionary models, which we then compare with their stellar dynamical masses for comparison in the stellar mass range $0.1-1.3 \> M_\odot$. Models with a moderate-to-low fraction of cold stellar spots ($f=17\%$) most accurately reproduce the dynamical stellar masses ($100\%$ of the targets agree within $\pm1\sigma$). While a higher spot coverage ($f=34\%$) provides similar stellar mass predictions similar to magnetic equipartition models, larger fractions ($f\geq51\%$) significantly disagree with dynamical masses. Magnetic equipartition models overestimate stellar masses up to a factor $\sim20\%$, whereas non-magnetic models underestimate them up to $\sim12\%$. For some models, there is evidence that the stellar mass discrepancies are anticorrelated with dynamical stellar masses. When stellar dynamical mass priors are considered in HR diagram fitting, the median age of a single source can change up to $\sim25\%$, while the median ages inferred across different tracks become consistent, with the age scatter decreasing by $\gtrsim77\%$. These results provide strong empirical constraints for testing and developing evolutionary models of pre-main sequence stars.

We present a JWST analysis of intracluster light (ICL) in XLSSC 122 at z = 1.98, currently the most distant known strong lensing galaxy cluster with an evolved member population. Using deep JWST imaging complemented by HST data and careful control of systematics, we robustly detect diffuse emission extending to several hundred kpc from the brightest cluster galaxy (BCG) down to about 29 mag arcsec^-2. Multi component PSF convolved Sersic modeling separates the surface brightness profiles into three components: a BCG core, a BCG envelope, and an ICL component, with stable Sersic indices across wavelengths. Nearly flat color profiles indicate minimal radial variation in the stellar populations of the BCG envelope and the ICL. The median ICL fraction measured across seven bands is about 17 percent, demonstrating that the buildup of intracluster stars in massive halos was already well underway by z about 2. The ICL fraction peaks near 5000 Angstrom in the rest frame, resembling the behavior observed in dynamically active clusters. We also detect a southern excess of ICL relative to the best fit Sersic model and quantify it using wavelet based modeling, providing additional support that this system is dynamically active. The BCG + ICL light distribution and strong lensing mass map show strong morphological agreement within about 100 kpc. These findings establish the ICL as an early forming and dynamically informative component of massive halos.

We present the first systematic census of quenching mechanisms using kinematic asymmetries in a large sample of $\sim$6,700 galaxies from the MaNGA survey, providing a unified view of what halts star formation in the local Universe ($z<0.2$). We quantify stellar and nebular gas disturbances through the higher-order terms of a Fourier series expansion. These asymmetries serve as powerful diagnostics, as different quenching mechanisms leave distinct kinematic signatures on gas and stars. Our analysis reveals that the most effective quenching pathways leave minimal kinematic imprints by the time galaxies are fully quenched. This "kinematic regularity" points toward slow-acting processes (>3 Gyr) such as starvation and maintenance feedback. A striking finding emerges from our mass-matched analysis: quenched symmetric satellites are significantly more compact than their asymmetric counterparts ($3.4\sigma$), a trend that is even more pronounced for symmetric centrals ($12.3\sigma$). Our results suggest that environment drives the dominant satellite quenching pathway through rapid gas stripping followed by long-term starvation. These compact, kinematically undisturbed satellites (the most representative case within our sample) have undergone intense gas stripping and central compaction, creating bulge-like structures with old, metal-rich stellar populations. Combined with halo gas cut-off and the prevention of cosmological accretion due to starvation, this creates an irreversible quenching path. Conversely, the larger sizes of disturbed, quenched centrals are consistent with merger-driven growth. Internal processes, likely driven by the AGN cycle over 1-3 Gyr that prevents hot halo gas cooling, sustain quenching maintenance in this population. The absence of asymmetric satellites in the star-forming regime suggests environmental quenching operates without significant kinematic perturbation.

Galaxy mergers have long been proposed as a mechanism for funneling gas toward galactic centres, potentially triggering accretion onto supermassive black holes (SMBHs) and igniting active galactic nuclei (AGN). While simulations often support this scenario, observational studies have yielded conflicting results regarding the AGN-merger connection. In this study, we analyze 31 galaxies from cosmological zoom-in simulations spanning redshifts $0.5 < z < 3$. We identify mergers using detailed merger trees based on six-dimensional dark matter particle information and identify AGN activity through SMBH accretion histories. To bridge the gap between simulations and observations, we generate mock JWST-like images and extract non-parametric morphological parameters. Employing a $k$-nearest neighbours (KNN) classifier in a five-dimensional space (four morphological parameters and redshift), we identify mergers in the mock-observed dataset. Our analysis reveals a statistically significant enhancement of AGN activity in merging systems, particularly at lower redshifts ($0.5 < z < 0.9$), where central gas reservoirs are more depleted. This supports the view that mergers contribute more significantly to AGN triggering in environments with low internal gas reservoirs, while their impact may be less pronounced in gas-rich systems. However, when relying solely on morphological classifications from mock observations, the observed AGN-merger connection weakens, especially at higher redshifts. This underscores the challenges in detecting merger-induced AGN activity observationally and highlights the importance of combining simulations with realistic mock observations to fully understand the AGN-merger relationship.

Galaxy clusters offer powerful insights into the large-scale structure of the Universe and the physics of baryons in hot state. Their scientific exploitation, however, hinges on our ability to accurately measure key thermodynamic properties. In this work, we aim to assess the reliability of current analysis techniques in reconstructing these properties, with particular focus on samples similar to those observed in the Cluster HEritage project with XMM-Newton (CHEX-MATE). We develop a suite of dedicated end-to-end simulations of CHEX-MATE-like clusters selected from large scale hydrodynamical simulations, and processed through a newly developed realistic XMM-Newton simulator. We apply a full X-ray data analysis pipeline to the mock datasets, including imaging, spectral fitting, and profile reconstruction. The gas density profiles can be robustly recovered across a wide radial range, when using azimuthal mean surface brightness profiles. Our reconstruction techniques are able to reproduce the intrinsic density profile with the correct scatter, with deviations of at most 10% between 0.1 and 1xR500c. The gas mass is reconstructed with better than 1% accuracy. Accurate measurement of temperature profiles is more challenging and possibly subject to biases, particularly in the presence of azimuthal variations and multi-temperature gas along the line of sight, which dominate over projection effects. Our results highlight the need for caution in interpreting cluster temperature measurements and underscore the value of tailored mock observations for understanding observational systematics. These findings also suggest that biases in X-ray temperature measurements may alter the interpretation of the thermodynamical state of the intra-cluster medium, an outlook particularly relevant in light of recent low velocity measurements from the XRISM mission.

Stage-IV photometric galaxy surveys are designed to measure the position and shapes of billions of galaxies. Their aim is to characterise the large-scale distribution of matter in the Universe using galaxy clustering and weak gravitational lensing. As a byproduct, stage-IV surveys are expected to detect more than a hundred thousand strong gravitational lenses. In this article, we propose the use of weak-lensing perturbations to strong lenses, specifically their line-of-sight (LOS) shear, as a cosmological probe. This new observable allows us to define three new correlation functions: the LOS shear with itself, with galaxy positions, and with galaxy shapes, thereby promoting the standard $3\times 2$pt correlation method to a $6\times 2$pt scheme. We design estimators for these new correlation functions and determine their expectation values as a function of the matter power spectrum. We then derive the analytical expression for the full covariance matrix of the $6\times 2$pt correlation scheme. Considering various scenarios for the stage-IV strong-lensing samples, we demonstrate that the cosmological information carried by the LOS shear of strong lenses will be detectable with a very high signal-to-noise ratio, even in the most pessimistic of cases. Strong lenses are thus extremely promising cosmological probes, whose synergy with galaxy positions and shapes should also contribute to mitigating systematics in stage-IV surveys.

The classical picture that planet formation occurs in protoplanetary disks that are isolated from their environment is undergoing a major shift toward a more connected picture. An increasing amount of evolved disks are found to be actively interacting with their environment, often showing various types of spiral structures. In this work, we aim to investigate if these spirals can be a direct result of ongoing late infall using the grid-based 3D hydrodynamics code FARGO3D. We perform a detailed analysis of the spiral properties and appearance in scattered light and CO line emission using the radiative transfer code RADMC3D. In scattered light, we find both well-defined spirals with few arms (m=2) and more flocculent structures: The gradual accretion of gas remnants after a major accretion event has the most success in the former, whereas active accretion via streamers favors the latter. The m=2 spirals we find have a very low pattern speed, making them easily discernible from spirals caused by a perturber. We also find spiral patterns in the $^{12}$CO residual motions, but their morphology does not match the one found in scattered light. The disk perturbations are strongest in the upper layers (z>4H), which is reflected by the reduced amplitude of the residual motions in the more optically thin $^{13}$CO emission. Moreover, we find that the formation of m=2 spirals is not promoted in disks with lower mass, despite being more susceptible to deeper kinematic perturbations. While the late-infall streamers impact planet formation directly through the delivery of fresh material, we show that the midplane remains unperturbed unless the infalling mass is of the same order of magnitude as the disk mass. Planet formation can therefore only be impacted by late infall through secondary mechanisms that lead to dust trapping or the generation of turbulence starting from surface-level perturbations.

Analyses of baryon acoustic oscillations (BAO) commonly employ template-based methods to extract compressed parameters from the clustering of dark-matter tracers, which are then interpreted in terms of ratios of the sound-horizon scale and cosmological distances relative to a fiducial cosmology. A small mismatch between the sound-horizon scale derived from the standard analytic formulation (integral over the sound speed) and the effective scale imprinted in clustered matter can, however, introduce a systematic bias in cosmological inference. We extend previous work to a broader class of cosmological models, quantify this bias for surveys with DESI-like precision, and propose strategies to correct for the effect. We find that the induced bias becomes a significant fraction of the statistical uncertainty for deviations from the fiducial cosmology, at the level of $|\Delta \Omega_m| = 0.03$ and $|\Delta N_\mathrm{eff}| = 0.3$, and for very precise data corresponding to a forecasted Year-5 DESI survey (or other stage IV dark energy galaxy surveys). We present several ways to correct for this effect, suitable for a variety of applications. We therefore recommend that analyses exploring such parameter regimes either apply the proposed corrections or include an appropriate systematic error budget.

The precision era of multi-messenger astronomy, together with modern astrophysical, cosmological, and gravitational wave observations, increasingly points toward the existence of a ``dark" sector that cannot be explained within the framework of the Standard Model of particle physics and General Relativity. In this review, we explore extensions of standard physics and examine how observational data can be used to probe new particles and interactions. We consider a wide range of scales, from Solar System tests to galactic and cosmological observations, and investigate both conventional dark matter candidates, such as weakly interacting massive particles, and alternative scenarios including ultralight fields and primordial black holes. We discuss constraints derived from compact objects such as neutron stars, black holes, pulsars, and magnetars observations as well as from high-energy astrophysical phenomena. In addition, we analyze extensions of General Relativity involving additional scalar fields and their impact on gravitational wave signals and stochastic backgrounds from primordial black holes. We also study the capture and accumulation of dark matter in compact objects, which can alter properties such as mass, radius, and tidal deformability, and consider scenarios in which dark matter decays into Standard Model particles. While current observations already place significant limits on dark matter and modified-gravity models, upcoming experiments and observatories are expected to further probe or discover such new physics by improving constraints on particle masses and interaction strengths.

The STIS team re-derived on-orbit sensitivity curves for the echelle modes for post-servicing mission 4 observations using the standard DA white dwarf G 191-B2B. These new updates relied on the recent CALSPECv11 models, which introduced improvements in the fluxes of the primary standard stars of the order of~1-3% depending on the wavelength of interest. As part of this effort, the team also released new blaze shift coefficients and echelle ripple tables. We present a detailed description of the procedures followed in the derivation of these new throughputs and the accompanying updates.

Constraining stellar models using asteroseismic and spectroscopic observations is a powerful method for precisely determining the fundamental properties of stars in different kinematic components of our galaxy. We use spectroscopy and individual oscillation mode frequencies to perform a homogeneous modeling study of eight evolved metal-poor stars enhanced in $\alpha$-elements. We compare a full treatment of $\alpha$-enhancement against an ad hoc correction to the total metallicity and show that the stellar properties inferred from asteroseismic modeling using both sets of models agree with each other. Additionally, we find that the uncertainties on stellar parameters derived from the both $\alpha$-enhanced modeling methods are comparable. This is in qualitative disagreement with existing works showing red-giant ages constrained by only the global asteroseismic parameters to depend strongly on the opacities and abundances assumed in 1D modeling. We also show that the observed frequency of maximum oscillation power ($\nu_{\text{max}}$) is larger than the value predicted from applying the $\nu_{\text{max}}$ scaling relation to the masses, radii, and temperatures inferred from the detailed modeling. This discrepancy is pronounced at low metallicities, consistent with recent findings indicating a breakdown of the $\nu_{\text{max}}$ scaling relation for metal-poor stars. Understanding the extent to which the $\nu_{\text{max}}$ scaling relation fails for low-metallicity solar-like oscillators through detailed modeling will enable more accurate mass and age determinations for hundreds of giant stars in the Galactic Halo for which only global asteroseismic parameters are available.

We report the discovery of significant X-ray polarization from the dynamically confirmed black hole X-ray binary (BHXB) GS 1354-64 during its 2025-2026 outburst, obtained with the Imaging X-ray Polarimetry Explorer (IXPE). The observation, obtained shortly after a bright X-ray flare, captures the source in an intermediate state following a stalled (failed) state transition. We discover significant 2-8 keV polarization at the ~4% level with high statistical support--14-sigma significance from frequentist analysis and log Bayes Factor 283+/-1 in a Bayesian framework--measuring PD 4.0+/-0.2% and PA-1+/-2 degrees (90% credible interval). The PD exhibits a statistically significant increasing trend with energy--the strongest such increase yet observed by IXPE in a BHXB--going from 2.1+/-0.3% in the 2-3 keV band to 11+/-3% in the 6.5-7 keV band, while the PA appears stable across both energy and time to within statistical uncertainties. Timing analysis of the IXPE data reveals a ~5 Hz Type-C quasi-periodic oscillation. IXPE + NuSTAR spectropolarimetric modeling suggests that the data can be described by polarized thermal disk and Comptonized components with PAs differing by ~90 degrees, or by a dominant Comptonized polarized component whose effective PD increases across the IXPE bandpass--the inferred component-level polarization levels are therefore model-dependent. In either picture, GS 1354-64 retains a strong coronal component during the transitional period observed by IXPE. These results illustrate how X-ray polarimetry can provide a sensitive diagnostic of the accretion state and geometry in black hole X-ray binary accretion flows, exploring a liminal phase at the cusp of state transition.

In this work, we present comprehensive observational constraints on the time-varying vacuum $\Lambda(t)$CDM cosmology using the latest baryon acoustic oscillation (BAO) data from Dark Energy Spectroscopic Instrument (DESI) Data Release 2 measurements in combination with cosmic chronometer $H(z)$ data. We explicitly quantify the impact of vacuum dynamics on the expansion history, the total effective equation of state parameter of the unified cosmic fluid, the effective dark energy equation of state parameter, and the deceleration parameter in the spatially flat $\Lambda(t)$CDM model. We perform a full Markov Chain Monte Carlo (MCMC) analysis and statistical model comparison, providing a consistent assessment of the $\Lambda(t)$CDM model relative to the standard $\Lambda$CDM framework. Our results demonstrate that $H(z)$ and BAO observations strongly constrain deviations from the $\Lambda$CDM model, driving the vacuum dynamics parameter toward $\alpha \simeq 0$, while significantly reducing parameter degeneracies and alleviating the Hubble tension.

The noise in bias frames for all four readout amplifiers in the Advanced Camera for Surveys (ACS) Wide Field Channel (WFC) is dependent on row number. This is because dark current accumulated during readout increases across the detector, influencing and increasing the read noise as a function of row number. In this report, we investigate bias frames taken with the ACS/WFC to explore the column dependence of read noise for each of the amplifiers for different anneal periods. Analyzing the data, we find that there is no column dependence of read noise and that the read noise values for the physical pre-scans are approximately 0.5 e$^-$ lower than in the science arrays because there is no readout dark accumulated in this area. We further investigate 1) the evolution of read noise over an anneal period, 2) a linear decrease in read noise within the initial columns per amplifier, and 3) pixels in elevated read noise columns. We conclude that 1) there is no visual trend of read noise over an anneal period, 2) amplifiers A and C have an initial linear decrease of read noise in the science arrays, and 3) masking unstable hot pixels in a column will decrease its read noise values.

We present the spatial part of the point source signal extraction strategy for the upcoming CHORD galaxy survey. CHORD, the Canadian Hydrogen Observatory and Radio-transient Detector, is an under-construction drift-scanning compact interferometric radio telescope. CHORD comprises 512 six meter dishes and observes in the 300 to 1500 MHz frequency range. One of its science goals is producing a catalogue of galaxies detected by the neutral hydrogen (HI) 21 cm emission line. CHORD's highly redundant dish layout creates the problem of spatial aliasing, the effect where the same signal could be feasibly produced from sources at multiple locations on the sky. The search will be done with a matched filter in the visibility plane. This paper presents the search strategy and a prediction tool that can quickly estimate the matched filter response at a given sky position, allowing a prediction of alias locations and severity. This tool confirms that although aliases are impossible to distinguish in a single snapshot, they become possible to distinguish when combining data over a period of time. It predicts that aliases will be harder to distinguish for observations closer to the celestial equator, but that scanning with offset adjacent strips can remove this degeneracy. It predicts that the optimal strategy for a single offset to disambiguate aliases is to re-point the array in declination by about two degrees. A future paper will combine these findings with realistic noise estimates and galaxy population statistics to make forecasts of the population of galaxies that CHORD will detect.

We test whether parsec-scale stellar systems in the Milky Way follow the galactic radial acceleration relation (RAR) or the baryonic TullyFisher relation (BTFR). We analyse 5646 Gaia DR3 open clusters from the Hunt \& Reffert catalogue. Observed accelerations are derived from velocity dispersions and characteristic radii, and baryonic accelerations from stellar masses and characterisitc radii. The clusters are placed on the RAR and BTFR planes and compared with Newtonian and MOND expectations. Approximately 90 per cent of open clusters (those with $N_\star \leq 250$) lie close to the RAR, albeit with significant scatter. In a first-of-its-kind test, a smaller fiducial sample is consistent with a best-fitting acceleration scale $g_\dagger \approx 1.2 \times 10^{-10}\ \mathrm{m\,s^{-2}} \pm 0.5$ dex, compatible with canonical MOND values. More massive clusters approach the Newtonian virial expectation. No correlations are found between RAR residuals and galactocentric radii, distance to the Galactic disk midplane, age, or morphology. Tidal effects and unresolved binaries are insufficient to reproduce the observations without fine-tuning. Interpreted within a MOND framework, the alignment of most open clusters with the RAR and BTFR suggests that low-acceleration dynamics operate on parsec scales within the Milky Way. This implies that the Galactic gravitational field is not smooth on these scales and may include regions where the total gravitational acceleration falls below $a_0$, partially mitigating the external field effect, thereby motivating higher-resolution modelling of the Galactic potential and informing other small-scale gravity tests within the Galaxy.

We present a systematic analysis of 2544 cataclysmic variable systems and related candidates observed by the Transiting Exoplanet Survey Satellite (TESS), with the aim of compiling a comprehensive catalogue of orbital periods. Using 2-minute photometric time-series data, we applied an automated algorithm to generate Lomb-Scargle periodograms and identify the most significant coherent periodic signals, which were subsequently verified through visual inspection. This process yielded a confident sample of 910 sources exhibiting at least one periodic signal, hereafter referred to as the Cataclysmic Variable Confident Catalogue (CCC). For each object, we report the most likely orbital period together with additional periodic features such as spin modulations and/or superhump signals when present. To assess consistency with previously published measurements, we cross-matched the CCC with the Ritter \& Kolb catalogue, identifying 300 overlapping systems, of which 215 showed full agreement with the R\&K orbital periods, while 39 displayed discrepancies for which the R\&K values were revised based on our TESS measurements and supporting evidence from the literature. Overall, the CCC provides a means to validate known orbital periods, propose corrections where necessary, and offer new determinations for systems with previously unknown periods, thereby supporting a more complete and reliable characterisation of the cataclysmic variable population.

Thermal imaging of debris disks is resolving the vertical height in an increasing number of systems, enabling the use of the vertical structure to decode the dynamical state of the planetary system. In this study, we examine the multi-wavelength structure of the archetypical edge-on debris disk of $\beta$ Pic, extensive imaging of which across mid-infrared to millimeter wavelengths makes it the prime system to study the vertical height across different grain size populations. We non-parametrically modelled the radial profiles and constrained the vertical height at each wavelength while taking into account the vertical warping, finding the disk to be on average 1.5 times thicker vertically in the mid-infrared compared to the millimeter and the scale height to be relatively constant across radius. The decreasing scale height with wavelength is in contrast to predictions from collisional damping, and could be a result of the combined effect of radiation pressure and random collisions. We also show that the disk is warped at millimeter wavelengths and find tentative evidence for clumps in ALMA images which will require follow-up observations to confirm. The millimeter vertical warping is consistent with findings in scattered light and the secular perturbation interpretation due to the inner giant planets, which could also explain the relatively constant apparent scale height across radius, and potentially earlier findings of a non-Gaussian vertical profile which this study confirms.

A recently proposed generalized entropy by Luciano and Saridakis extends the standard Boltzmann-Gibbs and Bekenstein-Hawking framework through a microscopically motivated construction involving two independent entropic exponents. When applied within the gravity-thermodynamics correspondence, this entropy leads to a modified cosmological dynamics that can be interpreted as an effective dark energy sector of entropic origin, while recovering $\Lambda$CDM in appropriate limits. In this work, we perform the first observational confrontation of the resulting entropic cosmology at the background level. Focusing on the case $\alpha_\delta=0$, we constrain the model using Cosmic Chronometers, Pantheon$^+$ Type Ia supernovae calibrated with SH0ES, BAO measurements from DESI DR2 and compressed Planck 2018 CMB information. We find that the model yields a statistically robust fit to the combined data sets and can simultaneously satisfy Pantheon$^+$, SH0ES and CMB shift-parameter constraints, unlike $\Lambda$CDM. Although the entropic parameters remain close to their standard values, the $\Lambda$CDM limit is excluded at the $2\sigma$ level within the restricted parameter space considered. These results indicate that the Luciano-Saridakis entropic cosmology offers a viable extension of the standard model with the potential to alleviate the Hubble tension at the background level.

We report the discovery of a Lyman-alpha emitter (LAE) protocluster at z = 4.90 in the COSMOS field, comprising four distinct overdensity peaks spanning ~65 x 36 cMpc^2, with the primary concentration exhibiting a 4-fold surface density enhancement relative to the field within a 1.5 proper Mpc (pMpc) radius. Using SILVERRUSH narrowband survey data combined with JWST COSMOS-Web imaging, we perform a first systematic rest-frame optical and UV morphological comparison of protocluster versus field LAEs at this redshift using JWST NIRCam rest-frame UV (F150W, ~2540 Angstrom) and optical (F277W, ~4700 Angstrom) imaging. Sersic profile fitting for 16 protocluster members and 23 field LAEs reveals a significant size difference: protocluster LAEs are ~40% larger in rest-optical (median R_e = 0.81+0.26/-0.04 kpc vs. 0.58+0.11/-0.04 kpc, p = 0.041) with no significant difference in rest-UV (p = 0.51) or Sersic index. At fixed stellar mass, protocluster LAEs are offset by +0.12 dex (~31%) in rest-optical size from the field size-mass relation (68% CI: [+0.08, +0.21]; Mann-Whitney p = 0.033), with 75% exhibiting positive size residuals compared to 44% of field LAEs. This wavelength-dependent environmental signature suggests that protocluster environments at z ~ 5 preferentially affect extended stellar populations, possibly through tidal interactions or an earlier onset of star formation in the dense environment, with no significant environmental difference detected in rest-UV sizes, providing observational evidence for environmental influences on the structure of LAEs during the early build-up phase of cosmic star formation.

Accretion is the primary driver of protostellar evolution, regulating mass assembly and shaping the physical and chemical environments of young stellar objects. Quantifying accretion in the Class 0 protostellar phase is particularly important, yet remains observationally challenging due to high extinction toward the central protostars. In this paper, we present JWST NIRSpec and MIRI/MRS IFU data towards the Class 0 protostar L1527 IRS. We extract one-dimensional spectra and find emission from atomic and molecular hydrogen, water, OH, and several ionic species. The atomic hydrogen lines, Br$\alpha$, Pf$\alpha$, and Pf$\gamma$ are the most critical to this study since they can be used as accretion diagnostics. The existence of these atomic hydrogen lines viewed in scattered light indicates that accretion is likely occurring magnetospherically rather than through a boundary layer. Moment 0 emission maps show that the hydrogen emission is co-spatial with the scattered light continuum with a strong east-west asymmetry which is not due to outflow shocks. We additionally present moment 0 maps of other detected species and discuss their emission morphology. By primarily analyzing the Br$\alpha$ line, the strongest of our detected atomic hydrogen lines, we characterize the accretion onto L1527 IRS by estimating the accretion luminosity to be $0.4~\text{L}_\odot$ and the accretion rate to be around $1\times10^{-7}~ \text{M}_\odot \text{yr}^{-1}$. We lastly discuss the implications of our results with respect to both non-steady and asymmetric accretion possibly occurring in L1527 IRS.

Modern datasets have the precision necessary to uncover new information by including higher-order, non-Gaussian information into cosmological inference. The marked power spectrum offers access to such information while preserving the structure of two-point correlators. This approach to higher-order statistics has the advantage that many modeling questions can directly benefit from progress already made in standard cosmological analyses using the power spectrum and correlation function, while increasing the data vector size negligibly and retaining much of the degeneracy-breaking power of the bispectrum. In this work, we first restructure the marked power spectrum to isolate its higher-order information and demonstrate its ability to break parameter degeneracies. We then investigate the effect of survey geometry on the marked power spectrum and find that a treatment similar to that of the power spectrum is sufficient. Additionally, we investigate the perturbative modeling and covariance structure of the marked power spectrum, shedding light on its degeneracy breaking power and cross-covariance with the power spectrum. Finally, we demonstrate that the cosmology dependence of the marked power spectrum is smooth, indicating that cosmological inference is possible by modeling the cosmology dependence through interpolation rather than analytical modeling.

Cosmological reionization was a highly out-of-equilibrium event that affected every parcel of the intergalactic medium, making it a candidate for astrophysical generation of intergalactic magnetic fields. During reionization, the first stars and galaxies ionized the surrounding, largely neutral, medium in ever-expanding envelopes. Photoionization from sources on one side of the front, combined with the quadrupolar angular dependence of the photoionization cross section, leads to an anisotropic electron velocity distribution. We investigate instabilities in these reionization fronts as a mechanism to generate seed magnetic fields. The Weibel instability has the potential to create a magnetic field from these anisotropies. We calculate the magnitude of the isotropic and anisotropic distribution within a simulated reionization front. We find that the fractional anisotropy can grow to $6\times 10^{-3}$ toward the middle of the ionization front. We show that the linear growth timescale of the Weibel instability is fast compared to the crossing time of the ionization front ($\sim 2\times 10^5$ seconds). We briefly speculate on the possible non-linear evolution of the instability and the implications for cosmological magnetogenesis.

M dwarfs are prime targets for exoplanet searches due to their low masses and radii, which enable the detection of small planets in their habitable zones (HZs). However, the magnetic activity of M dwarfs can introduce signals in radial velocity measure- ments that may be mistaken for planetary signatures, making the understanding of stellar activity cycles crucial for accurate planet detection and characterisation. We aim to identify and characterise long-term magnetic activity cycles in M dwarfs using a homogeneous and extensive spectroscopic dataset in order to better understand their magnetic variability and its implications for exoplanet detection. We analysed 13 years of high-resolution spectra obtained with the SOPHIE spectrograph for two early M dwarfs known to host exoplanets. We simultaneously monitored chromospheric activity using two indicators, the H{\alpha} index and the Mount Wilson S -index. Long-term trends were modelled using both sinusoidal and low-order polynomial fits to robustly identify stellar activity cycles. For GJ 617A, we report a cycle of approximately 4.8 years, while for GJ 411, we find several characteristic timescales of variability of about 4.9 years. In addition, TESS photometric data reveal signs of short-term variability in GJ617A. The periods of the long-term variability detected for GJ 617A and GJ 411 do not coincide with any of the planetary signals previously reported, which reinforces the hypothesis that they are of magnetic origin. If indeed the variability is due to activity, the cycles detected would not be driven by the same mechanism: The cycle in GJ 617A is consistent with a solar-like dynamo, while the rotation seems to play a different role in the long-term cycles detected in GJ 411.

Neutrino astronomy provides another window to exploring the Universe, exemplified by the detection of a megaelectronvolt neutrino burst from the core-collapse supernova (CCSN) SN~1987A (refs.~\citenum{hir+87,bio+87}). Commonly discussed theories suggest that some CCSNe could produce neutrinos with energies a thousand times more than those of SN~1987A \cite{tm18}, which has been probed with new-generation facilities \cite{abb+12,aar+15,abb+23}. The interaction of SN ejecta with a dense circumstellar medium (CSM) or a jet, launched in a CCSN, being choked in the stellar envelope of the progenitor or an outside CSM are both well-accepted scenarios for the high-energy neutrino production. Here we report the detection of a high-energy neutrino flare at a 3.9$\sigma$ significance from SN~2017hcd, made by our analysis of the public track-like neutrino data taken by the IceCube Neutrino Observatory \cite{IceCube17}. A Type IIn SN with optical emissions arising from the ejecta--CSM interaction, SN~2017hcd's neutrino flare lasted $\sim$1--2 month, with its central time $\sim$14-day prior to the SN's optical discovery time. Its estimated isotropic neutrino energy (all flavors) is approximately two orders of magnitude higher than the energy ($\sim 10^{50}$\,erg) carried in the SN's ejecta, too high to be explained with the ejecta--CSM scenario. Thus, a choked jet may be the source of the neutrino flare.

Based on the Ca\,II H and K lines observed by LAMOST, we employ the photospheric ($R'_{\rm HK}$) and basal ($R^+_{\rm HK,L}$) flux calibrated chromospheric activity indices to examine the relationship between chromospheric activity and the stellar rotation rate. We identify the rotation periods of 11,108 stars observed by Kepler and TESS by cross-matching our chromospheric activity catalog with previous studies. Our statistical results show that chromospheric activity increases with the rotation rate until it reaches a saturation level. As the stellar effective temperature increases from 4950 to 5850 K, the saturation values of the rotation period ($P_{\rm rot}$) vary correspondingly from 4.38 to 1.23 days for $R'_{\rm HK}$ and from 9.88 to 1.33 days for $R^+_{\rm HK,L}$. Similarly, the corresponding saturation Rossby number Ro ranges from 0.200 to 0.032 for $R'_{\rm HK}$ and from 0.302 to 0.107 for $R^+_{\rm HK,L}$. The saturation is also found to be significant in stars with thick convective zones, whereas it is less apparent in stars with higher effective temperatures. For solar-like stars in the $T_{\rm eff}$ range of 4800 to 6000 K, The values of chromospheric activity indicators are saturated when $P_{\rm rot}<1.45 $ days (Ro$<$0.100) and $P_{\rm rot}<2.85 $ days (Ro$<$0.097) for $R'_{\rm HK}$ and $R^+_{\rm HK,L}$, respectively.

Determining precise stellar ages and masses for evolved giants is crucial for Galactic archaeology but challenged by spectral degeneracies. Gaia's low-resolution XP spectra offer a unique opportunity to infer these parameters on a massive scale using data-driven methods. We extend a transformer-based astronomical foundation model to evolved stars, establishing a unified framework to simultaneously predict atmospheric parameters ($T_{\mathrm{eff}}$, $\log g$, $[\mathrm{M}/\mathrm{H}]$) and evolutionary labels (mass, age) with physical consistency. Treating spectra as token sequences, we integrated mass and age into the model's vocabulary. The model is trained on Gaia XP spectra cross-matched with the APOGEE DR17 DistMass catalog. Our generative approach enables flexible input handling, including spectral inpainting and parameter-to-spectrum generation. On an independent test set, the model achieves a prediction scatter of $\sigma \approx 0.114 \, M_{\odot}$ for mass and $\sigma \approx 1.334$ Gyr for age. Beyond numerical accuracy, it successfully reproduces the giant branch's mass-luminosity relation and autonomously disentangles interstellar extinction from intrinsic temperature variations without explicit physical priors. It also robustly recovers missing spectral data and estimates reliable uncertainties. Validating that foundation models can internalize stellar physics from data, this physically-aware, probabilistic framework offers a powerful tool for unraveling Milky Way history using large-scale spectroscopic surveys.

A machine-learning-based method is developed to identify objects with unusual stellar spectra. The method employs an autoencoder, a neural network trained to compress spectral data into a low-dimensional representation and subsequently reconstruct it. Spectra that deviate significantly from the dominant patterns in the training dataset are identified using the reconstruction error as an anomaly score. The models are applied to selected datasets from the MaNGA Stellar Library, an empirical library of stellar spectra. Several spectra are flagged as anomalous: an object with likely instrumental and/or reduction issues, two carbon stars, and an oxygen-rich thermally pulsating asymptotic giant branch star. The sources of the large reconstruction errors are examined, and the effectiveness and limitations of autoencoder-based approaches for detecting anomalous stellar spectra are discussed.

We investigate the accretion dynamics of the black hole X-ray binary Swift J1727.8-1613 during its $2023-2024$ discovery outburst that lasted for $\sim10$ months. Insight-HXMT monitored the rising phase of the outburst of Swift J1727.8-1613 roughly continuously from 2023 Aug 25 to 2023 Oct 05. Strong signatures of type-C Quasi-Periodic Oscillations (QPOs) are observed during this phase of the outburst. In our recent paper, nature of the QPOs are studied with the propagating oscillatory shock (POS) model. In this paper, we report on the observation of both positive (or hard) and negative (or soft) time-lags in the $4-10$ keV (LE), $10-30$ keV (ME), and $30 -150$ keV (HE) bands, computed with respect to the $2-4$ keV reference band. We detect a clear transition from hard to soft lags as the outburst evolves. We show the evolution of QPOs and associated time-lags between different X-ray energy bands, correlated with changes in the QPO frequency, spectral state, and the size of the Comptonizing region. Our analysis reveals strong anti-correlations between the time-lags and both QPO frequency and photon index, and a strong positive correlation with the shock location. These evolving lag characteristics and their correlations provide crucial insights into the changing accretion geometry and the interplay of radiative processes, further supporting dynamic models like the POS in explaining the coupled spectro-temporal evolution in black hole X-ray binaries.

Clouds are the largest source of uncertainty in climate simulations. For exoplanets, cloud simulation is particularly challenging because of the lack of observational data to tune parameterized cloud models. Here we apply Community Aerosol and Radiation Model for Atmospheres (CARMA), a size-resolved bin cloud microphysics model, to the atmospheric global climate model Community Atmosphere Model (CAM6) and simulate exoplanets with a range of planetary rotation rates. CARMA produces fewer liquid clouds than the native CAM6 parameterized cloud microphysics scheme (Morrison-Gettelman two-moment microphysics, MG), more ice clouds, and a significantly different ice cloud size distribution. Overall, this leads to a decrease in the magnitude of the net CRE by 4-10 $W/m^2$, which is unlikely to change the determination of habitability from a climate perspective in most cases. The difference in ice cloud size distribution is likely to strongly affect transmission spectral retrievals. Our work confirms that the MG parameterized cloud microphysics scheme can produce reasonable climate simulation when extrapolated to some exoplanet contexts and highlights the value of resolved cloud microphysics for evaluating parameterized schemes and for interpreting observations.

We measure galaxy stellar mass functions (SMFs) for nine of the most massive galaxy clusters in the local universe ($0.07 < z < 0.11$) using deep and complete spectroscopy from the MAssive Cluster Survey with Hectospec (MACH). We construct the cluster SMFs down to $\log(M_*/M_\odot) \gtrsim 8.5$. For comparison, we measure the SMF for field galaxies, complete to $\log(M_*/M_\odot) \gtrsim 10.5$, based on Sloan Digital Sky Survey (SDSS) spectroscopy over the same redshift range. The mean MACH SMF shows a shape similar to that of the field SMF but with a significantly higher amplitude at $\log(M_*/M_\odot) < 11.4$. At $\log(M_*/M_\odot) > 11.4$, the MACH SMF shows a clear excess, indicating the contribution of massive galaxies, including Brightest Cluster Galaxies (BCGs). Based on homogeneous MACH spectroscopy, we compare SMF shapes for quiescent and star-forming members as a function of cluster-centric distance. The quiescent SMFs display a curved shape with a peak at $\log(M_*/M_\odot) \approx 10.5$; the star-forming SMFs decline monotonically with increasing stellar mass. We further compare the mean MACH SMF with SMFs derived from similarly massive clusters in the IllustrisTNG-300 simulations. The shape of the observed and simulated SMFs agree well overall. However, the MACH clusters contain roughly a factor of two more galaxies at $9.0 < \log(M_*/M_\odot) < 10.5$. These results demonstrate that constructing cluster SMFs from complete spectroscopic samples can test simulations and provide powerful constraints on galaxy formation and evolution in dense environments.

We study the impact of supermassive black hole (SMBH) growth, $\langle \dot{M}_\mathrm{SMBH}\rangle$, major and minor galaxy mergers, and gas processes, on the average gas metallicity of galaxies, with the aim to uncover which of these processes drive the scatter in the gas metallicity-stellar mass relation (MZR) at different redshifts in nodes, filaments and voids. At $z=5$, minor mergers produce the largest differential in $\log[Z_g/Z_\odot]$ for all environments, where the node population displays a maximum $0.38$ dex increase in the average $\log[Z_g/Z_\odot]$ compared to non-merging galaxies. The node population also displays a consistent $0.1$ dex reduction in $\delta \log[Z_g/Z_{\odot}]$ across all redshifts, whilst filament and void galaxies show a lower magnitude of reduction. Major mergers show little influence on these same properties. This suggests minor mergers regulate metallicity and contribute to over galaxy mass growth concurrently, accelerating chemical evolution post merger. Between $z=1-3$, a high $\langle \dot{M}_\mathrm{SMBH}\rangle$ leads to a reduction in $\delta \log[Z_g/Z_{\odot}]$ for all environments. Here, node galaxies show the largest reduction of approximately $0.25$ dex, suggesting that metal-rich outflows strongly drive the MZR at intermediate times. Finally, galaxies with low $M_{gas}/{M_{tot}}$ show increased $\delta \log[Z_g/Z_{\odot}]$ across all redshifts and environments, again a $0.25$ dex maximum for node galaxies. These galaxies also spike in $\delta \log[Z_g/Z_{\odot}]$ at late times, below $z=1$. At this time, galaxies in the nodes show negative $\langle \dot{M}_\mathrm{gas} \rangle$ whilst also showing the largest $\delta \log[Z_g/Z_{\odot}]$ values we observe of $0.2$ dex, suggesting the importance of the balance between gas accretion and starvation in driving MZR scatter at low redshifts.

We present a timing analysis of the high-cadence optical light curve of the high-redshift flat-spectrum radio quasar PKS\,0805$-$07 obtained during \textit{TESS} Sector~34 (MJD $59230.90$--$59239.90$). We search for short-timescale quasi-periodic oscillations (QPOs) using complementary time-series techniques, including the Lomb--Scargle periodogram (LSP) and the weighted wavelet $Z$-transform (WWZ), and evaluate their significance against red-noise variability using Monte Carlo simulations. The LSP reveals a dominant modulation at $f \approx 0.597\,\mathrm{d^{-1}}$ ($P \approx 1.7\,\mathrm{d}$) exceeding the $99.99\%$ confidence level, while the WWZ independently recovers a consistent timescale at the $\sim 99.9\%$ level and shows that the signal is temporally localized rather than persistent across the full light curve. The modulation spans $\sim 5$ coherent cycles, indicating a transient quasi-periodic feature. We discuss possible physical interpretations of the detected modulation. In a disk-based scenario, orbital motion of a hotspot near the innermost stable circular orbit implies a black hole mass of $M_{\rm BH}\sim7.2\times10^{8}\,M_\odot$, consistent with typical FSRQ values. Alternatively, magnetohydrodynamic kink instabilities in the relativistic jet can naturally produce day-scale variability for standard blazar parameters and account for the transient character of the signal. We conclude that the observed modulation is consistent with a compact, short-lived structure embedded within stochastic jet variability.

Sagittarius A* (Sgr A*) is a compact radio source at the Galactic center. Observations have confirmed that its mass is approximately (4.1)*10$^{6}$ M$_{\odot}$, and Sgr A* is generally believed to be powered by gas accretion onto a supermassive black hole. Multifrequency radio observations of the pulsar J1745-2900, about 0.12 pc away from Sgr A*, reveal an unusually large Faraday rotation. Combined with X-ray observations, this indicates that there is a strong magnetic field (greater than 8 mG) leading to a low $\beta$ plasma at large this http URL show that the gas starts to be captured by the black hole below tens of thousands of the Schwarzschild radii $r_S$, where the gas pressure starts to dominate. Assuming that the accretion rate along magnetic fields at large scales decreases with the distance to the black hole following a power law, it is shown that, with an accretion disk below tens of $r_S$, as revealed with the EHT observations, there should be a supersonic wind above such a small accretion disk, and the accretion flow may be convection-dominated from tens of $r_S$ to tens of thousands of $r_S$. Detailed modeling is warranted.

Ridges in galaxy density fields measured by photometric surveys are 2D projections of filaments in the cosmic web, and so should lens light from background galaxies. We report on a detection of this effect in Dark Energy Survey Year 3 data at high significance, though not independently of galaxy-galaxy lensing. We describe improvements to the existing subspace-constrained mean shift algorithm to locate these ridges efficiently at scale, and examine the dependence of the signal in simulations on cosmological and algorithmic parameters. We find that it depends primarily on $S_8=\sigma_8 \left( \Omega_m / 0.3 \right)^{1/2}$, and discuss improvements to our methodology that would be needed to allow precision parameter estimation.

Galaxy morphology encodes key information about formation and evolution. Large imaging surveys require automated, reproducible methods beyond visual inspection. Non--parametric indices provide an useful framework, but their performance must be quantitatively assessed. We present a homogeneous catalog of non--parametric morphological indices for DECaLS galaxies with effective radii larger than 2 arcsec. Our goal is to evaluate the reliability of indices in separating spirals and ellipticals, test their consistency with existing classification schemes, and establish their applicability for the upcoming surveys focused in the southern hemisphere. We developed galmex, a modular Python package for preprocessing images and measuring a variety of non--parametric indices. Using bona-fide spirals and ellipticals as control samples, we assessed the discriminatory power of each index, and compared them with CNN-based T-Types and Galaxy Zoo DECaLS labels. We use the indices as input for a Light Gradient Boosted Machine (LightGBM) to obtain probabilistic classifications. Concentration is the most reliable parameter from the Concentratiom + Asymmetry + Smoothness system (CAS), while asymmetry--based indices (A and S) are limited to detecting disturbed morphologies. MEGG indices (M20, Entropy, Gini, G2) provide stronger separation and trace a gradient with T--Type. By using a simple binary (0/1) label for ellipticals/spirals, classifiers trained on non--parametric indices achieve high accuracy and well--calibrated probabilities, dominated by entropy, concentration, and Gini. We release the first public catalog of CA[A_S]S+MEGG indices for DECaLS, together with galmex. We combine the non-parametric indices with machine learning framework to derive spiral/elliptical separation for galaxies below z~0.15 through a probabilistic approach.

This work presents a description of the scientific goals and objectives of OPTIMus (OPTical, Infrared, Millimeter survey of massive star-forming regions), a survey of massive star-forming regions in the optical, infrared, and millimeter wavelengths. The survey is aimed at constructing a comprehensive characterization of the multicomponent and structurally complex interstellar medium in the vicinity of young massive stars, combining both observational and theoretical aspects. Using multi-wavelength observational data, we will reconstruct the three-dimensional structure and determine the physical parameters of HII regions, photodissociation regions, and the surrounding molecular clouds. The paper describes the observational data obtained with the BTA 6-m and Zeiss-1000 telescopes of the Special Astrophysical Observatory of the Russian Academy of Sciences, the 2.5-m telescope of the Caucasian Mountain Observatory of the Sternberg Astronomical Institute of Moscow State University, and the 20-m telescope of the Onsala Space Observatory.

Gas accretion process can fuel both star formation and black hole activity, playing a critical role in galaxy evolution. The counter-rotating structures are believed to originate from gas accretion, serving as an ideal laboratory for studying its impact on galaxy evolution. Based on the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, we built a sample of 147 galaxies with counter-rotating stellar disks (CRDs). This is the largest CRD sample to date, accounting for $\sim$1.5% of the MaNGA survey. For a subset of 138 CRDs, global stellar mass ($M_\ast$) and star formation rate (SFR) were measured in reference. We constructed a control sample with similar $M_\ast$ and SFR but lacking counter-rotating structures. The CRDs relatively exhibit more bulge-dominated morphology, lower molecular gas mass fraction and reside in less dense environment, supporting the hypothesis that they primarily originate from gas accretion. We classified 96 out of 138 CRDs into four types based on their stellar and gas kinematics following the criteria from Bao et al. (2022). There are two additional CRD types: 8 CRDs show misalignment between both stellar disks and gas disk, indicating multiple gas accretion events with differing angular momentum directions; 34 CRDs lack ionized gas emission, showing the highest $M_\ast$ among all the CRD types, which may represent a final stage of CRD evolution. We compared the radial gradients of gas-phase metallicity and stellar population properties between CRD types, and found that the impact of gas accretion on galaxy evolution primarily depends on the abundance of pre-existing gas in progenitors.

We present a novel gradient regularization to completely eliminate the magnetic divergence error in meshless magnetohydrodynamics (MHD), which offers a high spatial resolution and conservative advantage, due to its Lagrangian nature. Comparing with the counterpart of constrained-gradient (CG) technique, we reform $\nabla \cdot \mathbf{B}=0$ by an implicit projection method to modify the magnetic-field gradients. The accuracy of modified-gradient (MG) method is verified and it achieves exact divergence-free results with round-off precision, by using tests of shock tube, 2D and 3D vortex, magneto-rotational instability, and especially, advection experiment, compared with CG method and the GIZMO code. It leads to noticeable improvement in pattern, amplitude and numerical dissipation of divergence error of magnetic field.

At the end of the Cassini mission, Saturn's rings have been claimed to be spectacularly young compared to the age of the Solar System: their unusual ice-rich composition corresponds to initially pure ice rings polluted by interplanetary dust particles for 100 to 400 Myr. Since then, this exposure age has been commonly accepted as the real age of the rings. In this paper, we review the processes that are involved in determining the exposure age. We aim to see how the exposure age depends on various parameters and how relevant it is to define the real rings age. First, a new expression for the gravitational focusing onto planar rings, important parameter but crudely defined in the literature, is derived. Then, an analytical formula describing how the dust fraction varies with time in static or viscously evolving rings is provided, including possible vaporisation at impact. Finally, we introduce a cleaning process from space weathering to possibly alter dust and reduce its amount to make rings look younger than they are. We first found that the gravitational focusing is 5 times less important than previously thought, which automatically increases the exposure age from 0.5 to 2 Gyr. Moreover, depending on the impact properties (vaporisation rate, space weathering efficiency), several billion years can easily be reached. Finally, we find that the dust fraction in the rings converges towards a finite value which, in particular with an efficient space weathering mechanism, can be close to the observed one in the current rings. In this case, neither the age nor the initial composition of the rings can be derived, and the measure of the dust fraction and bombardment rate only constrains the physical parameters of the impacts and the efficiency of the space weathering. As long as the latter parameters are not known, the exposure age argument in favour of young rings is completely undercut.

The New Robotic Telescope (NRT) is an international collaboration to build and operate a 4 m diameter fully robotic telescope. The telescope will take advantage of the superb atmospheric conditions at the Observatory of the Roque de los Muchachos (ORM). In conjunction with a large aperture, entirely robotic operation, quick response, and a set of versatile instrumentation in the optical and near-infrared this guarantees a high scientific impact focused mainly in the area of time domain astronomy. This paper presents the scientific motivation and the status of the project, discussing possible technical solutions under evaluation for the optics, mechanics and control system.

The supernova remnant (SNR) G315.4$-$2.3 (MSH 14$-$63 or RCW 86) exhibits strong emission across the electromagnetic spectrum. Radio polarization observations probe magnetic fields and will help to understand the evolution of the SNR. We aim to investigate the radio spectrum and magnetic field properties of the SNR. We observed G315.4$-$2.3 using the Australia Telescope Compact Array (ATCA), covering the frequency range of 1.1-3.1 GHz. We performed rotation measure (RM) synthesis on the $Q$ and $U$ frequency cubes to obtain polarized intensity and RM. The regular component of the line-of-sight magnetic field was estimated from RM. The fractional polarization versus wavelength squared was used to constrain the properties of the turbulent magnetic field. We obtained image cubes of Stokes $I$, $Q$, and $U$, along with images of polarized intensity $P$, RM, and fractional polarization $p$. The radio spectra are very similar for different areas of the SNR. The foreground RM was estimated to be 55 rad m$^{-2}$, and the internal RM of most SNR areas is less than about 50 rad m$^{-2}$. The regular magnetic field along the line of sight was estimated to be about 1.4 $\mu$G in the southwest, much smaller than the total magnetic field. For most parts of the southwest and northeast, $p$ is less than 8% and is nearly constant with $\lambda^2$. We estimated the ratio of turbulent to regular magnetic field to be larger than about 3. The scale of the turbulent magnetic field for some area in the northwest might be smaller than about 0.4 pc. The radio characteristics, including spectrum and turbulent magnetic field, are very similar in the northeast and southwest, even though the evolution is quite different for these two regions based on the current models. These should be taken into account for future modeling of the evolution of the SNR.

We explore what fraction of delta sunspots in which the polarity inversion line (PIL) is sharp in photospheric magnetograms are made from a writhe kink in an emerging twisted flux rope. We searched simultaneous full-disk magnetograms and continuum images from Helioseismic and Magnetic Imager (HMI) on Solar Dynamics Observatory (SDO) to find 28 random sharp-PIL delta sunspots that are born well on the disk. Only one of these is made from a single newly emerged bipolar magnetic region (BMR) and therefore is a candidate for being made from a single emerging writhe-kinked flux rope. That outcome indicates that few, if any, sharp-PIL delta sunspots are made by a single emerging writhe-kinked flux rope. That is the main new finding of this paper. Each of the other 27 is made by merging of two or more emerging or emerged BMRs. We name delta-sunspot genesis from a single BMR Type I genesis. We identify another three genesis types among the other 27 delta sunspots: Type II, Type III, and Type IV. We present an observed example genesis for each of the four genesis types, and for each example present schematic drawings depicting our scenario(s) for the cause of that example genesis. The core idea of each scenario is that the delta sunspot is made by packing opposite-polarity magnetic flux together by advection into a convection downflow.

NRT is an international project to build and operate the world's largest robotic telescope. The telescope will have a segmented primary mirror with an equivalent diameter of 4 m, a set of simultaneously mounted optical and near-infrared instruments, and a response time of less than 30 seconds. The project builds on the experience gained with the successful twenty-year operation of the Liverpool telescope, and with the GTC optics and control system. All of the above together with the excellent conditions for astronomical observation of La Palma, represents a solid base and guarantees that NRT will be one of the leading facilities in the field of time domain astronomy. This contribution will analyze the current status of the project with special emphasis on the development of its optics, and the plans for its construction and operation.

Alfvén waves are widely believed to play an important role in the transport of energy from the solar photosphere to the corona through the partially ionized chromosphere. In previous work, the properties of torsional Alfvén waves were theoretically studied using a multi-fluid model. Here, we compare those multi-fluid results with those obtained using the single-fluid magnetohydrodynamic approximation, as a way to assess the performance of the latter in the context of Alfvénic waves in the lower solar atmosphere. We consider a broadband photospheric driver that excites torsional Alfvén waves with frequencies ranging from 0.1 mHz to 300 mHz. These waves propagate upwards to the corona along a magnetic flux tube expanding with height. For both models, we compare the energy flux, chromospheric reflection, transmission and absorption coefficients, and the associated heating rates. In general, the results are almost identical in the two models, with the exception of two minor differences: (1) the net energy flux reaching the corona is approximately 5% larger in the single-fluid model, mainly owing to the higher reflectivity found in the multi-fluid model for wave frequencies exceeding 10 mHz; and (2) in a narrow region around 500 km above the photosphere, the single-fluid model underestimates the plasma heating rate due to ion-neutral damping by about a factor of two compared with the multi-fluid model. Both discrepancies arise from the approximate treatment of the ion-neutral drift in the single-fluid model and are expected to have a very limited impact on practical applications.

In the excursion-set formalism, the mass distribution of primordial black holes (PBHs) is derived from the first-passage time of a random walk describing the density contrast as the coarse-graining scale varies. We address two recent criticisms that have been raised about this approach. First, it was argued that the random walks are subject to colored (i.e. correlated over time) noise, making the first-passage-time problem cumbersome. We show that this arises from an incorrect separation of drift and noise when sampling on the Hubble-crossing surface: if Fourier modes are uncorrelated, the noise is strictly white. Moreover, sampling along the Hubble-crossing surface precludes using the density dispersion as a time variable, explaining the reported pathologies. Sampling instead on a synchronous surface removes both issues. This requires solving a first-passage-time problem with a moving barrier, for which we provide an efficient numerical framework. Second, it was suggested that cloud-in-cloud (i.e. that large black holes may engulf smaller ones) is irrelevant for PBHs and that the excursion set is therefore not needed. While valid for widely separated scales, this statement fails for broad power spectra with enhanced continua of modes. We further show that Press-Schechter estimates neglecting boundary evolution can break down even without cloud-in-cloud effects. Our results establish the robustness and necessity of the excursion-set formalism in realistic PBH formation scenarios.

Surface runoff shapes planetary landscapes, but global hydrological models often lack the resolution and flexibility to simulate dynamic surface water bodies beyond Earth. Recent studies of Mars have revealed abundant geological and mineralogical evidence for past surface water, including valley networks, crater lakes, deltas and possible ocean margins dating from late Noachian to early Hesperian times. These features suggest that early Mars experienced periods allowing liquid water stability, runoff and sediment transport. To investigate where surface water could accumulate and how it may have been redistributed, we developed a global high-resolution (km-scale) surface hydrological model. The model uses a pre-computed hydrological database that maps topographic depressions, their spillover points, hierarchical connections between basins, and lake volume-area-elevation relationships. This database approach greatly accelerates simulations by avoiding repeated geomorphic processing. The model dynamically forms, grows, merges and dries lakes and putative seas without prescribing fixed coastlines, by transferring water volumes between depressions according to their storage capacities and overflow rules. We explore model behavior over the present-day Mars' topography measured by MOLA (Mars Orbiter Laser Altimeter) topography for a range of evaporation rates (from 0.1 m/yr to 10 m/yr) and total water inventories expressed as Global Equivalent Layer (from 1 mGEL to 1000 mGEL). 48 Simulations are iterated to reach the steady state. The model outputs the extent and depth of surface water bodies and identifies main drainage pathways using overflow fluxes as runoff indicators. Results show a transition toward a contiguous northern ocean between low (1-10 m) GEL values and increasing concentration of water in northern lowlands and major impact basins at higher GEL.

We present the first spatially resolved at $\sim20$ pc scale application of AGN-specific metallicity diagnostics for the nearby Compton-thick Seyfert 2 galaxy Mrk 573 ($z = 0.017$). We use Hubble Space Telescope narrow-band imaging, MUSE integral-field spectroscopy and apply AGN strong-line metallicity diagnostics based on [O III], [S II], H$\beta$, H$\alpha$, and [N II] emission lines. We construct maps of $12 + \log$(O/H) for two different metallicity calibrations and two different N/O-O/H scaling relations out to $\sim1$ kpc and down to $\sim20$ pc scales. Our analysis reveals metallicity enhancement in AGN-dominated regions, with oxygen abundances reaching up to few times Solar. The metallicity shows a patchy spatial distribution, varying on $\sim100$ pc scales, appears to trace the high Seyfert/LINER index (SLI) value regions and the VLA 6 cm jet/radio lobe emission. These spatial correspondences and the lack of evidence for star formation in the bicone region suggest that the enrichment originates from metals transported from the nuclear AGN regions by winds, outflows, or jets.

Since 2019, three scintillator detectors of the EEE collaboration have been continuously measuring cosmic muon rates at 78.9°N at the Ny-Ålesund Research Station (Svalbard). The resulting six-year time series reveals a pronounced annual modulation, driven primarily by seasonal atmospheric variations. Utilizing routine radiosonde profiles collected above the same site, we applied several established techniques --along with a tailored analysis approach-- to investigate the relationship between muon rate and atmospheric temperature. The temperature-corrected muon-rates are analysed using the Lomb-Scargle periodogram technique in order to investigate the presence of remaining periodic structures. Finally, the temperature corrections coefficients of our analysis are compared with measurements in other stations located at lower latitudes.

Wave-optics effects in strongly lensed gravitational waves (GWs) provide a new interferometric probe of dark matter substructure. We compute the full diffraction integral for GWs propagating through statistically generated cold dark matter subhalo populations and quantify the resulting frequency-dependent amplification in the Laser Interferometer Space Antenna (LISA) band. We show that realistic galaxy-scale lenses generically produce percent-level amplitude and phase distortions in strongly magnified images, primarily induced by subhalos in the mass range $10^4$-$10^7\,M_{\odot}$. These signatures arise naturally within the standard cold dark matter paradigm and should be detectable in high signal-to-noise LISA events. Strongly lensed GWs thus offer a direct and complementary window on dark matter structure at subgalactic mass scales inaccessible to electromagnetic measurements.

We present a novel method for the tomographic reconstruction of the bias-weighted mean electron pressure $\langle bP_e \rangle$ and star formation rate density $\langle b\rho_{\mathrm{SFR}} \rangle$, by simultaneously modelling the contribution from the thermal Sunyaev-Zel'dovich (tSZ) effect and the Cosmic Infrared Background (CIB) to the cross-correlation between photometric galaxy samples and multi-frequency Cosmic Microwave Background (CMB) maps. The resulting measurements are independent of the galaxy clustering properties and robust against cross-contamination between tSZ and CIB. Applying this method to publicly available data, we reconstruct the cosmic evolution of $\langle bP_e \rangle$ and $\langle b\rho_{\mathrm{SFR}} \rangle$ out to $z\sim1$, making our measurements publicly available. Our measurements of both quantities are broadly compatible with predictions from the fiducial FLAMINGO hydrodynamical simulation, although we observe a lower gas pressure at low redshifts, in agreement with other measurements.

Strong gravitational lensing provides an independent and powerful probe of cosmic expansion by directly linking observables to cosmological distances. Upcoming surveys such as LSST will discover large number of galaxy-galaxy strong lensing systems, offering a new route to precise cosmological constraints. In this paper, we propose a Fisher-like sensitivity factor to map how the cosmological information of strong-lensing distances changes across the lens-source redshift plane. Applying such factor to the distance ratio $D_{ls}/D_s$, the time-delay distance $D_{\Delta t}$, and the double-source-plane ratio, we determine the ``sensitivity valleys'' where an observable becomes insensitive to a given parameter. The realistically simulated LSST lens population, which largely lies outside the distance-ratio valleys, covers the most sensitive region for $(w_0,w_a)$ parameter space. We then develop a new hierarchical framework, which could calibrate the redshift evolution of lens mass-density slopes and constrain cosmological parameters simultaneously. Focusing on the LSST mock data, we demonstrate that ignoring mass-profile evolution can bias $\Omega_m$ by up to $\sim 10\sigma$, while modeling the lens evolution could perfectly recovers the fiducial cosmology and yield stringent cosmological constraints (e.g., $\Delta\Omega_m \simeq 0.01$ and $\Delta w \simeq 0.1$ for $\sim 10^4$ lenses).

We present an alternative measurement of the Baryon Acoustic Oscillation (BAO) using Legendre multipole representation of the Ly$\alpha$ forest correlation functions from the second data release (DR2) of the Dark Energy Spectroscopic Instrument survey. Compressing the auto- and cross-correlation functions into Legendre multipoles yields a positive-definite covariance matrix without any smoothing -- unlike the baseline DR2 analysis -- thanks to a significantly reduced data vector size. We introduce the statistical corrections required to debias the finite-sample covariance matrix estimate and demonstrate that monopole and quadrupole terms for both auto- and cross-correlations can be used even when the correlation functions are distorted by continuum errors and contaminated by metals. This formalism has slightly diminished the constraining power of the BAO scale, while considerably weakening constraints on nuisance parameters. We measure the isotropic BAO scale with $0.96\%$ precision at $z_\mathrm{eff}=2.35$, the Hubble parameter $H(z_\mathrm{eff})=(238.7\pm3.4)~(147.09~\mathrm{Mpc}/r_d) ~\mathrm{km~s}^{-1}~\text{Mpc}^{-1}$, and the transverse comoving distance $D_M(z_\mathrm{eff})=(5.79 \pm 0.10)~(r_d/147.09~\mathrm{Mpc})$~Gpc for a given value of the sound horizon ($r_d$). Our BAO results are entirely consistent with the baseline DR2 analysis.

We present near-infrared (NIR) observations of the 2021 eruption of the recurrent nova RS Ophiuchi. The dataset provides both pre- and post-eruption perspectives on the eruption, as well as NIR spectra at high cadence. The spectrum obtained in 2020 June (14.3 years after the 2006 eruption, and 428.1 days before the 2021 eruption), is that of the red giant secondary, on which are superimposed several emission lines which most likely arise in the red giant wind. Spectra obtained during the eruption consist of emission (including coronal) lines, superimposed on a bremsstrahlung continuum at 8900K. The temperature of the coronal gas is estimated to be $10^{6.0}$K on day 11.7, and $10^{5.9}$K on day 31.7. The high cadence observations, obtained on day 31.7 of the eruption, provide no conclusive evidence for rapid ($<\sim1$~minute) variations in the HeI 1.0833$\,\mu$m line. Data obtained about one year after the eruption show that there may have been changes in the spectral type of the secondary after the 2021 eruption.

Small-scale vortical motions in the upper solar atmosphere are abundant and occupy about 2.8% of the photosphere at any given time. Although considerable work has focused on the detection and analysis of individual solar vortices, the interconnected and multi-scale behaviour of these coherent structures remains largely unexplored. We present a methodology for studying this behaviour through vortex interactions, to improve our understanding of how small- and large-scale photospheric flows contribute to energy transfer into the upper solar atmosphere and to the driving of solar activity. We represent vortices as a network of interacting structures. We apply a community detection algorithm to derive an optimal reduced network composed of highly interconnected vortex groups. From the interaction patterns and group structure, we define three roles within each community: peripheral, connector and hub. We then track both vortex communities and their member vortices from the photosphere into the chromosphere and across their lifetimes. On average, vortices assigned to these roles persist to greater heights in the chromosphere and have longer lifetimes than unclassified vortices. This shows that community detection can identify vortices with greater dynamical influence on the upper solar atmosphere. We also find that 32% to 58.6% of vortex communities exhibit global periodic behaviour following a helical path. This collective vortical motion may indicate an enhanced mechanism for wave excitation. Solar vortical community detection, therefore, offers a new framework for studying solar vortices and a new perspective on the importance of collective vortex dynamics.

Methanol, the simplest complex organic molecule found in space, is considered a key compound necessary for the formation of chemical species of prebiotic interest. Methanol detections in protoplanetary disks remain scarce, even though it is frequently detected in the material surrounding other Young Stellar Objects. We investigate the presence of methanol in the protoplanetary disk around the HL Tau protostar, motivated by the detection of spatially resolved warm water emission. Given the similar volatility of methanol and water, thermally desorbed gas-phase methanol is expected to emit from the same region of the HL Tau disk where water vapour has been observed. Accordingly, we selected and imaged the most promising ALMA archival observations to search for rotational methanol lines. We found no methanol emission in the analysed archival datasets. Assuming optically thin emission and LTE, we derive stringent upper limits on the methanol column density for different excitation temperatures: < 7.2 x 10^(14) cm^(-2) at 100 K and < 1.8 x 10^(15) cm^(-2) at 200 K, assuming a circular emitting region with a radius of 17 au (~ 0.12''). Furthermore, we obtain a stringent upper limit on the methanol-to-water column density ratio (< 0.55 x 10^(-3) at 100 K and < 1.4x 10^(-3) at 200 K), which is, on average, an order of magnitude lower than the values measured for other Young Stellar Objects and Solar System comets. We argue that the most likely explanation for the methanol non-detection in HL Tau is the presence of optically thick dust in the central region of the disk, which obscures part of the methanol emission. The upper limit on the methanol-to-water ratio in the HL Tau disk is at least an order of magnitude smaller than most clouds, YSOs and comets, possibly due to radiative transfer and/or excitation effects, or due to a different chemical evolution compared to the other sources.

We present an automated framework to search for optical counterparts of LIGO-Virgo-KAGRA (LVK) gravitational wave (GW) superevents using public Zwicky Transient Facility (ZTF) alerts processed through the ALeRCE broker. The goal is to filter and identify optical transients potentially associated with binary black hole (BBH) mergers during the LVK O4a and O4b observing runs. Using the Automatic Learning for the Rapid Classification of Events (ALeRCE) infrastructure, we spatially query ZTF alerts within GW localization regions and apply machine learning classifiers, host-galaxy crossmatching, and temporal cuts within 200 days post-merger to isolate plausible candidates. Our search yielded one candidate in O4a and four in O4b, several consistent with the supernova or tidal disruption event regime. This work demonstrates that public alert brokers can establish a robust baseline for systematic searches for electromagnetic counterparts to GW superevents in current and future observing runs. Our algorithm provides a systematic approach to search for BBH counterparts for all significant LVK GW superevents using survey telescope alerts. The search, together with the accompanying analysis, demonstrates the significance of the counterpart candidates, with one candidate ultimately identified as a transient event consistent with a Bowen fluorescence flare in a now discarded active galactic nucleus (AGN).

The large population of broad-line Active Galactic Nuclei (AGN) observed with the James Webb Space Telescope (JWST) at $z \gtrsim 4$ opens a new window onto the black hole-galaxy connection in the first Gyr of cosmic history. We use the JADES survey-level dataset and develop a forward-modeling Bayesian framework that explicitly accounts for broad H$\alpha$ detectability, ensuring that selection effects are incorporated into the likelihood function. With this approach, we constrain the black hole-stellar mass ($M_{\mathrm{BH}}$-$M_\star$) relation to be $\log M_{\rm BH} = -4.06^{+0.50}_{-0.51} + 1.17^{+0.06}_{-0.06}\,\log M_\star$, with an intrinsic orthogonal scatter of $\sigma_{\rm int} = 0.63^{+0.14}_{-0.11}$ dex. The slope and normalization are consistent with local determinations, indicating that the average scaling was already established by $z \sim 4$-6. This suggests that the primary evolution of the relation occurs in its dispersion rather than in its mean normalization. In contrast, the substantially larger intrinsic scatter relative to the nearby Universe reveals a wider diversity of black hole-galaxy growth histories, likely driven by bursty accretion, delayed feedback, and differences in merger or seeding histories. Future JWST samples will be crucial to test whether this increased scatter is a persistent feature of the high-redshift Universe.

The study of Stark broadening of neutral helium lines, despite significant advances over recent decades, has not led to updated large grids of helium line profiles relevant to the spectroscopic study of helium-rich stars. While the semi-analytical approach based on the standard Stark broadening theory is efficient for generating such grids, it presents challenges in incorporating additional physical effects into the model. Motivated by recent studies that highlight potential issues with line profiles in the context of white dwarf stars, this paper leverages advances in computer simulations to create a new grid of line profiles for 13 neutral helium lines in the optical range. These profiles cover densities ranging from 10^14 to 6 x 10^17 cm^-3 and temperatures from 10,000 K to 40,000 K, with the exception of the narrower He I 4713 line, for which the profile grid begins at 10^15.5 cm^-3. The primary goal of this research is to present the new grid and compare it with both the semi-analytical approach and other simulation results. By doing so, corrections to the previous grid will be explored, providing a foundation for future studies that utilize this updated grid. We also examine the impact of these new profiles on the determination of physical parameters for a range of astrophysical objects, including DB white dwarfs and other helium-rich stars.

We present a novel method for constructing high-accuracy, time-varying disk-integrated stellar absorption line profiles that isolate the effects of granulation alone. This framework provides an effectively unlimited supply of physically consistent training data, offering a unique opportunity to study granulation-driven velocity variability with no contamination from other stellar processes or instrumental systematics. Our interpolation scheme enables accurate profile generation at arbitrary limb angles and successfully reproduces observed disk integrated solar bisector shapes from IAG spectra. Using four Fe I lines (525.0, 615.2, 617.3, and 627.1 nm), we produce 1000 model star disk-integrated realisations per line and find an isolated granulation-induced RV scatter of 0.16-0.21 m s^-1. Using our synthetic profiles and assuming infinite signal-to-noise, we find strong correlations between various line-shape metrics and convective blueshift, demonstrating that line-shape diagnostics can, in principle, trace granulation effects. Equivalent width proves the strongest diagnostic, achieving up to 60% scatter reduction. However, the strength of all simple line shape diagnostics rapidly diminishes once photon noise is injected. Even when artificially boosting the signal to represent a spectrum containing ~1000 spectral lines, the achievable improvement with these metrics remains below 10% at typical signal-to-noise ratios. Our results highlight the need for more robust, noise-resilient diagnostics and position our synthetic dataset as a valuable testbed for developing and benchmarking such methods.

The discovery rate of optical transients will explode to 10 million public alerts per night once the Vera C. Rubin Observatory's Legacy Survey of Space and Time comes online, overwhelming the traditional physics-based inference pipelines. A continuous-time forecasting AI model is of interest because it can deliver millisecond-scale inference for thousands of objects per day, whereas legacy MCMC codes need hours per object. In this paper, we propose SELDON, a new continuous-time variational autoencoder for panels of sparse and irregularly time-sampled (gappy) astrophysical light curves that are nonstationary, heteroscedastic, and inherently dependent. SELDON combines a masked GRU-ODE encoder with a latent neural ODE propagator and an interpretable Gaussian-basis decoder. The encoder learns to summarize panels of imbalanced and correlated data even when only a handful of points are observed. The neural ODE then integrates this hidden state forward in continuous time, extrapolating to future unseen epochs. This extrapolated time series is further encoded by deep sets to a latent distribution that is decoded to a weighted sum of Gaussian basis functions, the parameters of which are physically meaningful. Such parameters (e.g., rise time, decay rate, peak flux) directly drive downstream prioritization of spectroscopic follow-up for astrophysical surveys. Beyond astronomy, the architecture of SELDON offers a generic recipe for interpretable and continuous-time sequence modeling in any time domain where data are multivariate, sparse, heteroscedastic, and irregularly spaced.

A complete picture of the gas thermodynamics around galaxies is imprinted on the cosmic microwave background (CMB). Indeed, the thermal, kinematic, and relativistic Sunyaev-Zel'dovich effects (tSZ, kSZ, rSZ) measure the gas density, temperature, pressure of baryonic feedback and bulk velocity around galaxies, along with the gravitational potential it sits in. This full thermodynamic picture promises to constrain galaxy formation models and gas related uncertainties in the impact on galaxy lensing. Recent kSZ measurements around galaxies suggest that the gas may be more extended than anticipated, pointing to powerful feedback processes and large baryonic corrections to lensing. How robust are these conclusions about the galaxy-halo connection, including satellite fraction and high-mass outliers, or to 2D projection effects and large-scale velocity modes? In this paper, we give estimates for these effects using a simulated sample of DESI-like luminous red galaxies within the IllustrisTNG hydrodynamical simulation and the Abacus N-body simulation. We show that analyzing projected 2D profiles can lead to biases when computing quantities like the gas fraction. We also find that in the absence of spatial filtering, the 2-halo term is non-negligible for kSZ even at the smaller radii where the 1-halo term dominates. We show that a 1% uncertainty in the satellite fraction of galaxies can propagate into uncertainties of $\pm$1%, $\pm$3% and $\pm$5% in the 1-halo terms of the kSZ, tSZ, and rSZ signals, respectively. We show that masking the 2% most massive objects in the sample reduces the profile amplitudes by up to 10%, 40%, and 75% for the kSZ, tSZ, and rSZ signals, respectively. Finally, we show that naive simulations of the kSZ effect can be biased by an artificial Doppler term, which is automatically removed when high-pass or compensated aperture filtering is applied.

Among the most compelling cold dark matter candidates, the axion has recently been subject to a wide range of astrophysical studies aiming to constraints its properties. We present updated bounds on the isocurvature fraction, $f_{\rm{iso}}$, which parameterizes the contribution of isocurvature perturbations induced by post-inflationary produced axion-like particles (ALPs) to the ordinary power spectrum. We use new simulations based on the Sherwood-Relics suite to fit high-resolution Lyman-$\alpha$ forest flux power spectrum data. With the published noise model of the Lyman-$\alpha$ forest data, we find a tentative detection of $f_{\rm{iso}}$ = ${0.0064^{+0.0012}_{-0.0014}}$ (68% C.L), after accounting for the degenerate effect of IGM thermal evolution. With a more conservative modelling of the residual noise in the data, the upper bound is weakened to $f_{\rm{iso}}< 0.0084$ (95% C.L), which translates into an ALP temperature-independent mass $m_a > 1.73 \times 10^{-18}$eV. Our constraints are stronger than bounds derived from large-scale structure probes at higher and lower redshifts and are competitive with those derived from UV luminosity function data. Interestingly, the best current Lyman-$\alpha$ forest data prefers a non-zero contribution from isocurvature modes.

Carbonaceous Chondrites have special significance in the stellar evolution and in particular in the evolution of life on earth. The carbonaceous meteorite that fell in Mukundpura village, Jaipur, Rajasthan on 6th June 2017 is one such rare CM2 (Carbonaceous Chondrite) carbonaceous meteorite. We carried out high resolution scanning and transmission electron microscopic (TEM) studies on typical thin sections, showing abundant grains of iridium (Ir), pentlandite (NiS), and more interestingly crystalline carbon (C). These crystallite carbon grains resemble nanodiamond like signature in the freshest Mukundpura meteorite. The high-resolution Raman spectroscopic measurements are carried out on the crystalline carbon grains, showing well resolved three distinct peaks with a vibrational mode at 1315 cm-1, with the onset of a weak vibrational mode at 1150 cm^-1, substantiating the observation of nanocrystalline diamond in Mukundpura meteorite. The broad peak centered at 1360 cm^-1 and 1575 cm^-1 (as an average), suggest the presence of graphitic carbon as well together with apparent presence of nanocrystalline diamond. The average size of nanocrystalline diamond is ~ 3-5 nm. High iridium content in this meteorite supports the meteoric impact related iridium anomaly in geological stratigraphic boundaries (this http URL-Tertiary boundary) that has caused mass extinction of flora and fauna.

Because of the special angular distribution of excited electrons by the photoelectric effect, the Gas Pixel Detector (GPD) is effective in measuring keV X-ray polarization of astrophysical events (e.g. gamma-ray bursts), by capturing ionization tracks of excited electrons as polarized images. Traditionally, the emission angles of photoelectrons are extracted from polarized images first, and statistics are then performed on these angles to infer the polarization direction and intensity. However, observation with the wide field of view requires the incident angle of X-rays not directly attainable through the traditional analysis process. In this paper, we propose using the generalized sliced Wasserstein (GSW) distance, projected by neural networks with random weights, as a completely data-driven approach to analyze X-ray polarization based on two-dimensional polarized images. We find the structures of the randomized neural networks matter when focusing on different aspects of the polarized images, and take advantage of the discrimination abilities by different neural network structures. The proposed method, named the structured GSW distance, successfully distinguishes polarized images with different configurations of incident angles and polarization directions. Furthermore, we build a simplified statistical model based on the von Mises distribution and the circular Wasserstein distance and compare the model against the proposed method, showing their high consistency. The computational method reported in this paper may benefit GPD-based polarimetry in astroparticle experiments and also pattern analysis on raw data from pixel detectors.

We will address the existence of a new symmetry for an imperfect fluid by introducing local four-velocity gauge-like transformations for the case when there is vorticity. A similar tetrad formulation as to the Einstein-Maxwell spacetimes formalism presented in previous manuscripts will be developed in this manuscript for the imperfect fluids. The four-velocity curl and the metric tensor will be invariant under these kind of four-velocity gauge-like local transformations. While the Einstein-Maxwell stress-energy tensor is locally gauge invariant under electromagnetic gauge transformations, the perfect fluid stress-energy tensor will not be invariant under four-velocity gauge-like local transformations. We will dedicate our analysis to the imperfect fluid stress-energy tensor that will be invariant under local four-velocity gauge-like transformations when additional transformations are introduced for several variables included in the stress-energy tensor itself. We will also pay special attention to the construction of a vorticity stress-energy tensor invariant under local four-velocity gauge-like transformations. An application on neutron stars will be developed in order to show the simplifications brought about by these new tetrads.

The non-linearly coupled Einstein-matter field equations predict the existence of shadows with well-defined boundaries around black holes. We prove that, in spherically symmetric hairy black-hole spacetimes whose matter fields satisfy the weak energy condition, the radii of these shadows are bounded from below by the dimensionless relation $r_{\text{sh}}/r_{\text{H}}\geq 3\sqrt{3}/2$, where $r_{\text{H}}$ is the horizon radius of the central hairy black hole. The characteristic shadow of the (bald) Schwarzschild black-hole spacetime saturates the analytically derived lower bound.

Typical gravitational production of relics from amplification of inflationary perturbations assumes Bunch-Davies initial conditions, i.e. a vacuum with initially no particles. In this paper we investigate the impact of non Bunch-Davies initial conditions to the final abundance of relics, with particular attention to the parameter space where the total dark matter abundance is reproduced. We present a general framework for any initial condition, through which we show their non-trivial effect on both spectrum and late-time abundance. We argue that for particles whose source of conformal symmetry breaking comes only from a mass term (spin-1/2 fermions and conformally coupled scalars), the choice of initial conditions has little impact on the mass range relevant to dark matter. For other particles, e.g. the longitudinal mode of spin-1, we see a large deviation from the standard computation. We exemplify and quantify our results with an initial thermal state and a two-stage inflation scenario, highlighting that the total dark matter can be obtained for a wide range of masses.

Axions and other very weakly interacting slim particles (WISPs), with masses below 1 GeV, arise naturally in many extensions of the Standard Model of particle physics. In particular, they could offer a new framework to explain the nature of dark matter and may help address a range of puzzling observations in astrophysics and particle physics. This review provides an overview of ongoing WISP searches and outlines the prospects for the next decade, spanning their theoretical motivation, indirect signatures in astrophysical observations, and dedicated laboratory experiments. It is based on the work carried on by the EU-funded COST Action ``Cosmic WISPers in the Dark Universe: Theory, astrophysics, and experiments'' (CA21106, this https URL). This network plays a key role in coordinating and supporting WISP searches across Europe, while also contributing to the development of a roadmap aimed at securing European leadership in this research area. It is emphasized that Europe is currently pursuing a rich, diverse, and cost-effective experimental program, with the potential to deliver one or more transformative discoveries.

Dark matter direct detection experiments involving electron recoils are beginning to test highly-predictive, thermal-relic milestones for sub-GeV dark matter models. Due to the Lee-Weinberg bound, thermal dark matter candidates in this mass range necessarily require comparably-light mediator particles to achieve a suitably large annihilation cross section. Here we present new thermal-relic milestones for sub-GeV dark matter candidates that couple to vector mediators. In these models, the mediators are massive gauge bosons of anomaly-free abelian extensions to the Standard Model, including the dark photon, gauged $L_i - L_j, B-L$, and $B-3L_i$ models, where $B$ is the baryon number, $L$ is the lepton number, and $i,j$ index the lepton families. Since the same interactions that govern cosmological production also govern electron scattering, the targets we present are firmly predictive and allow for these models to be robustly discovered or falsified. Furthermore, since the mediators we study exhaust the minimal anomaly-free U(1) extensions to the Standard Model, our results offer a complete list of predictive milestones for sub-GeV dark matter coupled to vector mediators.

Low-thrust electric propulsion missions are often designed under simplifying assumptions such as constant thrust or fixed specific impulse, neglecting the strong coupling between trajectory dynamics, spacecraft power availability, and propulsion performance. In deep-space environments with reduced solar irradiance, these assumptions can lead to suboptimal or infeasible designs, underscoring the need to simultaneously optimize the trajectory and power subsystem. This paper presents a multidisciplinary design optimization (MDO) framework for the simultaneous design of low-thrust trajectories and spacecraft power systems, with explicit coupling to electric propulsion performance. The framework incorporates a high-fidelity variable-specific impulse model of the SPT-140 Hall thruster, in which thrust and efficiency are directly constrained by time-varying solar power availability and solar array degradation, rather than treated as fixed parameters. The coupled problem is posed as a time-optimal control problem and addressed using a framework built on top of OpenMDAO and Dymos toolchains, where Dymos employs a collocation-based direct-transcription approach for trajectory optimization. OpenMDAO provides accurate analytic partial derivatives, enabling efficient gradient-based optimization. A Fast Fourier Series shape-based method is used to generate dynamically feasible initial guess trajectories, and the resulting nonlinear programming problem is solved using IPOPT. The proposed framework is demonstrated through a low-thrust orbit insertion scenario around asteroid 16-Psyche, a regime in which reduced solar irradiance makes power-aware trajectory design particularly critical. Simulation results demonstrate the framework's ability to capture key power-propulsion-trajectory trade-offs, highlighting the importance of integrated power optimization for realistic electric propulsion mission design.

Recently, an exact rotating black hole solution in a parity-violating theory of gravity was obtained via a conformal transformation of the Kerr solution in general relativity, with parity-violating effects encoded in the conformal factor. We study the quasinormal modes (QNMs) of a test scalar field minimally coupled to gravity on this conformal Kerr background, treating the parity-violating effects perturbatively while allowing for arbitrary black hole spin, from the non-rotating case to the near-extremal regime. For low spin, we derive a perturbative formula for the QNM frequencies that includes the leading-order parity-violating correction. For high spin, particularly in the near-extremal regime, we find sizable deviations from the Kerr QNM frequencies. Our results point to a new avenue for probing parity-violating physics in the strong-gravity regime through black hole QNMs.

We investigate the principles of quantum field theory using a stiff de Sitter space. We demonstrate that a non-unitary Lagrangian on a Euclidean AdS geometry can produce the perturbative expansion of late-time correlation functions to all orders. This discovery greatly simplifies perturbative computations while also allowing us to prove fundamental features of these correlators, which are part of a Euclidean CFT. This allows us to construct an OPE expansion, limit the operator spectrum, and deduce the analytic structure of the spectral density that captures the conformal partial wave expansion of a late-time four-point function. In general, the standard CFT concept of unitarity does not apply to dimensions and OPE coefficients. Rather, the positivity of the spectral density represents the unitarity of the de Sitter theory. This assertion is non-perturbative and does not depend on the use of Euclidean AdS Lagrangians. In a scalar theory, we compute tree-level and entire one-loop-resummed exchange diagrams to demonstrate and verify these characteristics. In the spectrum density, an exchanged particle shows up as a resonant characteristic that may be helpful in experimental searches.

Modeling inviscid compressible flows with shocks and vortex dominated dynamics remains challenging for particle grid methods due to moving discontinuities, cell crossing noise, and quadrature degradation under strong deformation. Building on a FLIP/APIC framework with vorticity aware tensor artificial viscosity, we identify a long time RTI failure mode: particle depletion at spike heads degrades quadrature and particle grid coupling, producing nonphysical, void-like dents. Standard mitigations (CPDI lite and subcell-jittered seeding) reduce but do not eliminate this artifact. We therefore add two sampling-aware controls: (i) conservative split resampling that replenishes depleted cells while exactly conserving mass, momentum, and internal energy; and (ii) a soft-switch that attenuates only the APIC affine term when local support is insufficient. Tests on the Sod shock tube and single/multi mode RTI show that the method removes spike head voids in long-time RTI while preserving vortex roll up, and matches reference Euler growth metrics

A new phenomenological dark energy model, originally associated to the large-scale structure formation and considered as a solution to the fine-tuning and coincidence problems related to the cosmological constant, was analyzed within the framework of General Relativity in a Friedman-Robertson-Walker spacetime and its model parameters were estimated using cosmic chronometers and recent DESI data. It turns out that the proposed model can serve as an alternative evolving dark energy model with a novel equation of state function, apart from other popular propositions in the literature. Due to the form of this phenomenological energy density ansatz, which starts to rise with the nonlinear structure growth in the universe and falls with the domination of cosmic voids, we prefer to call it structure-induced dark energy. Observational constraints show that it is not only a suitable solution for the fundamental problems such as coincidence or fine-tuning problems, it gives flexibility, when considering the cosmic tensions and presents a new perspective on the evolving dark energy models.

In this letter, we explore the phenomenological impact of inflationary gravitational particle production in the physics of Dark Matter (DM). Large-scale DM fluctuations generated during inflation behave as gravitational particles upon their post-inflationary horizon reentry and alter the conventional Boltzmann dynamics of DM with a non-conserving source term, thereby producing significant phenomenological consequences. Within this framework, we analyze four distinct types of DM classified according to their production mechanisms. Dark matter may be completely non-interacting with the thermal bath, behaving as Inert Dark Matter. Alternatively, depending on the strength of its interactions with bath particles, DM may exhibit WIMPy, UFOy, or FIMPy behavior, sharing characteristics with their conventional counterparts. The late-time enhancement of the DM number density, driven by the successive horizon reentry of gravitationally produced low-momentum modes, enlarges the viable parameter space for both thermal and non-thermal DM scenarios. Remarkably, this expanded parameter space remains consistent with current constraints from $\Delta N_{\rm eff}$ and Lyman-$\alpha$ bound.

We study slow-roll inflation driven by a scalar field non-minimally coupled to gravity within the effective framework of Loop Quantum Cosmology (LQC), including inverse-volume corrections. We consider two physically motivated classes of potentials, a Higgs-like quartic potential $V\propto\phi^{4}$ and string-inspired fractional monomial potentials $V\propto\phi^{p}$ with $p<1$. Working at first order in the slow-roll expansion, we derive analytic expressions for the inflationary observables, namely the scalar spectral index $n_s$, the tensor-to-scalar ratio $r$, and the running $\alpha_s\equiv dn_s/d\ln k$, and then solve the corrected background dynamics numerically to obtain quantitative predictions. Confronting these results with current observational constraints from Planck 2018 and ACT DR6, we find that the model can lie within the allowed region of the $(n_s,r,\alpha_s)$ parameter space, including a mild preference for slightly larger $n_s$, as suggested by recent ground-based measurements. We also compute the probability of achieving sufficient slow-roll inflation in this setting. Although effective LQC replaces the initial singularity with a nonsingular quantum bounce, the likelihood of a sufficiently long inflationary phase depends on the pre-inflationary dynamics and on the inflaton potential. Using the canonical Liouville measure on the effective phase space, we determine the fraction of post-bounce trajectories that yield sufficient inflation and find that the non-minimal coupling parameter $\xi$ substantially enlarges the phase-space volume of favorable initial conditions relative to the minimally coupled case, exhibiting an attractor-like enhancement that saturates at large $\xi$.

Two-dimensional particle-in-cell simulations with a realistic mass ratio reveal the generation mechanisms of the out-of-plane magnetic field in magnetic islands/flux ropes during magnetic reconnection. In the absence of an initial guide field, reconnection produces a large electron temperature anisotropy (around 4.5) inside magnetic islands that drives the Weibel instability. Strong out-of-plane magnetic fields (Bz/B0 around 0.4, greatly exceeding the Hall field) with a regular bipolar structure grow inside islands. A space-time analysis reveals a one-to-one correspondence between the temperature anisotropy and the development of the Weibel magnetic field. The instability relaxes the anisotropy, but island merging leads to anisotropy reemergence and re-excitation. In the presence of a strong ambient guide field (Bg/B0 = 0.5), the electron outflow from the X-point deflects along the separatrices and forms a circular current loop wrapping the flux ropes. This flux-rope separatrix current generates an out-ofplane magnetic field that reinforces the ambient guide field, reaching Bz/B0 around 1.4. The current can, in some cases, drive the electron Kelvin-Helmholtz instability, which produces electron vortices and strengthens the magnetic field. Mergers significantly broaden the islands and further strengthen the field. These self-generated out-of-plane magnetic fields scatter electrons and reduce their temperature anisotropy, which can potentially affect electron heating via Fermi reflection. The simulation results are supported by spacecraft observations suggesting that ambient guide fields can be enhanced within flux ropes in Earth's magnetotail.

Hermiticity is usually treated as a foundational axiom of quantum mechanics, guaranteeing real spectra and unitary time evolution. In this work we argue that Hermiticity is more naturally understood as a symmetry law arising from the global conservation of an inner product current. We show that in spacetimes admitting complete Cauchy surfaces without boundary flux this conservation reduces to the familiar Hermiticity condition of the canonical inner product. However, in the presence of causal horizons, most strikingly in black hole geometries, this conservation law becomes obstructed for restricted observers. Tracing over inaccessible degrees of freedom then inevitably yields completely positive trace preserving dynamics with an effective non-Hermitian generator. Using quantum thermodynamics and the monotonicity of relative entropy, we demonstrate that the generalized second law may be reinterpreted as an entropy balance that compensates precisely for the flux of inner product charge through the horizon. The structure of Einstein equations, through the Bianchi identity and the Raychaudhuri focusing equation, provides the geometric mechanism underlying this balance. We also show that black hole ringdown can serve as a realistic observational probe of this idea and may provide quantitative upper bounds on the strength of horizon-induced inner product flux. In this way gravity, entropy production, and effective non-Hermiticity are unified under a single structural principle, with Hermiticity emerging as the special case of globally conserved inner product symmetry.

We measure the cross-correlation between cosmic shear from the third-year release of the Dark Energy Survey, thermal Sunyaev-Zel'dovich (tSZ) maps from Planck, and X-ray maps from ROSAT. We investigate the possibility of developing a physical model able to jointly describe both measurements, simultaneously constraining the spatial distribution and thermodynamic properties of hot gas. We find that a relatively simple model is able to describe both sets of measurements and to make reasonably accurate predictions for other observables (the tSZ auto-correlation, its cross-correlation with X-rays, and tomographic measurements of the bias-weighted mean gas pressure). We show, however, that contamination from X-ray AGN, as well as the impact of non-thermal pressure support, must be incorporated in order to fully resolve tensions in parameter space between different data combinations. We obtain simultaneous constraints on the mass scale at which half of the gas content has been expelled from the halo, $\mathrm{log}_{10}(M_c)=14.83^{+0.16}_{-0.23}$, on the polytropic index of the gas, $\Gamma=1.144^{+0.016}_{-0.013}$, and on the ratio of the central gas temperature to the virial temperature $\alpha_T=1.30^{+0.15}_{-0.28}$.

Globular clusters (GCs) are some of the oldest bound structures in the Universe, holding clues to the earliest epochs of star formation and galaxy assembly. However, accurate age measurements of ancient clusters are challenging due to the age-metallicity degeneracy. Here, we report the discovery of 36 compact stellar systems within the 'Relic', a massive, quiescent galaxy at $z=2.53$. The Relic resides in an overdensity behind the Abell~2744 cluster, with a prominent tidal tail extending towards two low-mass companions. Using deep data from the UNCOVER/MegaScience JWST Surveys, we find that clusters formed in age intervals ranging from 8~Myr up to $\sim2$~Gyr, suggesting a rich formation history starting at $z\sim10$. While the cluster-based star formation history is broadly consistent with the high past star formation rates derived from the diffuse host galaxy light, one potential discrepancy is a tentative $\sim2-3\times$ higher rate in the cluster population for the past Gyr. Taken together with the spatial distribution and low inferred metallicities of these young-to-intermediate age clusters, we may be seeing direct evidence for the accretion of star clusters in addition to their early in-situ formation. The cluster masses are high, $\sim10^6-10^7~M_{\odot}$, which may explain why we are able to detect them around this likely post-merger galaxy. Overall, the Relic clusters are consistent with being precursors of the most massive present-day GCs. This unique laboratory enables the first connection between long-lived, high-redshift clusters and local stellar populations, offering insights into the early stages of GC evolution and the broader processes of galaxy assembly.

Using open astronomical multifrequency databases, we constructed light curves and developed a comprehensive visualisation and sonification analysis for the blazars Mrk~501, Mrk~1501, Mrk~421, BL~Lacerta, AO~0235+164, 3C~66A, OJ~049, OJ~287, and PKS~J2134-0153. This study employed Musical Instrument Digital Interface (MIDI) and Parameter Mapping Sonification (PMSon) techniques to generate waveforms, spectrograms, and sonifications. These representations demonstrate that data visualisation and sonification are powerful tools for analysing astronomical objects like blazars, providing insights into their multifrequency variability. This work highlights how sonification and visualisation can aid in identifying potential patterns, power variations, regularities, and gaps in the data. This multimodal approach also underscores the importance of inclusivity in scientific communication, offering accessible methods for exploring the complex behaviour of blazars.

(119951) 2002 KX14 is a large classical TNO with limited previous observations and unresolved questions regarding its physical properties. Five stellar occultations by 2002 KX14 were observed from 2020 to 2023, involving multiple telescopes across different locations in Europe and the Americas. The five occultations resulted in 15 positive chords, accurately measuring the 2002 KX14's shape and size. The projected ellipse has semi-major and semi-minor axes of $241.0 \pm 7.2$ km and $157.1 \pm 5.2$ km, respectively, corresponding to an average area-equivalent diameter of $389.2 \pm 8.7$ km. The geometric albedo was estimated at $11.9 \pm 0.7\%$.

Surface release of radiatively active particles, with high infrared- (IR-)to-visible extinction ratios, has been proposed as a method of warming Mars. However, to warm Mars using aerosols, particles released locally must disperse globally. Here we provide an initial reference study in a plume tracking, dry Martian atmospheric model to address this question. The winds that transport aerosols respond to the aerosol's IR forcing, implying strong radiative-dynamical feedbacks (RDF). We investigate RDF from surface release of two particle compositions: carbon (graphene) and metal (Al). Self-lofting helps particles rise and spread locally and regionally, and the Hadley cell strengthens under warming, aiding latitudinal mixing. Within our model, Mars RDF enable engineered-aerosol warming. Warming is slightly greater for three-dimensional vs. 1D-models and also depends on spectral resolution of radiative transfer. We assess implications for Mars warming. Many open atmospheric science questions remain, including the role of agglomeration, dry-deposition rate uncertainty, and modeling water cycle feedbacks.

Galaxy model fitting is widely employed to estimate properties such as galaxy shape, size, and color. Understanding how the outputs of galaxy model fitting respond to weak-lensing shear distortions is crucial for accurate shear estimation and mitigating shear-related systematics in weak lensing image analyses. In this paper, we investigate how the fitted parameters - specifically flux, size, and shape - respond to weak-lensing shear distortions within the AnaCal framework. To achieve this, we introduce quintuple numbers, a novel algebraic system inspired by dual numbers from automatic differentiation. Quintuple numbers enable the propagation of shear response information throughout the entire model-fitting process by linking analytical pixel shear responses to those of the fitted parameters. We integrate quintuple numbers into the AnaCal framework to derive the shear responses of shapes estimated with model fitting and validate the pipeline using image simulations that include realistic blending. Our results demonstrate that the multiplicative bias remains below 0.003 for ground-based, oversampled images.

We investigate the influence of large-scale cosmic web environments on galaxy quenching using a volume-limited, stellar mass-matched galaxy sample from SDSS DR18. Galaxies are classified as residing in sheets, filaments, or clusters based on the eigenvalues of the tidal tensor derived from the smoothed density field. The quenched fraction increases with stellar mass and is highest in clusters, intermediate in filaments, and lowest in sheets, reflecting the increasing efficiency of environmental quenching with density. A flattening of the quenched fraction beyond $\log_{10}(M_\star/M_\odot) \sim 10.6$ across all environments signals a transition from environment-driven to mass-driven quenching. In contrast, the bulge fraction continues to rise beyond this threshold, indicating a decoupling between star formation suppression and morphological transformation. At the high-mass end ($\log_{10}(M_\star/M_\odot) \gtrsim 11.5$), both quenched and bulge fractions bifurcate, increasing in clusters but declining in sheets, suggesting a divergent evolutionary pathway where massive galaxies in sheets retain cold gas and disk-like morphologies, potentially sustaining or rejuvenating star formation. The AGN fraction also increases with stellar mass and is somewhat higher in sheets than in clusters, indicating enhanced AGN activity in low-density, gas-rich environments. The high-mass trends are independently corroborated by our analysis of specific star formation rate, $(u-r)$ colour, concentration index, and D4000 in the stellar mass-density plane, which show that massive galaxies in sheets remain bluer, younger, more star-forming, and structurally less evolved than their cluster counterparts. Our results highlight the cosmic web as an active driver of galaxy evolution.

LuSEE-Night is a pathfinder radio telescope on the lunar far side employing four 3-m monopole antennas arranged as two horizontal cross pseudo-dipoles on a rotational stage and sensitive to the radio sky in the 1-50 MHz frequency band. LuSEE-Night measures the corresponding 16 correlation products as a function of frequency. While each antenna combination measures radiation coming from a large area of the sky, their aggregate information as a function of phase in the lunar cycle and rotational stage position can be deconvolved into a low-resolution map of the sky. We study this deconvolution using linear map-making based on the Wiener filter algorithm. We illustrate how systematic effects can be effectively marginalised over as contributions to the noise covariance and demonstrate this technique on beam knowledge uncertainty and gain fluctuations. With reasonable assumptions about instrument performance, we show that LuSEE-Night should be able to map the sub-50 MHz sky at a ~5-degree resolution.

Lensed quasars are powerful probes of cosmology, the co-evolution of supermassive black holes and their host galaxies, and the distribution of dark matter. We cross-match 1,724 previously identified candidates from KiDS, HSC, and the DESI Legacy Imaging Surveys (DESI-LS) with DESI DR1, obtaining 937 DESI spectra for 677 unique systems. Combining DESI spectroscopy with observations from the Palomar 200-inch Double Spectrograph (P200/DBSP), we confirm two lensed quasars with source redshifts $z_s={1.93,3.23}$ and Einstein radii $\theta_{\rm E}={0.39,1.07}$ arcsec, respectively. We further identify 12 likely lensed quasars that are well reproduced by a simple SIE model, exhibit lens galaxies in the image modeling, and have at least one available spectrum; across these, $\theta_{\rm E}$ spans $045$-$2.34$ arcsec and $z_s$ spans $1.13$-$2.88$. In 9 of the 12 cases, the systems already satisfy our lensing criteria except that only one quasar image currently has a spectrum; obtaining a second spectrum for the other image would enable immediate confirmation. Moreover, we report eight new static strong lenses spanning galaxy- to group- scale lenses. These results provide valuable targets for follow-up studies and underscore the efficiency of wide-field spectroscopic surveys such as DESI in confirming gravitationally lensed quasars and galaxies.

Wide-separation lensed quasars (WSLQs) are a rare subclass of strongly lensed quasars produced by massive galaxy clusters. They provide valuable probes of dark-matter halos and quasar host galaxies. However, only about ten WSLQ systems are currently known, which limits further studies. To enlarge the sample from wide-area surveys, we developed a catalog-based pipeline and applied it to the CatNorth database, a catalog of quasar candidates constructed from Gaia DR3. CatNorth contains 1,545,514 quasar candidates with about 90% purity and a Gaia G-band limiting magnitude of roughly 21. The pipeline has three stages. First, we identify groups with separations between 10 and 72 arcsec using a HEALPix grid with 25.6 arcsec spacing and a friends-of-friends search. We then filter by intra-group color and spectral similarity, reducing the 1,545,514 sources to 14,244 groups while retaining all known, discoverable WSLQs. Finally, a visual check, guided by image geometry and the presence of likely foreground lenses, yields the candidate list with quality labels. We identify 333 new WSLQ candidates with separations from 10 to 56.8 arcsec. Using available SDSS DR16 and DESI DR1 spectroscopy, we uncover two new candidate systems; the remaining 331 candidates lack sufficient spectra and are labeled as 45 grade A, 98 grade B, and 188 grade C. We also compile 29 confirmed dual quasars as a by-product. When feasible, we plan follow-up spectroscopy and deeper imaging to confirm WSLQs among these candidates and enable the related science.

We study how the stellar velocity dispersion within the Scorpius-Centaurus OB association (Sco-Cen) has evolved over approximately 20 million years, from its formation to the present day, by investigating 32 stellar clusters in Sco-Cen. Using data from the Gaia mission along with supplementary stellar radial velocities, we identified a surprising sequence of abrupt jumps and intervening plateaus in the evolution of velocity dispersion correlating with times of star formation bursts. We find that the association is almost isotropically expanding and that star formation propagated from inside-out with a speed of about 5-6 km/s. We measure a present-day expansion rate of about 10-12 pc/Myr and observe that younger star clusters within the association exhibit higher velocities compared to older ones. This result, along with the stepwise increase in velocity dispersion over time, suggests a structured and sequential star formation process rather than a random one. This phased evolution suggests that stellar feedback is the primary driver of Sco-Cen's star formation history, expansion, and eventual dispersal. Our findings emphasise the value of precisely characterising stellar populations within OB associations, particularly through the creation of detailed, high-resolution age maps.

We perform 3 dimensional moving-mesh hydrodynamical simulations of bubble nebulae around ultraluminous X-ray sources, using state-of-the-art software AREPO. We use a Monte-Carlo method to inject outflows with uniform mass outflow rate and momentum, in a conical funnel with a specific half opening angle. Simulation results show that the morphology of the bubble is determined by the initial momentum of the outflows, while the mechanical power of the outflows only influences the size of the bubble without changing its shape. Low mechanical power also results in a short cooling timescale of the system, leading to an early collapse of the bubble shell. The half opening angle of the outflows and the viewing angle of the system determine the observed bubble eccentricity together. Compared with the observational morphology of the ULX bubble sources NGC 55 ULX-1 and NGC 1313 X-2, our simulation favors the fact that the high velocity outflows of the accretion disks in these two systems are confined in a narrow funnel region.

We investigate the properties of extreme emission line galaxies (EELGs) at $z=4-9$ and their role in reionization. Compact, low-mass galaxies with intense optical emission lines are linked to elevated specific star formation rates (sSFRs) and recent bursts of star formation. Feedback in these systems may enable the leakage of ionizing radiation into the intergalactic medium. Using JWST/NIRSpec spectroscopy from the CAPERS, CEERS, and RUBIES surveys, we compile 160 NIRCam-selected EELGs in the EGS field. These galaxies show extreme rest-frame equivalent widths (EWs), with a median EW([O III]+H${\beta}$)=1616Å and EW(H${\alpha}$)=763Å. They are low-mass (median log(M$_{\star}$/M$_{\odot}$)=8.26) with high sSFRs (median 43 Gyr$^{-1}$), above the $z\sim6$ main sequence. UV slopes are diverse, with a mean $\beta=-2.0$, and only 7% have extremely blue continua ($\beta<-2.6$). Emission-line diagnostics suggest stellar populations as the primary ionizing source, although an AGN fraction of 14% is cannot be entirely ruled out. These galaxies are efficient ionizing photon producers, with mean log($\xi_{\rm ion}$ [Hz erg$^{-1}$])=25.37, exceeding typical values at similar redshifts. Escape fractions, however, are heterogeneous: 16% of EELGs at $z<7$ show escape fractions $>$5% for both Ly${\alpha}$ and LyC photons, while 82% lack detectable Ly${\alpha}$ emission. The median inferred LyC escape fraction is modest (5%) but enhanced in compact super-Eddington systems with sSFR >25 Gyr$^{-1}$. These results indicate that EELGs contribute approximately 16-40% of the total ionizing emissivity required to sustain hydrogen reionization. EELGs are extremely compact, with a median effective radius of 0.49 kpc, and exhibit a recent star-formation burst. Our analysis indicates that sSFR and star-formation rate surface density are the primary drivers of their extreme emission line strengths.

We present an analysis of the gas-phase mass-metallicity relationship (MZR) and fundamental metallicity relationship (FMR) for $65$ star-forming galaxies at $2 < z < 8$ from the JWST/EXCELS survey. We calculate gas-phase metallicities (12 + log(O/H)) using strong-line calibrations explicitly tested against the EXCELS sample, and report direct-method metallicities for $19$ galaxies. Our sample spans $8.1<\log(\rm M_\star/M_\odot)<10.3$ and $0<\log(\rm SFR/M_\odot \, yr^{-1})<2$, consistent with main-sequence star-forming galaxies at the same redshifts. We find a clear MZR at both $2<z<4$ ($\langle z \rangle = 3.2$) and $4<z<8$ ($\langle z \rangle = 5.5$), with consistent slopes and mild evolution in normalization of $\simeq 0.1 \, \mathrm{dex}$, matching trends from simulations and recent observations. Our results demonstrate rapid gas-phase enrichment in the early Universe; galaxies at $z \simeq 3$ (within the first $\simeq 15$ per cent of cosmic time) are enriched to $\simeq 40$ per cent of the metallicity of equivalent mass galaxies at $z=0$. We find tentative evidence for SFR-dependence in the MZR scatter, though results remain inconclusive and highlight the need for larger high-redshift samples. Comparison with locally derived FMRs reveals a clear offset consistent with other $z > 3$ studies. We discuss potential drivers of this offset, noting that high-redshift samples have significantly different physical properties compared to local samples used to define the $z=0$ FMR. Our results confirm that low-mass, high specific star-formation rate galaxies common at high redshift are inconsistent with the equilibrium conditions underlying the local FMR, and highlight the rapid chemical enrichment at early cosmic epochs.

The adsorption of volatile molecules onto dust grain surfaces fundamentally influences dust-related processes, including condensation of gas-phase molecules, dust coagulation, and planet formation in protoplanetary disks. Using advanced ab-initio density functional theory with r$^2$SCAN+rVV10 van der Waals functionals, we calculate adsorption energies of H$_2$, H$_2$O, and CO on carbonaceous (graphene, amorphous carbon) and silicate (MgSiO$_3$) surfaces. Results reveal fundamentally different adsorption mechanisms: weak physisorption on carbonaceous surfaces ($|\Delta\epsilon_{\rm ad}|\sim 0.1-0.2~{\rm eV}$) versus strong chemisorption on silicates ($|\Delta\epsilon_{\rm ad}|\sim 0.5-1.5~{\rm eV}$) via coordination bonds. Kinetic Monte Carlo simulations incorporating these energies demonstrate divergent surface evolution: carbonaceous grains exhibit distinct condensation radius compared to silicates, while the cocrystal of H$_2$O and CO significantly increases the desorption temperature of CO. The actual radii of gas-phase molecule depletion could thus be a comprehensive result of temperatures, chemical compositions, and even evolution tracks. Meanwhile, silicates maintain chemisorbed molecular coatings throughout most disk regions. Such dichotomy in surface coverage could also provide a natural mechanism for carbon depletion in inner planetary systems.

BL Lac has entered an active state since 2020, with multiwavelength observations revealing intense flares. In this study, we conducted 12-night multicolor optical monitoring using an 85 cm telescope from 2020 September to 2024 June and collected long-term broad-band archived data from radio to $\gamma$-rays. Intraday variabilities were detected on four nights, and most of them exhibited a bluer-when-brighter behavior. Both clockwise and counterclockwise spectral hysteresis loops were found within a single night. However, no reliable intraband time lag was detected for the intranight variabilities. On long timescales, the cross-correlation analysis shows that the variations of the optical, X-ray, and $\gamma$-ray bands do not reveal an obvious time delay, while the variations in the radio bands lagged them by about 370 days. The measured time lags suggest two distinct emission regions respectively responsible for the optical to $\gamma$-ray radiation and for the radio radiation, with a spatial separation of approximately $4.50\times10^{19}\ \rm cm$. We modeled the broad-band spectral energy distributions during four flaring epochs and one quiescent epoch, and found evidence for the possible persistent existence of a very high energy emission region. We also confirmed a spectral evolution of the source from an intermediate synchrotron peaked BL Lac object to a low synchrotron peaked BL Lac object.

In the wideband paradigm of pulsar timing, the time of arrival of a pulsar pulse is measured simultaneously with the corresponding dispersion measure from a frequency-resolved integrated pulse profile. We present a new method for performing wideband measurements that rigorously accounts for measurement noise. We demonstrate this method using observations of PSR J2124$-$3358 made as part of the Indian Pulsar Timing Array experiment using the upgraded Giant Metre-wave Radio Telescope, and show that our method produces more realistic measurement uncertainty estimates compared to the existing wideband measurement method.

Nuclear star clusters (NSC) are dense and compact stellar systems, of sizes of few parsecs, located at galactic centers. Their properties and formation mechanisms seem to be tightly linked to the evolution of the host galaxy, with potentially different formation channels for late- and early-type galaxies (respectively, LTGs and ETGs). While most observations target ETGs, here we focus on the NSC in M~74 (NGC~628), a relatively massive, gas-rich and star-forming spiral galaxy, part of the PHANGS survey. We analyzed the central arcmin of the PHANGS-MUSE mosaic, in which the NSC is not spatially resolved. We analyzed the NSC stellar populations in a point spread function (PSF) aperture, and compared it to the host galaxy. Within the PSF size, the NSC is contaminated by the host-galaxy light. We performed a two-dimensional spectro-photometric decomposition of the MUSE cube, employing a modified version of the C2D code, to disentangle the NSC from its host. This method provided different data cubes for the NSC and the host galaxy, allowing for both their comparison in a PSF aperture, and the spatially resolved analysis of the host. Our results show a very old and metal-poor NSC, in contrast to the surrounding regions. While similar properties were found in NSCs hosted by galaxies of different masses and/or morphological types from M~74, they are somewhat unexpected for a relatively massive star-forming spiral galaxy. The spatially resolved stellar populations of the host galaxy display much younger (light-weighted) ages and higher metallicities, especially in the central region (${\sim}500$~pc) surrounding the NSC. This suggests that this NSC formed a long time ago, and evolved passively until today, without any further growth. No significant amounts of gas would have reached the very central region in the last 8~Gyr.

We present 4,110 strong gravitational lens candidates, 3,887 of which are new discoveries, selected from a sample of 5,837,154 luminous red galaxies (LRGs) observed with the Dark Energy Spectroscopic Instrument (DESI). Candidates are identified via the presence of background ionized oxygen [O II] nebular emission lines in the foreground LRG spectra which may originate from the lensing of higher redshift star-forming galaxies. Using the measured foreground redshift, background redshift, and integrated flux of the background [O II] doublet, we integrate over impact parameters to compute the probability that each candidate is a lens. We expect 53% of candidates to be true lenses with Einstein radii ranging from 0.1'' to 4'', which can be confirmed with high-resolution imaging. Confirmed strong lenses from this sample will form a valuable cosmological dataset, as strong gravitational lensing is the only method to directly measure dark matter halo substructure at cosmological distances. We independently recover the host of the multiply imaged gravitationally lensed type Ia supernova iPTF16geu. Monitoring these lenses for future multiply lensed transients will enable (a) H0 measurements via time-delay cosmography and (b) substructure measurements via flux ratios.

The slitless spectroscopy mode of the NISP onboard Euclid has enabled efficient spectroscopy of objects within a large FoV. We present a large and homogeneous sample of bright quasars identified from the Euclid Quick Data Release (Q1) by combining high-purity candidate selections from Gaia and WISE with the NISP spectra. Through visual inspection of the Euclid spectra of these quasar candidates, we identify approximately 3500 quasars with reliable redshifts at $0<z\lesssim 4.8$. Of these, 2686 are new spectroscopic identifications relative to existing public compilations. We generate the first Euclid composite spectrum of quasars covering rest-frame NUV to NIR wavelengths without telluric lines, which will be pivotal to NIR quasar spectral analysis. We obtain an empirical spectroscopic depth of $J_{\rm E}\lesssim 21.5$ and $H_{\rm E}\lesssim 21.3$ at the sensitivity of the Wide Field Survey, beyond which the number of securely identified quasars declines sharply. We analyse VIS morphologies using Sersic and CAS metrics, and a deep-learning PSF fraction to track nuclear dominance. At low redshift ($z<0.5$), obvious host structures are common and a single Sersic model fits about half of the sources; at intermediate redshift ($0.5<z<2$), the nuclear component dominates, with 90% of the Sersic fits saturating at the upper index limit. In this intermediate redshift regime, $f_{\rm PSF}$ is available, and we use it as a more reliable compactness measure than the single-Sersic and CAS parameters to quantify nuclear versus host emission. We also explore the novel Euclid NIR colour space and discuss the role of these quasars in refining AGN selection techniques for future Euclid data releases. The spectroscopic bright quasar catalogue of this work, and the composite quasar spectrum, will be available at this https URL. (abridged)

$N$-body simulation serves as a critical method for modeling cosmic evolution and poses a significant challenge in high-performance computing. We present CUBE2, an open-source cosmological $N$-body code emphasizing memory efficiency, computational performance, scalability and precision. The core of its algorithm utilizes multi-level Particle-Mesh (PM) method to solve the Poisson equation for matter distribution, leveraging the well-optimized Fast Fourier Transform (FFT) for computational efficiency. Precision is ensured by the optimized Green's function that seamlessly bridges gravitational interactions between multi-level PM and Particle-Particle (PP) calculations. The program design enhances per-core/node efficiency in processing $N$-body particles, while the Information Optimized Storage (IOS) addresses memory constraints for large particle counts. Using CUBE2, we run two cosmological simulations with particle counts of $6144^3$ on the Advanced Computing East China Sub-center (ACECS) to test performance and accuracy.

The alignment between the radio-based International Celestial Reference Frame (ICRF) and the optical Gaia Celestial Reference Frame (Gaia-CRF) is critical for multi-waveband astronomy, yet systematic offsets at the optical bright end (G<13) limit their consistency. While radio stars offer a potential link between these frames, their utility has been restricted by the scarcity of precise Very Long Baseline Interferometry (VLBI) astrometry. In this study, we present new VLBI astrometry of 11 radio stars using the Very Long Baseline Array (VLBA), expanding the existing sample with positions, parallaxes, and proper motions measured. All 11 radio stars were detected, for 10 of which parallaxes and proper motions can be estimated, achieving median uncertainties better than 0.1 mas and 0.1 mas/yr, respectively. These new samples greatly contribute to the link between ICRF and Gaia-CRF at the optical bright end.

Warm and self-interactive dark matter cosmologies have been proposed as nonbaryonic solutions to the tensions between the $\Lambda$ cold dark matter model and observations at the kiloparsec scale. In this paper, we used the dark matter-only runs of the \textsc{aida-tng} project, a set of cosmological simulations of different sizes and resolutions, to analyze the macroscopic impact of alternative dark matter models on the abundance, radial distribution, and clustering properties of halos. We adopted the halo occupation distribution formalism to characterize the evolution of its parameters $M_1$ and $\alpha$ with the mass and redshift selection of our sample. By dividing the halo population into centrals and satellites, we were able to study their spatial density profile. We found that a Navarro-Frenk-White model is not accurate enough to describe the radial distribution of subhalos and that a generalized Navarro-Frenk-White model is required instead. Warm dark matter models, in particular, present a cuspier distribution of satellites, whereas self-interacting dark matter exhibits a shallower density profile. Moreover, we found that the small-scale clustering of dark matter halos provides a powerful tool for distinguishing among alternative dark matter scenarios, in preparation for a more detailed study that fully incorporates baryonic effects and for a comparison with observational data from galaxy clustering.

We present upgraded Giant Metrewave Radio Telescope (uGMRT) observations of 0731+438, an \mbox{FR II} radio galaxy at a redshift of 2.429 with two lobes separated by 82 kpc. A blueshifted, faint and broad \mbox{H{\sc i}} 21 cm absorption line with velocity full width at half maximum (FWHM) $\sim 600\,\rm km\,s^{-1}$ is detected against the southern radio lobe that is 47 kpc from radio core, indicating a neutral hydrogen outflow associated with jet-gas interaction beyond the host galaxy. The outflow has a mass outflow rate of $\sim\,0.4T_{\rm s}\Omega\rm\, M_\odot\,{\rm yr}^{-1}$, which could increase to $\sim\,4.0T_{\rm s}\Omega\rm\,M_\odot\,{\rm yr}^{-1}$, corresponding to an energy outflow rate of $2.4T_{\rm s}\Omega\times10^{40}$ -- $1.5T_{\rm s}\Omega\times10^{41}\,\rm erg\,s^{-1}$, where $T_{\rm s}$ is the spin temperature and $\Omega$ is the solid angle of the outflow. Previous optical observations identified an extended emission line region aligned with the radio axis, ionized by the central Active Galactic Nucleus (AGN). Within this region, a warm and ionized outflow with a mass outflow rate of $\sim\,50\rm\, M_\odot\,{\rm yr}^{-1}$ and an energy outflow rate of $\sim1.7\times10^{43}\,\rm erg\,s^{-1}$ was detected. We propose that both the extended emission line region and the optical outflow are results of synergistic effect between jet and AGN radiation. The AGN likely exerts negative feedback on the host galaxy, as evidenced by the gas expulsion by the jet and the high velocity dispersion of ionized gas observed optically. So far, detections of jet-driven neutral hydrogen outflows remain rare. The high redshift, large outflow radii, substantial mass outflow rate and energy outflow rate of the neutral hydrogen outflow in 0731+438 expand the known parameter space of such outflows.

The dominant radiation mechanism that powers the prompt $\gamma$-ray emission in gamma-ray bursts (GRBs) remains poorly understood. High quality, time- and energy-resolved linear polarization measurements of prompt $\gamma$-ray photons can distinguish between synchrotron and inverse-Compton processes and provide crucial constraints on the outflow properties. This will be achieved by POLAR-2 that is proposed as a dedicated GRB polarimeter and successor to POLAR. The High-energy Polarimetry Detector (HPD) is one of the three instruments of POLAR-2 that features significantly improved sensitivity in the $(40-1000)$\,keV energy range and a detection area four times larger than that of POLAR. Here we demonstrate the capabilities of the HPD to constrain key physical model parameters by creating and fitting to synthetic sources using a time-resolved spectro-polarimetric theoretical model of prompt GRB emission. The time-resolved spectral and polarization fits are performed using a novel technique featuring maximum likelihood over an unbinned (in time and energy) list of detected events. The constrained model parameters directly relate to the underlying source physics that would reveal an accelerating, coasting or decelerating emission region. For a pulse fluence of $\mathcal{F}=10^{-5}\mathcal{F}_{-5}\,{\rm erg\,cm^{-2}}$ and higher we can constrain the time-integrated polarization degree to an absolute accuracy ($1\,\sigma$) of about $2.2\mathcal{F}_{-5}^{\,-1/2}$ per cent, as long as source photons dominate over the background. In bright GRBs, such unprecedented accuracy at these energies will allow to distinguish between different models for the prompt GRB emission mechanism and constrain the magnetic field geometry, jet angular structure and outflow composition.

The results of a study of the cosmic ray composition in individual events in the energy range (2-12.5) EeV using the muon correlation method is presented. The considered sample included showers with zenith angles less than 60 degrees recorded in the period 1974-2018. The existence of four separate groups of primary particles with different origins is confirmed. The obtained results have potential importance for understanding the composition of cosmic rays in the specified primary energy range.

It is still not clear which environmental processes operate in filaments. Given the ubiquity of filaments and their importance in feeding clusters, a proper understanding of these mechanisms is crucial to a more complete picture of galaxy evolution. To investigate them, we need large galaxy samples with spatially resolved information. As part of this effort, we analyse resolved H$\alpha$ maps of 685 galaxies inside and outside the filaments around the Virgo cluster in addition to extensive measurements of integrated physical properties. We create a pipeline to decompose the H$\alpha$ images into individual clumps that trace star forming regions. We find that the number and average size of clumps in a galaxy are well-defined functions of distance and angular resolution. In particular, the power-law relation between the number of clumps and the distance of a galaxy is consistent with a fractal structure of star forming regions. We formulate an algorithm to compare filament and non-filament galaxies after removing observational differences. Although we do not have any conclusive evidence for a difference in clump size distributions between filament and non-filament galaxies, we do find that filament galaxies have slightly more peripheral clumps than their non-filament counterparts.

The fractal dimension $D$ is used to map the large-scale galaxy distribution in the Universe by color types: blue, green and red. Using a $NUVrK$-complete COSMOS2020 subsample of 618,952 galaxies observed up to $z=4$, number densities were derived and plotted against two cosmological distance measures, the luminosity and comoving (galaxy area) distances, in order to estimate $D$ for each galaxy color type in two redshift intervals: $z\gtrless1$. We found a general gradient $D_{\mathrm{blue}}> D_{\mathrm{red}}>D_{\mathrm{green}}$ with $D=1.40-2.03$ for $z<1$. For $1<z\leq4$, the gradient changes to $D_{\mathrm{blue}}>D_{\mathrm{green}}>D_{\mathrm{red}}$, and the fractal dimension values are lower, $D=0.03-0.44$. These results suggest that the fractal dimension is a sensitive diagnostic for how galaxy populations trace the evolving cosmic web, and confirm the fractal dimension as a useful tool for observational mapping of large-scale structure by galaxy color.

We present a multi-wavelength study of nuclear outflows in the nearby dwarf Seyfert galaxy NGC~4395, which hosts an intermediate-mass black hole. Using \textit{JWST}/NIRSpec and MIRI IFU spectroscopy (1.66--28.6~$\mu$m), together with ALMA and Gemini/GMOS data, we probe the ionised and molecular gas on parsec scales. The JWST nuclear spectra reveal 134 emission lines, including H\,\textsc{i}, He, numerous fine-structure lines, H$_2$ rotational/ro-vibrational transitions, and several PAH bands. Modelling of the H$_2$ rotational lines reveals three warm/hot molecular components ($T\!\approx\!580$, 1480, and 2900~K), along with a cold ($<50$~K) phase traced by ALMA CO(2--1). Outflow signatures are detected in cold and warm/hot molecular gas, in H\,\textsc{i}, and in 36 fine-structure lines spanning ionisation potentials of 7.6--300~eV. Ionised outflow velocities range from 127 to 716~km\,s$^{-1}$, with blueshifted and redshifted components consistent with a stratified biconical geometry. The cold molecular gas shows a mass outflow rate nearly 1--2 orders of magnitude larger than that of the warm/hot molecular and ionised phases. The kinetic coupling efficiency is 0.003--0.12\% for the coronal-line gas and 0.4--1.4\% for the H\,\textsc{i} outflow, indicating that only the low-ionisation gas significantly impacts the surrounding ISM. Outflow velocity and the fraction of flux in the outflowing component increase with ionisation potential, implying that the most highly ionised gas originates closest to the AGN and is most efficiently accelerated.

We analyse CMB constraints on stable networks of cosmic strings and domain walls using for the first time full Planck 2018 data together with BICEP/Keck 2018 $B$-mode measurements. The defect-induced anisotropies are computed using the Unconnected Segment Model for Nambu-Goto and Abelian-Higgs strings, as well as for stable domain walls, and included in a full Markov Chain Monte Carlo analysis jointly varying all $\Lambda$CDM parameters, the tensor-to-scalar ratio, and the string/domain wall tension. No statistically significant evidence for defects is found, although we observe a mild preference for non-zero cosmic string tension. Our results improve previous constraints on the defect power spectrum by up to a factor of two. In the particular case of strings, the improvement is driven by the $B$-mode data, and is especially pronounced for Abelian-Higgs strings. We also present forecasts for the Simons Observatory and the LiteBIRD satellite: the former will improve string tension constraints by about a factor of three, while the latter those on the domain wall tension by about a factor of ten. Finally, we assess the impact of Nambu-Goto string loops on CMB anisotropies in light of both current and future observations.

We present data from the ALMA Central Molecular Zone Exploration Survey (ACES) Large Program, which provides broad spectral-line and 3 mm continuum coverage of the Central Molecular Zone (CMZ) at a spatial resolution of 0.1 pc. The survey delivers homogeneous, wide-field mosaics that enable direct comparisons of the physical and chemical conditions across diverse environments in the Galactic center. In this data release paper, we present the CS(2-1), SO(2_3-1_2), CH3CHO(5_1,4-4_1,3), HC3N(11-10), and H40a lines observed simultaneously within two broad spectral windows. These lines reveal pronounced spatial and chemical variations across the CMZ, tracing distinct components of molecular gas, shock-affected regions, and ionized structures. The high angular resolution and multi-line capability of the ACES dataset make it a powerful resource for future studies of gas dynamics, star formation activity, and the physical connection between the CMZ and Sgr A*.

Dedispersion is the computational process of correcting for the frequency-dependent time delay affecting a radio signal that propagates through the interstellar and intergalactic media. It is a crucial component of transient search pipelines that maximises the signal-to-noise ratio, especially when targeting highly dispersed signals: for instance, pulsar emissions making their way through a dense cloud of ionised gas, and fast radio bursts travelling cosmological distances. This paper introduces Streaming high Time-Resolution Imaging DEdispersion (STRIDE), a novel dedispersion algorithm to generate per-pixel dedispersed time series from high time and frequency resolution interferometric images. Unlike straightforward approaches to image dedispersion, STRIDE does not involve expensive manipulation of the input data layout, such as explicitly building dynamic spectra or shifting images. Furthermore, it is the first dedispersion algorithm to partition a dispersive sweep over the time dimension, in addition to frequency. As a consequence, images corresponding to the entire time span of the target dispersive delay are not required all at once. Instead, the algorithm works with an arbitrarily-sized subset of images at a time, adopting an incremental, streaming-based approach to dedispersion. In evaluating STRIDE on the presented test case, it is shown that the minimum memory requirement is reduced by 97.9%, going from 684.5 GB to 14.4 GB. As current and future generations of widefield interferometers increasingly turn to imaging techniques for detection and localisation of radio transients, STRIDE positions itself as a strong alternative to traditional dedispersion methodologies. It arguably is the only viable option for imaging-based searches with low-frequency instruments such as the Murchison Widefield Array (MWA) and low-frequency Square Kilometre Array (SKA-Low).

Dark matter may accumulate in neutron stars given its gravitational interaction and abundance. We investigate the influence of strongly-interacting dark matter, described by a QCD-like one-flavor $G_2$ gauge theory, on neutron stars. This choice allows to test, for the first time, a first-principles-determined non-Abelian dark matter equation of state, which supports composite fermionic dark matter and thus a Fermi-pressure-stabilized dark matter component. The ordinary matter part of the mixed star is described by available model-agnostic equations of state that interpolate between the low-density regime and high-density regime. We find that strongly-interacting dark matter has a similar impact on neutron stars as other model equation of states and confirm that strongly-interacting dark matter can be accommodated by constraints from neutron star observations within our uncertainties.

In the single-field case, Horndeski provides the most general scalar-tensor theory with second-order field equations. By contrast, systematic multi-field extensions remain incomplete: while the general field equations for the bi-Horndeski case are known, a general action has not been established, and for cases with three or more fields, neither a general action nor general equations are available. We characterize Horndeski by two mild axioms: closure under invertible pure disformal transformations and the requirement that the theory includes the minimal Horndeski theory. Under this characterization, we recover the standard single-field action up to boundary terms and obtain a practical path to multi-field constructions. In particular, we show that antisymmetric structures, such as those identified by Allys, Akama, and Kobayashi, appear within this framework, and indicate that this viewpoint has the potential to account for features captured by known bi-Horndeski equations.

We experimentally demonstrate a novel interferometric architecture for next-generation gravity missions, featuring a laser ranging interferometer (LRI) that enables monoaxial transmission and reception of laser beams between two optical benches with a heterodyne frequency of 7.3 MHz. Active beam steering loops, utilizing differential wavefront sensing (DWS) signals, ensure co-alignment between the receiving (RX) beam and the transmitting (TX) beam. With spacecraft attitude jitter simulated by hexapod-driven rotations, the interferometric link achieves a pointing stability below 10 urad/$\mathrm{\sqrt{Hz}}$ in the frequency range between 0.2 mHz and 0.5 Hz, and the fluctuation of the TX beam's polarization state results in a reduction of 0.14\% in the carrier-to-noise-density ratio over a 15-hour continuous measurement. Additionally, tilt-to-length (TTL) coupling is experimentally investigated using the periodic scanning of the hexapod. Experimental results show that the on-axis LRI enables the inter-spacecraft ranging measurements with nanometer accuracy, making it a potential candidate for future GRACE-like missions.

We discuss graviton dark radiation produced by the decay of a secluded scalar field that couples to the Standard Model (SM) only through gravity. Such scalar fields are long-lived, and their decay channels generically include gravitons. If such particles existed and dominated the early universe, a sizable branching ratio into gravitons would yield non-negligible dark radiation that significantly alters the subsequent thermal history of the universe. In this work, we focus on the dark glueball as a representative secluded hidden scalar and compare the decay rates into SM particles via a non-minimal coupling to gravity with those into gravitons, paying attention to how the breaking of conformal invariance affects the amount of graviton dark radiation. We find that decays into the SM are dominated by two-body decay channels into Higgs bosons and gluons. In particular, when the Higgs field has a large non-minimal coupling to gravity, the production of graviton dark radiation is naturally suppressed in the metric formalism, and the SM sector is preferentially reheated and energy transfer to other hidden sectors is suppressed. Finally, we present the expected gravitational-wave spectrum resulting from dark glueball domination.

We investigate the shadow cast by a rotating traversable wormhole in the Teo class endowed with a general redshift function, with particular emphasis on the emergence of cuspy structures. The shadow boundary is the common envelope of two critical orbit families: unstable circular orbits outside the throat and orbits at the throat itself. The formation of cusps, marking the transition between smooth and cuspy shadow boundaries, only becomes possible when the redshift parameter $\lambda$ is allowed to vary. Moreover, we uncover a universal critical value $\lambda_c$ that signals the onset of the cusp. A phase diagram characterized by the spin and redshift parameters reveals four distinct morphologies: smooth, cuspy, ears touching, and throat drowning. The morphology of the wormhole shadow may provide observational diagnostics for the different compact objects in future high-resolution imaging observations.

The Sun is a promising target for dark matter (DM) searches due to its ability to accumulate DM particles via scattering and catalyze their self-annihilation. However, at low DM masses, DM particles can also "evaporate" due to subsequent collisions with the hot thermal plasma of the Sun. While several modeling studies have calculated the competitive dynamics of DM evaporation and annihilation, observational studies have typically assumed a fixed 4 GeV "evaporation limit", below which DM evaporates before it can annihilate. In this paper, we consider the competitive effects of DM evaporation and annihilation on spin-dependent DM nucleon cross-section limits, finding that Solar observations can continue to exceed terrestrial constraints by between 1-5 orders of magnitude for DM masses between 2-4 GeV, and can even provide world leading constraints below 0.2 GeV where direct detection is limited.

We present theoretical and experimental evidence for a new phase of matter in neutron-rich systems that we call multimodal superfluidity. Using ab initio lattice calculations, we show that the condensate consists of coexisting s-wave pairs, p-wave pairs in entangled double pair combinations, and quartets composed of bound states of two s-wave pairs. We identify multimodal superfluidity as a general feature of single-flavor spin-1/2 fermionic systems with attractive s-wave and p-wave interactions, provided the system is stable against collapse into a dense droplet. Beyond neutrons at sub-saturation densities, we demonstrate that this phase appears in generalized attractive extended Hubbard models in one, two, and three dimensions. We elucidate the mechanism for this coexistence using self-consistent few-body Cooper models and compare with Bardeen-Cooper-Schrieffer theory. We also derive the form of the effective action and show that spin, rotational, and parity symmetries remain unbroken. Finally, we analyze experimental data to show that p-wave pair gaps and quartet gaps are present in atomic nuclei, and we discuss the consequences of this new phase for the structure and dynamics of neutron star crusts.